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Abstract 

The isotropy and homogeneity of the cosmic microwave background (CMB) favors 
"scalar driven" early Universe inflationary models. However, gauge fields and other 
non-scalar fields are far more common at all energy scales, in particular at high energies 
seemingly relevant to inflation models. Hence, in this review we consider the role 
and consequences, theoretical and observational, that gauge fields can have during 
inflationary era. Gauge fields may be turned on in the background during inflation, or 
may become relevant at the level of cosmic perturbations. There have been two main 
class of models with gauge fields in the background, models which show violation of 
cosmic no-hair theorem and those which lead to isotropic FLRW cosmology, respecting 
the cosmic no-hair theorem. Models in which gauge fields are only turned on at the 
cosmic perturbation level, may source primordial magnetic fields. We also review 
specific observational features of these models on the CMB and/or the primordial 
cosmic magnetic fields. Our discussions will be mainly focused on the inflation period, 
with only a brief discussion on the post inflationary (p)reheating era. 
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1 Introduction 



The prospects of research in cosmology has been in an ever accelerating expansion since the 
discovery of Hubble expansion in late 1920's. The cosmological and astrophysical observa- 
tions since then has been in favor of the Big Bang model which has been used to explain the 
presence of Cosmic Microwave Background (CMB), cosmic abundances of light nuclei and to 
reconstruct the thermal history of the Universe especially after the Big Bang Nucleosynthesis 
(BBN) when the temperature of the Universe was around 1 MeV. Gaining a detailed knowl- 
edge and picture of the history of Universe before that is limited since the main messengers 
from the early University are the electromagnetic waves reaching to us, and that they can 
only (almost) freely propagate after the recombination time, when the CMB released. 

Despite of the successes, some basic questions about the early Universe remain unan- 
swered within the (hot) Big Bang model. These questions which are horizon problem, flatness 
problem, relic problem, the source for CMB anisotropy (CMB temperature fluctuations) and 
seeds for the large structures, can flnd a suitable and natural answer once the hot Big Bang 
model is augmented with an inflationary era, a period of accelerated expansion [U [2j. This 
constitutes the Standard Model of cosmology. In this review we will focus on the inflationary 
period. 

Many different models of inflation have been proposed and studied in the literature. In 
these models inflation is generically driven by the coupling of one or more scalar flelds to 
gravity, and the dynamics during inflation is such that generically the potential energy of the 
flelds dominate over their kinetic term and the potential is flat enough to ensure the so-called 
slow- roll inflation (to be deflned in section [2]). The current observations indicate that the 
Hubble parameter during inflation H has an upper bound, H < 10~^Mpi, corresponding to 
inflaton energy density of order p^f < lO^^^Mpj ~ (lO^^GeV)"^. (Mpi is the reduced Planck 
mass. Our conventions and notations are summarized in Appendix [Al) The inflaton energy 
density is at least two orders of magnitude smaller than Mpi and of order of the energy scale 
of the Grand Unifled Theories (GUT's) of particle physics. One would hence naturally expect 
that inflation should be formulated within the proposed particle physics models working in 
the same energy range. 
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Models of particle physics are (chiral) gauge field theories the matter content of which 
generically include gauge fields, (chiral) fermions and scalar (Higgs-type) fields, which are in 
various representations of the gauge group. The scalar field(s) of these models are used for 
spontaneous symmetry breaking (SSB) and hence giving masses to chiral fermions through 
Yukawa couplings. Therefore, their potential is generically tuned for the sjTumetry breaking 
purposes. This latter condition is, however, at odds with the slow-roll inflation requirement 
demanding a lower scalar masses. In other words, generically, the same potential used for 
the Higgs-type fields of (beyond) particle physics standard model is not fit for infiationj^ 

The above general features would hence prompt the idea of using fields other than scalars 
in particle physics models for inflationary model building. Turning on non-scalar fields 
during inflation appears at first to be incompatible with isometry and rotational symmetry 
of the observed Universe; anisotropy at cosmological scales, if exists, should remain small 
enough to be compatible with the observations. This article is partly devoted to reviewing 
different ways proposed in the literature so far, in which gauge fields can be turned on in the 
inflationary background such that the observational anisotropy constraints are also satisfied. 
We will also review the observational signatures and features of these models. 

Effects of gauge fields may arise at the level of perturbations, without being turned on 
at the background level, or gauge fields may be turned on at the background level, as well 
as arising through perturbations. The former group, is what we first analyze in section [3] 
and the latter is what we study in the rest of the review. In the latter group of models with 
vector gauge fields in the background, one can distinguish two general classesH The first class 
includes models where contribution of the gauge field to the energy budget during inflation 
is small, i.e. although gauge fields are turned on in the inflationary background, inflation is 
mainly driven by other scalar fields. These models can have a controUed/controlable, but yet 
non-exponentially suppressed, anisotropy at both background and perturbation levels. The 
other class mainly consists of gauge-flation and chromo-natural inflation models where the 
contribution of the gauge field sector has the dominant contribution to the energy budget of 
the Universe during inflation, while the background geometry is isotropic. 

This review is organized as follows. Section [2] is a very quick review of basics of in- 
flationary cosmology. In this section we also fix our conventions and notations, as well as 
summarizing the current observational data which is used to constrain inflationary models. 
We also mention CMB anomaly which motivated people to study the statistical anisotropy 
or the effect of gauge fields on inflation. In section [31 we analyze models with gauge fields 
at perturbation level, which can be used to provide seeds for primordial magnetic fields 

^There are, however, ways to get around this general argument in specific models by relaxing the as- 
sumptions this result is based on. For example the above argument was based on the implicit assumption 
that gravity at the inflation scale is described by the Einstein GR and that Higgs-type fields are minimally 
coupled to gravity [3]. Another way to get around this argument is to explore the supersymmetric models, 
which have a huge wealth of scalar fields as superpartners of standard model fermions, and more options 
for fulfilling the flatness of the potential, e.g. see [4]. In this review our discussions will mainly be within 
Einstein GR with minimally coupled flelds. 

^In this review we will not consider the vector inflation models [S]. Having a vector field with standard 
kinetic term, without the gauge symmetry, i.e. when the vector field is not a gauge field, we expect to have 
ghost instability [6] (unless the gauge symmetry is broken via SSB mechanism) and hence not generically a 
theoretically viable framework. 
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and may also source non-Gaussianity (NG) on the CMB power spectrum and primordial 
gravitational waves. In section HJ we discuss anisotropic inflationary models, with gauge 
fields in the background. In these models contribution of the gauge field to the infiationary 
energy budget remains small. Consequently, the anisotropy in the expansion is also small. 
Nonetheless, observational effects of gauge fields may not be negligible. We explain how 
anisotropic interaction induces interesting phenomenology. We also discuss the possibility 
of testing anisotropic inflation. In section |5l we review the gauge-flation model, a model 
in which inflation is driven by a non-Abelian gauge theory minimally coupled to Einstein 
GR. We review and study inflationary background and study the stability of the inflationary 
trajectories in this model. We also discuss cosmic perturbation theory in gauge-flation and 
confront the model with the observational data. In section [6], we analyze an extension of 
the gauge-flation model, where besides the non-Abelian gauge flelds we have an axion fleld, 
"chromo-natural inflation". We present both background and perturbation theory of the 
chromo-natural model. The last section is devoted to discussion and outlook. To make the 
review self-contained in some appendices we have gathered discussions and analysis which 
are related to the topic of this article. Appendix |X] contains our conventions. In Appendix 
[B] we present a quick review of 6N formulation for cosmic perturbation theory analysis. 
Appendix [C] contains a review of Bianchi cosmologies relevant to the anisotropic models. 
In Appendix [Dl we present Wald's cosmic no-hair theorem [7j and the way it should be 
extended in inflationary cosmology [8]. 

2 Preliminaries of infiationary cosmology 

Inflationary paradigm has secured its status as the leading candidate with the currently 
available cosmic data. Moreover, inflation has the appealing feature that it can be formulated 
within the standard existing Einstein GR and quantum fleld theory frameworks. In this 
section we will briefly discuss basics of FLRW cosmology and review some simple models of 
inflation, through this we also flx the notations we will be using throughout the article. To 
this end, we flrst review background (classical) slow-roll inflationary trajectories, and then 
review cosmic perturbation theory. We also present a summary of the current observation 
data, coming from combined CMB and related astrophysical data. 

Our discussions in this section will be very brief. For more detailed discussion the reader 
is encouraged to consult the books and reviews on inflation, an incomplete list include [2]. 
Throughout this note we use natural units where the reduced Planck mass Mpi is 1. 

2.1 Infiationary backgrounds 

The standard cosmology is based on the Einstein general relativity which states that the 
dynamical evolution of our universe is described by a 4-dimensional geometry, g^^, governed 
by the Einstein equations, as 
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Here i?^^ is the Ricci tensor, R = R^^g^'^ is the Ricci scalar and T^j, is the energy-momentum 
tensor of the matter in the Universe. Furthermore, modern cosmology starts with the as- 
sumption of cosmological principle [2]: Universe is homogeneous and isotropic on cosmolog- 
ical large distances, 100 Mpc and larger. Thus, the line element of the space-time is given 
by the Friedman-Lemaitre-Robertson- Walker (FLRW) metric 

ds"" = -de + a\t) ( + r'(de' + sin^ Odcj)'')] , (2.2) 
\1 — Kr'^ J 

where t is the cosmic time, a{t) is the scale factor and K is the curvature constant and takes 
and —1 values respectively corresponding to close, fiat and open Universes. Due to 
the spatial isotropy and homogeneity of the metric (12.21) , the most general form of T^,y is a 
perfect fluid 

T^u = {p + P)u^u^ + Pg^,,, (2.3) 

where p and P are the energy and pressure densities. The Einstein equations for the isotropic 
homogeneous FLRW background reduces to the Friedman and Raychaudhuri equations, 
respectively 

H' = ]-p-^, (2.4a) 

O (1 

- = -l(p + 3P), (2.4b) 
a 

where a dot represents derivative with respect to cosmic time and H = ^ is the Hubble 
rate of the expansion. Ordinary matter (e.g. dust, radiation) have positive energy density 
and non-negative pressure and satisfy strong energy conditiorl^l (p + 3P > 0) which makes 
the RHS of equation fl2.4bl) positive. Thus, in conventional Big Bang theory which assumes 
p + 3P > 0, the Universe always exhibits a decelerated expansion (a < 0) [2]. 

Standard Big Bang theory, however, is not capable of explaining the current state of the 
Universe as we observe it today. Some of the most important problems with observations 
are the flatness problem (why our Universe is so nearly spatially fiat), the horizon problem 
(why the temperature of the CMB on the whole sky is so accurately the same), the monopole 
or heavy relic problem. Inflation, an accelerated expansion, a > 0, phase was proposed to 
overcome these problems [T]. For instance, inflation offers an explanation of flatness problem 
without any fine-tuning, since due to inflation the spatial curvature term, ^ in f l2.4ap . will be 
negligible. Furthermore, the fact that our entire observable Universe might have arisen from 
a single early causal patch can solve the horizon problem. The inflationary paradigm not 
only solves the above cosmological puzzles, but more importantly, also offers an explanation 
for the origin of large scale structures as well as CMB temperature anisotropy. See [2] for 
more detailed discussions. The simplest and most extensively studied inflation scenario is 
a quasi de Sitter inflation, almost exponential expansion, where equation of state of cosmic 
fluid is P ~ — p, H is almost a constant, and the Friedman equations has the solution 

a(t) ~ a(0)exp(i7t) . (2.5) 



^Strong energy condition states that for all time-like t'^, we have (T^y — \g^iuT)P't" > 0. 
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Based on the matter which drives accelerated expansion, inflation may have many re- 
alizations. The simplest inflationary model is a single scalar fleld, called inflaton fleld, 
minimally coupled to Einstein gravity 

S = l d'xV^i^ - lid.^r - V{^)). (2.6) 

Neglecting spatial gradients, the energy and pressure density are given by 

P=l^' + V, P=^-^'-V. (2.7) 

Also, the Friedman equation and the equation of motion of the scalar fleld are respectively 

3H' = ^cp' + Viip), (2.8) 
y3 + 3Hip + V^ = 0, (2.9) 

where = Note that since the spatial curvature term, ^ in f l2.4ap . damps exponentially 
during the inflation, we dropped it. As we can see from (12. 7p . in the limit that 0^ <^ V, the 
the equation of state satisfles p ^ —P, and an inflationary phase is possible with 

^ ^V{^). (2.10) 

In fact, in this simple model inflation can only occur if the scalar fleld satisfles the slow-roll 
conditions. Slow-roll inflation is quantifled in terms of the slow-roll parameters, e and r] 

and requires e, |?7| ^ 1. Smallness of e and r] is usually needed to make sure that inflation 
lasted long enough to solve the horizon and flatness problems. A useful quantity to measure 
the amount of inflation is the number of e-folds A^^^ 

Ar^ = ln(^) = / H{t)dt, (2.12) 



ao Jto 

where to cind tf denote the initial and the flnal time of inflation respectively, and Oq is 
the value of the scale factor at the beginning of the inflation. As an often-stated rule of 
thumb, in order to solve the horizon and flatness problems, we need about 60 e-folds [2]. 
Thus, Ne = 60 is taken as a standard minimum number of e-folds for inflationary models. 
However, the precise value of the required Ne depends on the energy scale of the inflation 
and the details of the reheating era [2]. Furthermore, perturbations (temperature anisotropy) 
observed in the CMB have left the horizon Nqmb — 40 — 60 e-folds before the end of inflation. 

Slow-roll parameters deflned in (12. lip are based on the time derivative of the Hubble 
parameter H during inflation. For speciflc models of inflation, however, it is often more 
useful to deflne the slow-roll parameters directly in terms of the inflaton fleld Lagrangian. 
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For the simple single scalar field model given in (12. 6p . one can easily show that slow- roll 
parameters (12. lip can also be written as 



e = -^, v = — —■ (2.13) 
2H^' ' Hip ^ ' 



It is also common to express the slow-roll parameters in terms of the potential V{lp) as 
Finally, infiation ends when e{ipend) = 1 and number of e-folds Ne is given as 

Ne^ / TT^y^^ / -T^dif. (2.15) 

-''Pund f J Vend V ^6 

Depending on the form of the potential, and also possibly the kinetic term which can 
be non-canonical in general, we may have different single-field models. A useful, but not 
exhaustive, classification of single-field inflationary models is as follows. 



• Large field models: The initial value of the infiaton field is large, generically super- 
Planckian, and it rolls slowly down toward the potential minimum at smaller ip values. 
For instance, chaotic infiation is one of the representative models of this class. The 
typical potential of large-field models has a monomial form as 

V{^) = \/o<^". (2.16) 

A simple analysis using the dynamical equations reveals that for number of e-folds Ne 
larger than 60, we require super-Planckian initial field valued, (p>Q > 3Mpi. For these 
models typically e rj 

• Small field models: Infiaton field is initially small and slowly evolves toward the po- 
tential minimum at larger cp values. The small field models are characterized by the 
following potential 

V{^) = Vo{l-i^r), (2.17) 
/i 

which corresponds to a Taylor expansion about the origin, but more realistic small 
field models also have a potential minimum at (y? 7^ which the system falls in at the 
end of infiation. A typical property of small field models is that a sufficient number of 
e-folds, requires a sub-Planckian infiaton initial value. For this reason they are called 
small field models. Natural infiation is an example of this type |12j . 

"'Note that e and 77 are called the Hubble slow-roll parameters while ey and ryy arc called potential 
slow-roll parameters. During the slow-roll these parameters are related as e ~ and rjv = e + rj. 

^In the presence of another natural cutoff A in the model, smallness or largeness of the infiaton field 
should be compared to A; A could be sub-Planckian and in general A < Mpi. For a discussion on this see 

mm- 
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Hybrid inflation models: These models involve more than one scalar field while infiation 
is mainly driven by a single infiaton field 0. Infiaton starts from a large value rolling 
down until it reaches a bifurcation point, (p = (pe, after which the field becomes unstable 
and undergoes a waterfall transition towards its global minimum. Its prime example 
is the Linde's hybrid infiation model with the following potential [I3] 



V{<P, X) = ^ix' - + + lm'<P'. (2.18) 

During the initial infiationary phase the potential of the hybrid infiation is effectively 
described by a single field while infiation ends by a phase transition triggered by 
the presence of the second scalar field, the waterfall field x- other words, when the 
effective mass squared of a waterfall field becomes negative, the tachyonic instability 
makes waterfall field roll down toward the true vacuum state and the infiation suddenly 
ends. 



Number of e-folds A^e is given as 



M 

a 

instability, x = becomes unstable and gets negative. 



where 0e = y is the critical value of the infiaton below which, due to tachyonic 



K-inflation: This is the prime example of models with non-canonical Kinetic term we 
discuss here. They are described by the action [1] 



S = J dV^(| + P(y,,X)), (2.20) 

where ip is a scalar field and X := — |(9^y9)^. Here, P plays the rule of the effective 
pressure, while the energy density is given by 

p = 2XP^x - P- (2.21) 
Thus, the slow-roll parameter is given as 

^XP^x 
^ 2XP^x-P' 

The characteristic feature of these models is that in general they have a non-trivial 
sound speed cj. for the propagation of perturbations (c/. our discussion in section [2^ 

= — . (2.22) 

P,x + 2XP,xx ^ ' 

Finding K-inflation actions P{ip,X) which are well-motivated and consistently embed- 
ded in high-energy theories is the main challenge of this class of models [9] . Nonethe- 
less, DBI inflation is a special kind of K-inflation, which is well-motivated from string 
theory with the action [15] 

S = l rf'xv/^[| - ((v^ - 1) + V(/)) ] , (2.23) 

where D = 1 - 2f{^)X. 
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2.2 Cosmic perturbation theory 

The FLRW Universe is just an approximation to the Universe we see, as it ignores all the 
structure and other observed anisotropies e.g. in the CMB temperature. One of the great 
achievements of inflation is having a naturally embedded mechanism to account for these 
anisotropies. The idea is based on the fact that inflaton(s) and metric are indeed quantum 
flelds and hence have quantum fluctuations. These quantum fluctuations are, however, 
"virtual" in the sense that they do not carry energy. Nonetheless, the cosmic expansion and 
existence of (cosmological) horizon makes it possible for some of these quantum fluctuations 
become "real" and classical. 

Quantum fluctuations appear in all wavelengths and comoving momentum numbers k. 
For the suhhorizon modes k ^ aH these fluctuations essentially appear as they are in 
flat spacetime, while for superhorizon modes k < aH fluctuations do not show oscillatory 
behavior, they freeze and behave as classical perturbations on the FLRW background. Given 
that k is constant and scale factor a{t) expands (almost exponentially), if we wait long 
enough any mode becomes classical. Nevertheless, inflation ends and not all the modes 
have had time to cross the horizon during inflationj^ Moreover, termination of inflation 
also opens the crucial possibility that some of the modes which have become superhorizon 
modes during inflation to become suhhorizon again in a later time after inflation ended and 
reenter the horizon. In particular, some of these modes, which are the modes left the horizon 
during 40 — 60 e-folds before the end of inflation, have reentered our cosmological horizon 
after the surface of last scattering and have been imprinted on the CMB. These classical 
modes are responsible for both structure formation and the CMB anisotropy. The precision 
measurements on the CMB anisotropies can then be used to restrict inflationary models. In 
this section, we review cosmic perturbation theory for inflationary models and extract data 
from the models which can be compared with the CMB observations to be reviewed in the 
next subsection. 

The standard cosmic perturbation theory starts with the assumption that during most of 
the history of the Universe deviation from homogeneity and isotropy at cosmological scales 
have been small, such that they can be treated as flrst-order perturbations [2]. Furthermore, 
the distribution of these perturbations (in the flrst order) is assumed to be statistically 
homogeneous and isotropic. Since the observable Universe is nearly homogeneous, and its 
spatial curvature either vanishes or is negligible until very near the present epoch, we will 
take the unperturbed metric to be flat FLRW (12.21) with K = 0. 

The most general perturbed FLRW metric can be written as 

ds^ = -(1 + 2A)dt'^ + 2a{diB + Vi)dx'dt + ((1 - 2C)6ij + 2dijE + 2d^iWj) + hij) dx'dx'J , 

(2.24) 

where di denotes a derivative with respect to x* where lower case Latin indices run over the 
three spatial coordinates. Due to the spatial isotropy and homogeneity of the unperturbed 
FLRW metric and energy-momentum tensor it is convenient to decompose the perturbations 

^ The above intuition, that modes become essentiahy classical at superhorizon scales, has been the basis 
for developing a powerful technique, the 5N formulation [16| . for dealing with the perturbations. We have 
reviewed this in Appendix [BJ 
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into scalars, divergence-free vectors, and divergence-free traceless symmetric tensors, which 
are not coupled to each other by the field equations or conservation equations, up to the 
linear order. In the perturbed metric, A,B,C and E are four scalars, while Vi and Wi are 
two divergence-free vectors and hij is a divergence-free traceless symmetric tensor. 

The most general first order perturbed energy-momentum tensor around the perfect fluid 
(12. 3p can be described as [2] 

=P5g,j + (^6,,{6P - ^VV^) + 9,,7r^ + d^nj + d.nY + nf^ , (2.25a) 

STio =PSgio - {diSq + Sq^) , (2.25b) 
SToo = - pSgoo + 6p . (2.25c) 

where p and P are the unperturbed energy density, pressure, respectively while 6p and 
6P are their corresponding perturbations. Furthermore, vr"^, vr^ and Sir^^j represent scalar, 
divergenceless vector and divergence- free, traceless tensor parts of dissipative corrections 
to the inertia tensor, respectively. In addition, 6q is the scalar part of the perturbed 3- 
momentum, while 6qY represents its divergence- free vector part. Note in particular that, 
the conditions vr'^ = tt^ = tt'^^- = describe a perfect fluid and Sq^^ = represents an 
irrotational flow. Being a perfect fluid or having irrotational flows are physical properties, 
thus their corresponding conditions are gauge-invariant. In other words, vr"^, vr^, vr-^^- and 6u\ 
are all invariant under infinitesimal space-time coordinate transformations. 

We have already used the rotational symmetry of the system for decomposing pertur- 
bations into scalar, vector, and tensor modes. The equations have also a symmetry under 
spatial translations which make it possible to work with the Fourier components of pertur- 
bations. Once we treat perturbations as infinitesimal, Fourier components of different wave 
number are decoupled from each other [2]. 

One important issue we need to care about is gauge degrees of freedom. In fact, when 
we define perturbed quantity, we have to make the following difference 

5Q{x\t) = Q{x\t) ~ Q{t) , (2.26) 

where Q{x^, t) is a real physical variable and Q{t) is a fiducial background variable. These two 
quantities reside in different spacetimes. Hence, we need to identify two points in different 
spacetimes to take the difference. Apparently, this is ambiguous, which is the source of 
the gauge degrees of freedom. Because of the presence of the gauge transformations which 
are generated by spacetime diffeomorphisms, not all metric and energy-momentum tensor 
perturbations are physically meaningful [2]. The gauge degrees of freedom may be removed 
by gauge-fixing (working in a specific gauge) or working with gauge-invariant combinations 
of the perturbations. In what follows we work out the gauge-invariant combinations of these 
modes. 

• Scalar modes 

Let us first focus on the scalar perturbations. This sector is the most involved and inter- 
esting one, involving eight scalars A, B, C, E, 6p, 6P, 6q and a^n^ . Of course, not all of 
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these quantities are gauge invariant. The gauge transformation in the scalar sector can be 
generated by the infinitesimal "scalar" coordinate transformations 

t ^i=t + 6t, 

. .. (2.27) 

where 5t determines the time slicing and 5x the spatial threading. In order to remove these 
gauge degrees of freedom, we can construct two gauge invariant combinations from the metric 
perturbations, the Bardeen potentials, 

m =C + a^HiE ) , (2.28a) 

a 

$ =yl - I {a\E - , (2.28b) 
and the four gauge invariant scalar parts of 5T^i, 

Spg =6p - pa^{E ) , (2.29a) 

(X 

B 

6Pg=6P-Pa^{E ), (2.29b) 

B 

Sqg =6q + {p + P)a\E ), (2.29c) 

ct 

and the anisotropic stres^ a^vr" [2], where 5q = {p + P)5u. 

Out of the ten perturbed Einstein equations, there are four scalars, two (divergence-free) 
vectors, and one massless tensor mode (gravitons). Among the four scalar perturbed Einstein 
equations, one is dynamical and three are constraints 

= $ - $, (2.30) 
5qg + 2{i! + H^) = Q, (2.31) 

5pg - 3H6qg + 2—* = , (2.32) 
6Pg + 6qg + 3H6qg + (p + P)$ = . (2.33) 

Note that the above equations do not form a complete set, unless the pressure P and 
anisotropic inertia a^vr** are given as independent equations [2]. For a general hydrody- 
namical fiuid with pressure P(p, S), the pressure perturbation 6P can be decomposed as 

6P = cl5p + 6S, (2.34) 

where 5S is the entropy perturbation and the sound speed is defined as = {^)s- For the 
adiabatic case (55" = 0, 5P is given directly by 5 p. 



^ Since the background energy- momentum tensor has the form of a perfect fluid, a^n^ is a gauge invariant 
quantity 
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It is common to construct two further gauge invariant combinations in terms of metric and 
energy- momentum perturbations. The comoving curvature perturbation "Jl^ and curvature 
perturbation on the uniform density hypersurfaces (: 

X = ^--^6qg, C = -^-^6pg. (2.35) 

p -\- 1 p 



Since 

2H k\ 
^ 3(p + P)a2 ' 

Jic and C become identical in the super- horizon scales ^ <^ 1. The crucial property of Jic and 
( is that in case of adiahaticity of the primordial perturbations, they are conserved outside 
the horizon. If cosmological fluctuations are described by such a solution during inflation, 
then as long as the perturbation is outside the horizon, "R^. and will remain equal and 
constant [2]. Later in this section we discuss more about the adiabatic and isocurvature 
fluctuations. 



• Vector modes 

We now consider the four divergence- free vector modes Vi, Wi, bq( and bi^^ . Again, not all 
of the above quantities are gauge invariant, but transform under infinitesimal vector gauge 
transformations induced by 

a;^ -> = X* + , (2.36) 

where di8x( = 0. The three gauge invariant divergence-free vector perturbations may be 
identified as 

= aW, - V, , (2.37) 

5q!{ and The perturbed Einstein equations have two independent vector equations, one 
constraint and one dynamical equation. These equations are 

di (^2aV}' - h^a^Zj)^ = , (2.38a) 

26qY + ^aZi = . (2.38b) 

Although not independent of the Einstein equations, it is useful to also present the vector 
part of the momentum conservation equation 

SqY + 3H5qY + Vvrf = 0, (2.39) 

which implies that, regardless of the matter content and the specific form of the anisotropic 
stress nY , dqj always damps like at large (superhorizon) scales. Thus, after horizon 
crossing, the flow always gets irrotational in inflationary systems. However, from (12.38^ we 

^Since the background fluid is a perfect and irrotational fluid, 5qY and ttY are gauge invariant quantities 
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see that in order to determine the dynamics of Zi we need to know irY ■ This term which can 
be non-zero in models of inflation involving gauge (vector) fields in general may change the 
usual picture. 

For the specific case of a perfect fluid in which nY = 0, equation f l2.38b ) implies that Zj 
as the only physical (observable) combination of metric components for vector perturbations, 
decays as ^. Consequently, 6qY decays like jj. In other words, vector perturbations are 
washed away by Hubble expansion, unless they are driven by a vector anisotropic stress 
irY ■ This source term is identically zero in scalar driven inflationary models. Due to their 
(exponential) decay, vector modes have not played a large role in these models. 



• Tensor modes 

Here we have two traceless divergence-free symmetric tensors hij and vr^^ and it is straight- 
forward to show that both are gauge invariant. The only field equation we have in this sector 
is the wave equation for hij (gravitational radiation), which is sourced by the contribution 

of TT^ij 

hij + 3Hhij + —hij = 2ttI . (2.40) 
a 

Being traceless and divergence-free, hij has 2 degrees of freedom which are usually decom- 
posed into plus and cross (-1- and x) polarization states with the polarization tensors e^'^. 
Since we assumed to have no parity-violating interaction terms in the action, the equations 
for both of these polarization have the same time evolution!^ and one may then introduce h 
variable instead 

= • (2-41) 

Then, in the computation of the power spectrum we consider this variable, treating it effec- 
tively as a scalar but multiply the power spectrum by a factor of 2 to account for the two 
polarizations. 

In the case that vr^ is zero, h becomes constant after horizon crossing (^ < if), then as 
long as the perturbation is outside the horizon, h will remain a constant [2]. 



2.2.1 Characterizing the primordial statistical fluctuations 

• Power spectrum of scalar models: The power spectrum of primordial curvature 
perturbation fl2.35p is one of the most practical and crucial statistical observables which 
may be used to distinguish models of inflation 

(ai,(k)3?,(k')) = (2vr)35(k + kOP^.(A;), ^ ^,P^Ak). (2.42) 

^ Cases with parity- violating contributions to tensor perturbations, which can lead to "cosmological bire- 
fringence" have been considered in the literature [171 HI] • We will return to this point in section |3l 



14 



where ( ■ ) denotes the ensemble average of the fluctuations. Furthermore, the scale- 
dependence of the power spectrum is described in terms of the scalar spectral index 

dlnA\ 

and the running of the spectral index by 

as ^ j^. (2.44) 

The special case with n<i = 1 and = corresponds to a scale-invariant spectrum. 

The power spectrum, Ag(A;), is often approximated by the following power-law form jT9] 

A2(fc) = A,(A;,)(-^)"^('=*)-^+^"4fc.)in(i)^ (2.45) 



where /c* is a pivot scale. 



• Adiabaticity of the power spectrum: Another quantity which can potentially offer 
an important test for models of inflation is the adiabaticity of the primordial perturbations. 
Roughly speaking, the adiabaticity may be defined as the following relation between the 
perturbations of Cold Dark Matter (CDM) or its various components and photons [2] 

^_3^ = 0, ) = 0). (2.46) 

Pni 4 Pr rir 

Then, any deviation from the above is defined as the isocurvature, or entropic, perturbation 

S„ ^ - Y-^. (2.47) 

Pm 4 Pr 

Since adiabatic and isocurvature perturbations lead to different peak structure in the CMB, 
they are distinguishable. CMB observations show that, if exists at all, isocurvature pertur- 
bations have to be subdominant [2U] . 

Whatever the constituents of the Universe, there is always an adiabatic solution^ of 
the field equations for which Jl^ and ( are conserved outside the horizon. If cosmological 
fluctuations are described by such a solution during inflation, then as long as the perturbation 
is outside the horizon, Jlc and ( will remain a constant [2]. 



• Non-Gaussianity: If the distribution of 3lc is Gaussian with random phases, then the 
power spectrum contains all the statistical information. The level of deviation from a Gaus- 
sian distribution is called non-Gaussianity and is encoded in higher order correlations func- 
tions of He- A basic diagnostic of non-Gaussian statistics is the bispectrum, the Fourier 
transform of the the three-point function of Jlc: 

{'Jl^{k)Jl,{k')Jl^{k")) = (27r)35(k + k' + k")Bji{k, k', k"). (2.48) 

^°Note that for these scalar modes, all individual constituents a of the Universe we have equal whether 
or not energy is separately conserved for these constituents. For this reason, such perturbations are called 
adiabatic. Then, any other solutions are called entropic. 
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which its momentum dependence is a powerful probe of the inflationary action. Note that 
here the 5-function is a consequence of translation invariance of the background and indicates 
that these three Fourier modes form a triangle. Different inflationary models predict different 
maximal signal for different triangle configurations and hence in principle the shape of non- 
Gaussianity can be a powerful probe of the physics of inflation |21j . 

One way to parameterize non-Gaussianity phenomenologically is through fj^i [22] 

D^,(x) = 3l,(x) + ^/^T' - {^»)] ■ (2.49) 

where 3ig(x) denotes the Gaussian curvature perturbation. Being local in real space, /^£'^' 
is a measure for the so-called local non-Gaussianity. 

In addition to /^£'^', the other commonly discussed parameter is the equilateral non- 
linear coupling parameter, /^^^*' which is defined through the power-spectrum normalized 
bispectrum for the equilateral configuration (/c ~ fc' ~ k") e.g. see [HI 1^ . 

Large non-Gaussianity can only arise if we have significant inflaton interactions during 
inflation. Therefore, in single-field slow-roll inflation, non-Gaussianity is predicted to be 
unobservably small [2U [25] , while it can be significant in models with multiple fields, higher- 
derivative interactions or nonstandard initial states [211 ISEl |27] . 

• Tensor modes: Power spectrum of the two tensor polarization modes h = h+ = hx , is 
defined as 

{h{k)h{k')) = (27r)35(k + k')P,ik), Al = —Pnik), (2.50) 
and the primordial gravitational waves spectrum quantify as 

A| = 2Al (2.51) 
The tensor scale-dependence is parameterized in terms of ut as 

tl 111 tv 

and the power spectrum itself is approximated as the following power-law form 

= A{k,){^r'^''\ (2.53) 

Finally, from the combination of fl2.45p and fl2.52p . we define another practical quantity, 
the tensor-to-scalar ratio r, as 

A2 

r ^ ^. (2.54) 

Prom the practical point of view, a tensor signal will be observable in CMB polarization in 
the upcoming observations if this ratio is bigger than 0.01 



^^The factor | in (|2.49p is conventional. Non-Gaussianity may be defined in terms of the Bardeen potential 
<&, $(x) = $g(x) + /^£'''($^(x) — ($^(x))), which during the matter era and at the linear order is <J> = 
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2.2.2 Computation of observable perturbations in inflation models 

Up to now we studied the cosmic perturbations without considering a specific model of 
infiation. At this point, considering the simple model of single-field infiation, we study the 
spectra of scalar and tensor perturbations generated in (12.61) . In terms of the gauge invariant 
Sasaki-Mukhanov variable v [221 EHj 

V = a{6y^ + ^C) , (2.55) 
perturbed scalar field equation of motion for a single scalar field takes the simple form 

v" + [k'' - —)v = 0, (2.56) 



where z = ^ and ' means derivative with respect to the conformal time r {dt = adr). Thus, 
the effective mass term ^ is approximately 

z" 

— ~ {aHf{2 + 3e - 2r/) ~ {aHf{2 + 5ey - 3r/y) . (2.57) 
z 

The general solution of this equation is then expressed as a linear combination of Hankel 
functions 



,(l+2.^W4(^^^^(l)+^^^(2)j^ (2.5^ 



where i/^^| + 2e — ?7 = | + 3ey — riy. Imposing the usual Minkowski vacuum state. 



2 1--./ 2 

-ifc 



in the asymptotic past ikr — )■ — oo), we obtain Ci = 1 and C2 = in (I2.58p . From the 
definition of the comoving curvature perturbation (I2.35P we find 

Jlc = — V. 2.60 

During the infiationary phase, fiuctuations in v generate scalar perturbations in the curvature 
Jlc. Then, the curvature perturbation Jlc becomes constant on the superhorizon scales and 
the power spectrum of the scalar fiuctuations at the end of infiation is 



A! ' 



2e (27r) 



(2.61) 

k=aH 



The spectral index of the scalar perturbations, to the leading orders in the slow-roll param- 
eters, is 

- 1 = -4e + 2r/ = -6ey + 2riv . (2.62) 

In the single- field slow-roll models non-Gaussianity is predicted to be small with f^L of 
order slow-roll parameters e and 77 [25]. Hence, we do not discuss them here. 
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Gravitational waves are also generated during inflation which introduce a tensor con- 
tribution to the fluctuations. For the simple scalar models the wave equation for tensor 
perturbations (12 .40^ reduces to a source-free graviton equation of motion, which is 



a" 



u" + ik' )m = 0, (2.63) 

a 



where u = ah and 

— = {aH){2-e). (2.64) 
a 

Using the slow-roll approximation and the Minkowski vacuum normalization, we have the 
same solution as (I2.58P with Ci = 1 and C2 = 0, with — \ + Finally, we get the power 
spectrum of the fluctuations in tensor components at the end of inflation 



(2.65) 

k=aH 



^-"(27r)2 

One can determine the tensor spectral index Ut which is 

nT = -2e. (2.66) 
Consequently, the ratio of tensor to scalar, r, is given as 

r - 16e = -Sut ■ (2.67) 
Recalling f l2.15p . we can relate it to the total inflaton fleld excursion during the inflation A(f, 

^~0(l)x(^). (2.68) 
Mpi ^ ' ^0.01^ ^ ' 

As mentioned before, a tensor signal is observable (in CMB polarization) only if r > 0.01. 
First showed by Lyth [31], equation f l2.68p implies that this level of tensor modes corresponds 
to super-Planckian Aip values. Using the effective theory of inflation, in [32], the Lyth-bound 
is generalized to all single-fleld models with two-derivative kinetic terms. They showed that 
the bound is always stronger than the above bound for slow-roll models. 



2.3 Observational constraints on CMB power spectrum 

We showed that the simplest (single-fleld) inflation models predict nearly Gaussian, scale- 
invariant and adiabatic scalar fluctuations. In this part, we briefly review the latest quan- 
titative constraints from the seven-year Wilkinson Microwave Anisotropy Probe (WMAP) 
observations [20]. Current observations provide values for power spectrum of curvature per- 
turbations A^, its spectral index n^, and impose an upper bound on the power spectrum of 
tensor modes A|,, or equivalently an upper bound on tensor-to-scalar ratio r. Moreover, it 
puts quantitative limits on physical motivated primordial non-Gaussianity parameters /]^£"' 
and /^^™'. These values vary (mildly) depending on the details of how the data analysis has 
been carried out. Here we use the best flt values of Komatsu et al. [20] which is based on 
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WMAP seven-years data combined with other cosmological data, within standard ACDM 
framework. Amphtude of curvature perturbations is obtained to be 



~ 2.5 X 10-^ 

and the data indicates a red tilted spectrum (n^ < 1) 

Us = 0.968 ±0.012 
as well small tensor-to-scalar ratio 

r < 0.24, 



(2.69) 



(2.70) 



(2.71) 



with no evidence of the running index, Furthermore, from the WMAP temperature 

fluctuation data we have the following constraints on the primordial non-Gaussianity param- 
eters 



- 9 < fj^l^' < 111, -151 < /X' < 253 . (2.72) 

The (non-)Gaussianity tests show that the primordial fluctuations are Gaussian to the 0.1% 
level [20], a strong evidence for the perturbative quantum origin of these perturbations. 

In [20], they explored the possibility of any deviations from the simplest picture: Gaus- 
sian, adiabatic, power-law power spectrum ACDM. However, they have not detected any 
convincing deviations from that. 
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Figure 1: In this flgure, we have the la and 2o observational contour bounds from the com- 
bined data of 7-year WMAP+BAO+HO. The points here represent the theoretical predic- 
tions of inflation models with monomial potentials, 'Vi}^) = Vo^p"' . The solid line represents 
the model with n = 4, the dashed line has n= 2 and the dotted line denotes the multi-axion 
fleld models with n= 2 and (3= 1/2 in [33], with e-folds Ng = 50 and 60 (from top to bottom). 
This flgure is taken from [2U] . 



Since the release of WMAP 3-year, the accessible parameter space for inflationary models 
has signiflcantly shrunk. In other words, the CMB data constrains inflationary models, and 
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even disfavored or totally ruled out some of them. For instance, in Fig{T] which is taken from 
[20] . the 7-year WMAP+BAO+HO observational contours on the spectral tilt, n^, and the 
tensor-to-scalar ratio, r, have been compared with the theoretical predictions of large field 
models with the monomial potential f l2.16p . The predicted point of the quadratic potential 
(n = 2) with Ne = 60 is on the boundary of the la contour, thus observationally preferable. 
However, predicting a large r, the quartic potential (n = 4) with N^. = 60 is totally out 
side of the 2a contour, unless iVg is sufficiently large. In fact, models with rather large 
tensor modes, such as X(j)^ chaotic inflation are disfavored now. Moreover, models with blue- 
tilted spectrum (n^ > 1), such as hybrid inflation models are ruled out now. Besides that, 
the possible parameter space for the models which can provide significant non-Gaussianity, 
such as K-infiation models with small sound speeds ^ 1 [M] are mainly shrunk by the 
constraints on the non-Gaussianity parameters. 

2.4 CMB anomaly and the eta problem 

As we mentioned, WMAP data strongly supports the infiationary scenario. It is believed that 
the primordial fiuctuations are statistically homogeneous, isotropic, and Gaussian. However, 
there seems to be some anomalies in the data prompting considering and analyzing various 
infiationary models. These "anomalies" include non-Gaussianity which we already discussed, 
a low power in the quadrapole moment [35l [36] , the alignment of the lower multipoles [37] , 
a five-degree scale cold spot with suppressed power [38], an asymmetry in power between 
the Northern and Southern ecliptic hemisphere [39], and broken rotational invariance [40] . 
These observational signature, although their statistical significance is still controversial [H] , 
stimulated activity in the research of statistical anisotropy [I2lllH|IMIl^ H ^ H7 t Hg | ll ^ [5U | 
\5T\. If the statistical anisotropy exists, a primary candidate of source of the anomaly should 
be gauge fields. Thus, from the observational point of view, it is well motivated to consider 
gauge fields in infiationary stage. Actually, as we will see in sections |3] and HI there are many 
phenomena induced by gauge fields. 

On the other hand, from the theoretical point of view, there is a reason to seriously 
explore a role of gauge fields during infiation. That is the eta problem, namely, scalar fields 
easily get radiative corrections which generically spoil slow roll conditions [52]. Hence, it is 
natural to explore a possibility to incorporate fields other than scalar fields to realize slow 
roll infiation. If we move beyond scalars, gauge fields are prime candidates. Indeed, as we 
will see below, there are at least two different mechanisms by which we can achieve slow roll 
inflation through the use of gauge flelds. 

3 Quantum gauge fields in infiationary background 

In this section, we explore possible effects of gauge flelds on inflation at the level of fluctu- 
ations. First, we review how to quantize gauge flelds coupled to the inflaton in inflationary 
background. We then discuss the old but important topic whether primordial magnetic 
flelds can be generated during inflation [531 El]- Recently, it turned out that there are 
other effects induced by gauge flelds. First of all, the presence of gauge flelds induces the 
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statistical anisotropy in curvature perturbations [55]. Indeed, the statistical anisotropy ap- 
pears both in the power spectrum and the bispectrum. Interestingly, gauge fields can induce 
detectable non-gaussianity in contrast to the single infiaton model with the conventional ki- 
netic term [25]. We review several mechanisms producing the non-gaussianity in primordial 
curvature perturbations from gauge fields. Moreover, gauge fields can generate primordial 
gravitational waves even in the low energy infiationary scenarios. 

3.1 Quantum gauge fields in cosmic background 

It is well known that gauge fields in 4-dimensions are conformally invariant at classical level. 
Due to the conformal invariance, gauge field fiuctuations, which are damped by the Hubble 
expansion as the vector perturbations discussed in the previous section, do not survive in 
the minimal set up. However, the non-minimal coupling described by the action 

S = -\j d^x^gf\(t))F,,F^^ , (3.1) 

may produce cosmological gauge field fiuctuations. Here, is the infiaton field and F^^y = 
dfj,Ay — dyA^ is a field strength of an Abelian gauge field A^. Let us consider gauge fields in 
isotropic and homogeneous background spacetime 

ds^ = a^{T) [-dr^ + 5ijdx'dx^] , (3.2) 

where r is the conformal time and a is the scale factor. We also assume the infiaton is 
homogeneous, cf) = </)(r) in the background. Therefore, /(</>) can be regarded as a time 
dependent coupling. Since the system is invariant under the gauge transformation — )■ 
Afj, + d^x with an arbitrary function x, we can work in the temporal gauge = and a 
gauge field can be expanded as 

= / (^^^("' '^^^"^ = S / (^^^^^)^^^k)e-- , (3.3) 

where polarization vectors (k) satisfy the transverse and normalization conditions 

A:,er(k) = 0, e^(-k)ef(k) = (5.., . (3.4) 

Here, a represents two independent polarizations. The wavenumber vector and polarization 
vector form a complete basis 

5^er(-k)eJ(k) = 5,,-Ml. (3.5) 
Substituting the expansion (13.31) into the action (13.11) . we obtain 

^ = \T.j - k'Al{r)A1,{r)\ , (3.6) 



21 



where a prime represents a derivative with respect to r. The coefficient A^, satisfies the 
equation 

Al" + 2^Al' + eAl = . (3.7) 
The canonical conjugate momentum is defined by 

< = ^ = /^^-. • (3.8) 

Now, we can quantize the system by promoting the fields to operators and imposing canonical 
commutation relations 

[yl-,vra=z5'^^-^5(k,-k2). (3.9) 
In terms of creation and annihilation operators, we can expand the operator as 

Al = Ukal + ula'-J, , (3.10) 
where [a^J, a^^^] = 6'^^'^^6{ki — k2). The mode function obeys the equation 

ul + 2^u'„ + k\k = . (3.11) 

In order to satisfy the canonical commutation relation fl3.9p . the mode function have to be 
normalized as 

d d i 

'^kg:^^! - ul—Uk = . (3.12) 

Once the mode function is determined, the vacuum can be defined by a^|0) = 0. Then, it is 
easy to calculate two point function 

{0\Aix)A\0)\0) = -j |^A:>,(r)|V'^- ^ j ^P^(A:, r)e^'^- , (3.13) 

where we have defined the power spectrum 

P^(A;,r) = -^K(r)^ (3.14) 

Similarly, we can define the power spectrum of the magnetic field Bi = eijkF^^ /2 by 

PB(fc,r) = -^|Mfc(r)|2 (3.15) 
and that of the electric field Ei = —A'- /a by 

PE{k,T) = -^Ju[{T)\' . (3.16) 
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The energy density of a gauge field can be expressed as 

f 2 f 

Pem{r) = Yj T + • ^^-^^^ 

Now, we need to determine mode function Uk- In the sub horizon hmit kr — ?■ — oo, we 
have a general solution 

As usual we choose the standard Bunch-Davis vacuum, i.e. to take the positive frequency 
mode functions as 

Uk ^e"^^^ . (3.19) 

fV2k ^ ' 



In the superhorizon limit kr — 0, we have a general solution 

dr 
1 



Ukr^di + d2 I ^ . (3.20) 



The constants di, d2 have to be determined by matching conditions. To perform this match- 
ing, we need to specify the time dependence of the coupling /(</)(r)). Here, we parametrize 
it as 



■2c 

a 



where c is a constant and a/ is the scale factor at the end of infiation. We assume that the 
background spacetime is de Sitter and the scale factor is given by a = 1/{—Ht). Then, with 
new constants, we obtain 

= + d^a"^"-^ . (3.22) 

Apparently, the two c > 1 /4 and c < 1/4 cases should be discussed separately. 

For c > 1/4, the second term is a growing mode. Hence, matching at the horizon crossing 
ttkH = k gives 

1 / 1 ^ 



d2 = ^ — . (3.23) 

Here, the suffix k represents the horizon crossing time of the mode k. Hence, fk should be 

\ -2c 



/. = ( ^ ) . (3.24) 



Thus, the power spectrum of the gauge field reads 
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Therefore, the power spectrum of the electric field is given by 

P.(A:,r,) = (4c-ir|;(^)"" , (3.26) 
and that of the magnetic field can be written as 

From these results, we see the gauge field and electric field survive during inflation for the 
parameter region c > 1. On the other hand, the magnetic field has no power on large scales 
unless c > 3/2. 

Next, we consider cases with c < 1/4. For these cases, the matching condition gives 

Thus, the power spectrum of the gauge field reads 

fajH 



The power spectrum of the electric field is given by 

PE{k,Tj) = 0, (3.30) 
and that of the magnetic field can be written as 

Thus, for the parameter region c < —1/2, gauge fields survive. In order for magnetic fields 
to survive, the parameter should be c < — 1. 



3.2 Primordial magnetic fields 

Based on the results in the previous subsection, let us first examine if primordial magnetic 
fields can be generated during inflation. Noting that the magnetic field has mass dimension 
two, namely IGauss = 10~^°GeV^, and assuming that Hubble during inflation is close to its 
current observational bounds, e.g. H ~ 10~^Mpi, just on dimensional grounds we can expect 

B^H^ ^ 10-i2^p2 _ io26GeV2 ~ lO^^Gauss . (3.32) 

After reheating, the energy density of magnetic flelds evolves as ~ 5^ oc 1/a^. As- 
suming the instantaneous reheating with maximal efficiency, i.e. all the energy of the 
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inflaton has turned into the thermal energy of a relativistic gas of particles [2], we have 
Ti = Ml,H\af/a,Y or 



Oo 



Tn 



1/2 



10 



29 



(3.33) 



where a/ and qq are the scale factor at the end of inflation and the present and Tq is the 
observed CMB temperature, Tq = 2.4 x 10~^^GeV. Thus, the expected magnetic fields at 
present should be 



5n = 5 X 10 



-58 



lO^^^Gauss. 



(3.34) 



To make the above estimation we have assumed the maximal efficiency for reheating. In 
reality, we should have ao/ct/ <C 10~^^. Hence, on dimensional grounds, we can obtain 
-Bo ~ 10~^ Gauss. Thus, it is a natural idea to generate primordial magnetic fields during 
infiation [M] • There are many works on primordial magnetic fields from infiation [521 ESI 
IFTl EHl ESI Eni ED E21 E3], although there exists no convincing model so far. 

In the previous subsection, we have discussed gauge fields in a fixed background. Here, 
we should note that previous analysis can be readily used in infiationary models. The most 
generic action for a single field inflation in this setup reads 



S 



-9 



M, 



pi 



2 i?- -9^09^0 -r(0)--/2(0)F^,F^^ 



In the slow roll regime, we need to solve 



mlH^ = v{<p) 



3F0 = -y^(0) 



Combining both equations, we get 



dloga 



pi y 



(3.35) 



(3.36) 



(3.37) 



which yields 

a oc exp 

Recalling (I3.2ip an appropriate choice for the coupling /(0) is 

f{4>) = exp 









V4> 



(3.38) 



2c 




[^Pi J 


v<t> 



(3.39) 



or more explicitly /(0) oc a For example, for a polynomial function V = 0", we obtain 



/(0) = exp 



(3.40) 
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Thus, with this choice for / we can directly use analysis of the previous subsection and 
inflation. 

Since / can be interpreted as an inverse of the coupling constant, large / means weak 
coupling. Hence, for positive c, the coupling starts from weak and becomes strong. We can 
normalize the coupling at the end of inflation. Then, the system stays in the weak coupling 
regime during inflation and the perturbative analysis is reliable. 

In the c > weak coupling case the electric fields dominates energy density of electromag- 
netic fields. On the other hand a scale invariant primordial magnetic fields requires c = 3/2. 
In this case, however, the energy density of electromagnetic fields grows as pem ~ {cLf/o.iY 
and we cannot neglect backreaction of the electric field produced through quantum fluctua- 
tions during inflation on the background inflationary trajectory. To avoid the backreaction 
we should set c = 1 and hence magnetic field cannot be scale invariant. In this case the 
primordial magnetic field has a blue spectrum 

B^y^^H'l^-^^ . (3.41) 

In the comoving scale 1 Mpc corresponds to k/{aH) ~ 10~^^ and hence 

B ~ lO^^Gauss (3.42) 

at the end of inflation and 

B ~ IQ-^^Gauss (3.43) 
at present. Apparently, this is too small. 

For negative c the effective coupling is a decreasing function of time. If we identify the 
coupling with the measured value at the end of inflation, the coupling during inflation is too 
strong and one may not rely on perturbative analysis. Despite the strong coupling issue, 
let us proceed with the order of magnitude estimate for the magnetic field. In the strong 
coupling cases c < the electric field can be negligible and to circumvent the backreaction 
problem we take c = — 1. Then, we obtain the primordial magnetic fields 

B ~ lO^^Gauss (3.44) 

at the end of inflation and 

B ~ IQ-^^Gauss (3.45) 

at present. Although the value of the primordial magnetic field in this case is much bigger 
than the c > case fl3.43p and close to the desired value, we should again stress that this 
result is not reliable due to the strong coupling problem. 

In conclusion, it is difficult to generate primordial magnetic fields if we want to avoid 
the strong coupling and backreaction problems. Later, we will discuss what happens in the 
analysis of primordial magnetic fields when we take into account the backreaction. Apart 
from a possible primordial magnetic field generation mechanism, one may explore the evolu- 
tion of the primordial magnetic fields after inflation and the effects they may have on other 
observables. Such analysis has been carried out in [Ml [651 166] . 
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3.3 Statistical anisotropy in power spectrum and bispectrum from gauge fields 

Besides the primordial magnetic fields, one may wonder if presence of gauge fields can in 
principle have observable effects on the CMB data. As discussed, for generic models of infla- 
tion vector fluctuations are generically suppressed by the cosmic expansion and to explore 
such a possibility we need to consider specific models of inflation with appropriate gauge 
field-inflaton couplings. Analysis of the previous subsection already provides a suggestion. 
To increase such effects one may consider models which allow a sharp change in the inflaton 
field, something similar to what happens in models of hybrid inflation, at the end of inflation 
(c/. discussions of section 12. ip . 

Let us consider a generic hybrid model with an inflaton 0, a waterfall field and a vector 
field yl^(/i = 0, 1, 2, 3) which couples to the waterfall field [55]. The action can be written as 

_1 1 d'x^gg^^g<^'^j\<l>)F,,F,, , (3.46) 

where F^j, = d^Ay — dyA^, ^(0, x, A^) is the potential of fields and an arbitrary function 
/(0) represents gauge coupling. Here, we do not restrict ourselves to the potential fl2.18p . 
The value of the inflaton at the end of inflation 0e is determined by parameters of this 
potential; in any case 0e and hence the coupling of the gauge field at the end of inflation 
/(0e) are fixed and do not fluctuate. As we discuss below, coupling other fields with the 
waterfall field x can, however, change the situation. To avoid standard problems in dealing 
with gauge field potentials, we assume that the vector field is massless and have a small 
expectation value compared to the inflaton. We hence neglect the terms coming from the 
coupling with the vector field in the background homogeneous equations of motion for the 
scalar fields and we treat the gauge field perturbatively. 

The curvature perturbation on superhorizon scales is given by the perturbation of e- 
folding number . In the standard single scalar inflation or hybrid inflation, in each 
Universe (causally connected Hubble patch), inflation ends when the inflaton reaches its 
critical value 0e which is determined by the inflationary potential. On the other hand, in 
the multi-component inflation, the critical value 0e may be different in each Hubble patch 
due to a light field other than inflaton 0. Hence, in such situation there is a possibility of 
generating the curvature perturbations because of the fluctuation of 0e. This mechanism is 
first proposed in [67]. We generalize their work and introduce a massless vector field A^ as 
another light field. 

From the above discussion it becomes apparent that to analyze the effects of fluctuations 
of 0e at superhorizon scales it is convenient to use the biN formulation [16]. Within the 
formalism the curvature perturbation generated at the end of inflation it = t^) can be 
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expressed as 

Cend 



ON 



dN d(f)e{A) 



5A' + - 



86, dA^dAi ^ d(hl dA' dA^ 



5A'5A^ , (3.47) 



where we set A^ = 0. That is, (pe can fluctuate due to the perturbation of the vector field 
A^. Let us take the hypersurface at the end of inflation t = to be a uniform energy density 
one. Then, the total curvature perturbation at the end of inflation t = is given by 



C(^e) Cinf ~l~ Cend 



where 



inf 



dN 



Id^N 



2 902 



(3.48) 



(3.49) 



is the conventional curvature perturbation generated by quantum fluctuations of the inflaton. 
Here, 0* represents a value at the initial hypersurface. 

The power spectrum of curvature perturbation At the leading order is given by 
(Ck.Ck.) = Pc(ki)5('Uki + k2) 

= iV.^P,(fc05(^)(k, + k,) + iV,2^^(M:(kOMi(k,)), (3.50) 

where we defined Pcj,{k) = H^/2k^, A'* = dN/dcj)^ and N^ = dN/dcpe- Here, we assumed 
that the scalar and the vector fields are statistically independent and Gaussian. The scale 
invariant power spectrum of the vector field {6A^{]<.i)6A^ {]<.2)) is obtained by setting c = 1 
in section 3.1 and given by 



(M:(k)Mi(k')) 



6'^ - 



(3.51) 



We note that to compute the above we have used canonically normalized field 5Ai = Ai/a 
where Ai is the coordinate basis field used in subsection 3.1. Using f l3.5ip . we obtain the 
expression for the power spectrum of curvature perturbations in the modified hybrid inflation 
model fl3.50p as 



Ae 



14 



P^{ki) , 



(3.52) 



where we have defined g,, = d(f)e/dA'', qij = d'^(f)e/dA^dAK The second term in the bracket 
above leads to statistical anisotropy in the power spectrum. It is usual to parameterize this 
anisotropy as [45], 



Pc{k) = pr{k) i+^7;3(q-k: 



(3.53) 
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where we used the notation q = q/ |q|, k = k/ |k| and the isotropic part 



/3) 



is separated. The coefficient of the anisotropic part reads 



(3.54) 



9P 



/5 



1 + /3 ' " \NJ, 

where we defined A''* = 1/ A/2e*. 

The bispectrum to the leading order is given by 

(Ck.Ck.Cka) = (27r)-3/2^^(ki,k2,k3)<5(=')(ki + k2 + k3) 

= (27r)^3/=^|iV2^„ [P^{ki)P^{k2) + 2 perms] ^^^^ (ki + ka + kg) 



(3.55) 



dA^ dAi \N,dA^^dA^^ dA^^ dA^^ 

X [{SAi{k,)SAiik2) ((5A^^^M^^)Jk3)) + 2perms]| ,(3.56) 

where -k in the last line denotes the convolution. From these expressions, we see the cur- 
vature perturbation has the direction-dependence due to the vector field. Substituting the 
expression (I3.5ip into (13. 56 p . we can deduce the bispectrum fl3.56p as 



(Ck.Ck.Cka) = (27r)-3/2^c(ki,k2,k3)5(^)(ki + k2 + k3 

5^(ki,k2,k3) = 



N,V ( \ 

mo 



Nee 



5«1 _ 



^2 



xP^{ki)P^{k2) + 2 perms 



(3.57) 



where we assumed that 5Aj is Gaussian. 



Recalling the anisotropy in the power spectrum (13.521) it is convenient to define the non- 
linear parameter /tvl as the bispectrum normalized by the isotropic part of power spectrum 
Pi'°{k) 



^/7VL(ki,k2,k3 
5 



5^(ki,k2,k3 



Pl'°{ki)Pi'°{k2) + 2 perms 



i^(ki,k2)fc| 



2 perms 



(3.58) 



F(k„k2) = ii + ^r^ + gf'' 



where we assumed the scale-invariant power spectrum and defined 

"l-(q-ki)^ 
g^i-(q-ki)^^ 
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Ne \qf 



l-(q-k2)^ 
- (q ■ k 



2;'^2 



(3.59) 



Strictly speaking, one should define the non-Gaussianity parameter /tvl normalizing the 
bispectrum with the full power spectrum instead of the isotropic one. However, since we 
expect the anisotropic part in the power spectrum to be small, the above definition provides 
a good approximation. We will discuss this point further in the end of this subsection. 

Under the slow-roll approximation, we obtain N^^:/N^ ~ N^e/N^ = 0{e). Hence, neglect- 
ing the first and second terms in the right hand side of the above equation, we have more 
simplified expression as 



i^(ki,k2 



g^^-(q-ki)A;^ 



«2 



(qk 



(3.60) 



where 



and the anisotropic part as 



hi2/ H ■ We can decompose the non-linear parameter into the isotropic part 



^/^z.(ki,k2,k3) = y'^i{i + r\^iMM)) 



(3.61) 



where the isotropic part reads 



5 



riso — /o^ „^ 



and the anisotropic part is deduced as 



r™(ki,k2,k3 



-q' U{(q ■ ki) + ^^2(q " ^2) + k\ki{(i ■ ki)(q ■ k2) 



^3 



(3.62) 



+2 perms 



(3.63) 

Here, we used the relation Ne = l/y/2€^. Taking a look at the above formula, we notice 
that the statistical anisotropy gives a non-trivial shape to the bispectrum even for the local 
model. 



Let us now consider a specific hybrid model f l2.18p with a charged waterfall field x iii the 
unitary gauge. The potential for this model is given as 



1 

2' 



1 

— / 
2 



(3.64) 



where Ai = 6ijA^ is the gauge field. The coupling constants are denoted by \,g,h, the 
infiaton mass is given by m, and the vacuum expectation value for x is represented by v. 
For this model the effective mass squared of the waterfall field x is given by 



Inflation ends when 

As discussed the critical value 0e depends on A*. We hence have 

1 / ^2^2 



q 



lAI 



q'' = {^^S'' + q'q' 



(3.65) 



(3.66) 



(3.67) 



}■ 
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and 



/3 



1 fh^\A\ 



(3.68) 



where we used the approximation ~ A^^e- We can now write down the power spectrum 



Pc(k) = pr{k) 



(q ■ k)^ 



(3.69) 



Notice that despite being direction-dependent, the power spectrum is scale invariant. The 
magnitude of statistical anisotropy is determined by the parameter (7/3. 

Let us now study the bispectrum which is more interesting. Substituting the expressions 
f l3.67p and f l3.68p into fl3.62p . we have the isotropic part of the /nl parameter 



6 . 

rzso 
'^J NL 



2 / 1 I 9 



2a2 



h^\A\ 



(3.70) 



where rj = V^^/V . Similarly, the anisotropic part fl3.63p reads 



r™(ki,k2,k3 



+ 1 



-,2 J,2 \ -1 



h^\A\ 




2 / . \ 2 

q ■ ki ) + q ■ ks 



q ■ ki ) q ■ ks ) + 



^,3 



,2 J,2 



2 perms 



9 



/l2|A| 



ki ■ ks ) q ■ ki I q ■ k 



1,3 



+2 perms 



(3.71) 



where as usual ki + k2 + ka = 0. From the above expression we see that the amplitude of 
the non-Gaussianity depends also on the magnitude of statistical anisotropy gp. For a 
large /3 (/3 ^ 1), (7^ ~ 1 and hence the statistical anisotropy appearing in the primordial 
power spectrum f l3.54p is large. On the other hand, for a small /3 (/3 <^ 1), the statistical 
anisotropy is small and the non-linear parameter /ati, is also small. 



Eg. (13.62]) provides the possibility of generating large non-Gaussianity even for the cases 
with the small statistical anisotropy in the power spectrum if we choose to be small 



and to be large while gp <^ 1. From f l3.6ip and (I3.7ip we see that for this case the 

anisotropic and anisotropic parts of the bispectrum are of the same order, in contrast to 
the power spectrum. Hence, it may be possible to detect the statistical anisotropy in the 
bispectrum with the future experiments and it will give us information about a new physics in 
the early universe associated with the violation of the rotational invariance. The mechanism 
explained in this subsection has been expanded and extended in many different ways which 
may be found inlMlESlEniiniESlESlEllESlIZS]. 
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3.4 Non-Gaussiantiy induced from gauge fields 



It is well known that a single field inflation with standard kinetic term cannot generate the 
sizable non-Gaussianity [25]. Presence of other flelds, including gauge flelds, may however 
change this result and generate detectable non-Gaussianity. In this subsection we consider 
two such setups in which gauge fields can enhance the bispectrum of the inflaton fluctuations. 
One such possibility is when gauge fields have an axionic coupling. The other one is when 
the gauge kinetic term has a non-trivial infiaton dependence. As we will discuss these two 
cases lead to two different shape non-Gaussianities, the former to an equilateral shape and 
the latter to a local shape. 



3.4.1 Non-Gausssianity from gauge fields with an axion coupling 



In the presence of axion coupling [T7], there exists tachyonic instability in the positive helicity 
polarization mode of gauge fields. The amplitude of gauge fields can hence be large, inducing 
large non-Gaussianity in curvature perturbations [771 EH [79]. To see explicitly how this 
happens, let us consider the Abelian U(l) gauge theory coupled to an axion 0: 



S 



d X^ 



-9 



Ml 1 

— R - -da 

2 2 ^ 



4 Oja 



(3.72) 



where fa is an axion decay constant and a is a constant. In the infiationary background, the 
gauge field in the Coulomb gauge obeys the following equations 



Ja 

To diagonalize equations we use circular polarization vectors e\ satisfying 

k = , ^ijkkj^^ = ^ike_^ , 
ei(-k) = 6i(k)% 6Uk)V,,(k) = 5,v. 

In terms of the polarization vectors, we can expand the gauge field as 

d^k 



(3.73) 



(3.74) 



A=± 



(27r)3/2 



[ax{k)AxiT,k)e\{k)e'''-^ + h.c] 



(3.75) 



where annihilation and creation operators aA(k),a^(k) satisfy the standard canonical com- 
mutation relations. The equations of motion fl3.73p boils down to 



k'T 



2ki 



A±(r,fc) = o, e 



acj) 



(3.76) 



Here, we neglected the deviation from de Sitter expansion and also assumed ^ is almost 
constant during slow-roll infiation. 
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In the interval (8^) ^ < k/aH < 2^, there exists tachyonic instabihty in the positive 
hehcity polarization modes where we can approximate the mode function as 



A+{T,k) 



k 



2k \2^aH 



1/4 



(3.77) 



Since negative helicity polarization modes do not suffer from the tachyonic instability, here- 
after we will neglect them. The exponential growth caused by the tachyonic instability is 
physically the same phenomenon as the usual mode-freezing at the horizon crossing and thus 
we have production of gauge fields due to axion coupling. 

The produced gauge fields in their own turn can source scalar fluctuations. To see this, 
let us decompose the infiaton fluctuations as 



d^k Qkjr) 
(27r)3/2 a{T) 



(3.78) 



For our analysis below it is enough to consider zeroth order in slow-roll parameter where 
a(r) = and hence the modes Qkij) satisfy 



+ k^ 



dr 



Qk{r) = Jkir) 



where the source term Jk, which is induced from the axion coupling, is 



Jzir) = a^(r) 



a 



d^k 



fa J (27r)3/2 



(3.79) 



(3.80) 



Because of the tachyonic instability, this source term can induce significant effect on the 
statistics of curvature perturbations. The solution consist of homogeneous Q™'^ and inho- 
mogeneous Ql ones 

Q,{r)=Qr{r) + Qi{T) . (3.81) 
The homogeneous solution can be expressed by the positive frequency mode function 



TX 



In order to obtain inhomogeneous solution, we use the retarded Green's function 
Gfc(r, r') = id{T - r') [Mrmr') - vlir)Mr')] 

satisfying 



,2 2' 



dr 



Using the Green function, we obtain 

QI{t)= f dT'G,{T,T')MT') . 



(3.82) 



(3.83) 



(3.84) 



(3.85) 
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Recalling the relation 



Ji^ = -E.9l , (3.86) 

a 



we can calculate the power spectrum 

(3i,(k03^e(k2)) = 4^ [{QT:QT:) + {QIQD] ■ (s-sr) 

Here, we used the fact that and Ql have no statistical correlation. Apparently, the first 
term is the conventional one 

(3ir(k03ir(k2)) = ^^^'(ki + k^) . (3.88) 

The contribution from the inhomogeneous terms is given by 

{Jiiik,)Jilik,)) = S^ [ dT'dT" ^'^J^l'^ {JkAn)JuAr,)) . (3.89) 

02 J a[T) a[T) 

For large ^, we can deduce 

P3^, = P [1 + P/2(e)e^'^T , (3.90) 
where /2(0 can be numerically calculated. For 2 < ^ < 3, the best fit becomes 

3 X 10-5 
J' 

Similarly, one may compute the correction to the bispectrum 



/2(0 = ■ (3.91) 



f Gkjr, n) Gkjr, t^) Gfc(r, rg) 

^ a(r) a(r) a(r) 
^{JkAri)JkAr2)Jk,{n)) ■ (3-92) 

Defining new variables 

|ki| = k , |k2| = X2k , Iksl = x^^k , (3.93) 



we obtain 



{^l{k,m{k2)^l{k,)) = -A(2vr)5/W<l±4^/3(e,X2,X3) ^^^i±i^^,(3.94) 

where f^{^,X2,x^) is a complicated function. Since the bispectrum should have equilateral 
configuration, we only need 

7.4 X 10"^ 

/3(e,l,l)= ^s., • (3.95) 
From the bispectrum, we can read off 

Taking into account the current observational bound —214 < /^l"'' < 266 obtained from the 
WMAP data, we get the constraint ^ < 2.65 [78]. 
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3.4.2 Non-Gaussianity from the gauge kinetic function 



Another possibihty for generating sizable non-Gaussianity from gauge fields discussed in the 
literature is through the gauge kinetic term [801 Elj- Let us consider the following model 



S 



M. 



<P-V{<P)--f{<i>)F,,F^^ 



(3.97) 



where instead of tachyonic instability, the gauge kinetic function generates gauge field quanta. 
Let us again consider the slow roll limit of infiation. Thus, the calculation is mathematically 
similar to the previous case and the difference is the source term. That is, the perturbations 
Qk{T) are governed by the equation fl3.79p but now the source term is 



J,(r) 



2 / J (27r)3/2 



(3.98) 



In the case of axion coupling, the exponential growth caused by tachyonic instability enhances 
the amplitude of the gauge field, which works at the horizon crossing. This is why the shape 
of the bispectrum is equilateral. Here, instead, the non-minimal coupling makes the gauge 
fields survive, which works on superhorizon scales and persists during infiation. Thus, the 
shape of the bispectrum should be local. 

The procedure for the calculation is exactly the same as the axion coupling cases. Defining 
the power spectrum 



27rJ 



Pa,,(A;i)5(ki-k2 



we obtain the result 



Pj,^ = P[1 + 192PiVcMB (iVtot - iVcMB)] 



(3.99) 



(3.100) 



where P is the conventional power spectrum, A^^cmb is the number of e-folds from the CMB 
scale to the end of infiation, and A^tot is the total number of e-folds of infiation. In order for 
the perturbative analysis to be valid, the correction must be subdominant 



192PiVcMB (iVtot - iVcMB) < 1 . 

Since P ~ 10^^, this can be easily satisfied. 
Next, we can calculate the bispectrum 

iXi^Mk^Wk-s)) = Bj,^{k,)5{k, + k2 + kg) 



(3.101^ 



(3.102) 



and the local non-Gaussianity can be characterized by /nl defined in fl2.49p . The result is 
intriguing 



fNhiki 



/•equiv. local 

/nl 



kf cos^(A;2, ks) + A;| cos'^^k^, ki) + kl cos'^{ki, /C2) 

i-S I 1,3 I t.3 
rv-^ ru2 \ 



(3.103) 
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where 

^equiv.iocai ^ isgOPiV^^B {N,,, - Ncmb) ^ 0.7 ^ (iVtot - iVcMe) (3.104) 

In spite of the constraint f l3.10ip . we can obtain sizable non-Gaussianity in the range 

10 < /nl < 100 . (3.105) 
Moreover, for ^3 ^ A;i ~ A;2, the fc-dependence of the bispectrum becomes 



3) 



oc coseK° + sineK° , e = tan"^ ~ 0.22 .(3.106) 
" 2^5 

That is, there is anisotropy in the bispectrum. 

3.5 Primordial gravitational waves from gauge fields 

Having a non-zero contribution from the gauge fields to the energy momentum tensor, they 
will appear as the source term for (primordial) gravity waves. This could be more pronounced 
when the gauge fields are produced due to the tachyonic instability induced by a time 
dependent axion, as discussed in section 3.4.1. So, let us consider the action (13. 72 p . It is 
straightforward to show that the gravity waves are governed by 



1 
2^ 



+ 2 
dr"^ a dr 



1 



hi, = -T72[EiEj + B,B,] , (3.107) 



where the source term should be calculated taking into account the tachyonic instability. As 
we see the two graviton polarizations appear with different source terms. This will cause 
a birefringent (chiral) gravity wave production which has its own interesting observational 
prospects [TTJ [791 Ell |83] . Although there are several mechanisms to generate birefringent 
gravitational waves [SU the tachyonic instability provides a natural and simple setup. In 
particular, it gives rise to primordial gravitational waves even if the energy scale of inflation 
is low. Hence, it is interesting to see the signature in the CMB [86] . 

Here we only review the effects of the source on the total gravity wave power spectrum, 
which is calculated to be 



2H^ fk 



GW 



7r2Mp2 \k^ 



TIT 



tt2 



(3.108) 



MO ^ ' ,5.7 • (3-109) 



where the approximate fitting function for the relevant range 2 < ,^ < 3 reads 

2.6 X 10"^ 
J' 

Here, the first factor is the standard amplitude for source-free gravity waves (I2.65P while the 
extra contribution is due to the gauge fields. This extra contribution could be large even in 
the low energy inflation. Besides the parity violation in gravity, one may also analyze non- 
Gaussianity induced by such gravity wave. Such an analysis has been carried out in [87[ [88]. 
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4 Anisotropic inflation with gauge flelds 

To be consistent with the current observations the degree of anisotropy at cosmological 
scales and on the CMB data should be quite small, of the order of 10~^ or less. Hence, the 
accelerated expansion during inflation should be practically isotropic. This point has been 
put as the cosmic no-hair conjecture [891 190]. The classic attempt in providing a proof for the 
conjecture was made by R. Wald [7] which goes under the name of cosmic no-hair theorem. 
This theorem, as discussed in [8] and reviewed in Appendix [Dl need not be obeyed by the 
inflationary models. So, it is interesting to explore possibility of having inflationary models 
with anisotropic hair. 

In general, since turning on vector gauge fields in the background gives rise to a preferred 
direction, one would naturally explore building anisotropic inflation models involving gauge 
fields. In the previous section, we discussed that gauge fields may be used to describe a 
possible statistical anisotropy in the CMB data. In our discussions, however, the gauge fields 
were turned on at the fluctuations level and their existence led to the anisotropic couplings 
between gauge fields and curvature perturbations. It is then natural to explore presence of 
gauge fields at classical background level and study how they can affect inflationary trajectory 
and other observables within anisotropic inflationary models. 

4.1 Historical remarks 

The first attempt to use vector fields as inflatons was proposed by Ford [91]. He used 
a potential of vector fields to realize exponential expansion of the universe. After that, 
there appeared other attempts to construct anisotropic inflation. However, because of the 
prejudice due to the cosmic no-hair conjecture, and the assumptions of Wald's theorem in 
view [^, people tried to make contrived models which generically break energy conditions. 
For example, higher curvature terms are used to violate the cosmic no-hair conjecture [92] . 
Recently, a non-minimal coupling between vector fields and curvature is adopted to make the 
vector fields slow roll [51in3]- Unfortunately, these models suffer from (ghost) instability, or a 
fine tuning problem [B]. After all, it turns out that violation of energy conditions, more than 
what is already required in inflationary models, is not necessary to realize anisotropically 
accelerating cosmologies [8] (see Appendix |D]) . The gauge kinetic term in supergravity is a 
necessary ingredient to realize anisotropically accelerating expansion. Indeed, anisotropic 
inflation could be ubiquitous in supergravity [9^ [95] . 

As shown in the previous section, gauge fields can play important roles in inflation. In 
these calculations, typically, the backreaction is not negligible. Indeed, the energy density 
of gauge fields is kept almost constant during "anisotropic" inflation models and anisotropic 
inflation provides a self-consistent framework for calculating various observables. A a matter 
of terminology, the inflationary setup in which the gauge kinetic term is /(0)^-F^ is named 
anisotropic inflation or "hairy inflation" 
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4.2 Inflation with anisotropic hair 



In this subsection, we study the bosonic part of the action in supergravity and explain how 
anisotropic inflation appears in the models with the action 



S 



(4.1^ 



where we set = 1 and V{(f>) and /(0) represent a potential of the inflaton and the gauge 
kinetic function, respectively. Note that the gauge kinetic function can contribute to the 
potential energy as is found in a different context [96l [97] . The picture we envisage is as 
follows. As we have seen in the previous section, gauge flelds can be produced at superhorizon 
scales and their spectrum could be red. Therefore, a coherent, space independent gauge fleld 
can be consistently turned on in the background and affect background evolution. We now 
explore the consequence of the coherent vector fleld. 

Without loosing the generality, one can take x-axis in the direction of the vector fleld. 
Using the gauge invariance, we can express the vector fleld as 



Af.dx" 



v{t)dx 



(4.2) 



Thus, there exists the rotational symmetry in the y-z plane. Given this conflguration, it is 
convenient to parameterize the metric as follows: 



ds^ = -J^{tfdt^ + e 



2a(t) 



(4.3) 



where e°, a and cr„ are an isotropic scale factor and spatial shears, respectively. Here, the 
lapse function 3\f is introduced to obtain the Hamiltonian constraint. With the above ansatz, 
the action becomes 



S 



dx-e 



3a 



1 ■ 



1 



2„-,2„-2a(t)+4a(t) 



(4.4) 



where an overdot denotes derivative with respect to the comoving time t. First, its variation 
with respect to (j_ yields 

cr_ = — 3d6'_ . (4.5) 

This gives (t_ oc e"^"", hence, the anisotropy in the y-z plane rapidly decays as the Universe 
expands. Hereafter, for simplicity, we assume (t_ = and set the metric to be 



ds'' = -de + e2°(*) [e-^-^Wrfx^ + e^'^^'\dy^ + dz' 

which is a Bianchi type-I model (see Appendix O for more discussions). 
The equation of motion for gauge fleld v is easily solved as 



(4.6) 



-a— 4a 



Pa 



(4.7) 
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where pa is a constant of integration. Taking the variation of the action with respect to 
[NT, Q!,(T and and substituting the solution (14.71) into them, we obtain the following basic 
equations 



a 



a 



1 



Pa 



-3a' + Vi<P) + '-^f{<P)-'e-'--'^ 
o 

Pa 



3 

= _3«0_K^+p^/(0)-3/^e 



-4«— 4fT 



(4.8) 

(4.9) 

(4.10) 
(4.11) 



where as before the subscript in denotes derivative with respect to its argument (f). 

Let us now check whether inflation occurs in this model. Using (14.81) and (14. 9p . the 
equation for acceleration of the cosmic expansion is given by 



3^ 3 



y _ Pa ^-2 -4«-4(t 

2 



(4.12) 



We see that to have (slow-roll) inflation the potential energy of the inflaton should dominate 
over the energy density of the vector field 



(4.13) 



and the shear S = 6". 



First, we need to look at the shear to the expansion rate ratio S/if, which characterizes 
the anisotropy of the inflationary universe. Notice that (I4.10p reads 



-3i7S + -p, . 



(4.14) 



If the anisotropy converges to a value, i.e. E becomes negligible, the terminal value should 
be given by 

^ = 3 ^0) ' ^^-''^ 
where we used (14. 8 p in the slow-roll approximation, i.e. 

1, 



a 



(4.16) 



In order to realize the above situation, p„ must be almost constant. Assuming the 
standard slow-roll approximation and that the vector field is subdominant in the evolution 
equation of the inflaton field (14. lip , one can show the coupling function /(0) should be 
proportional to e~^" to keep py almost constant. In the slow roll phase, the e-fold number a 
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is related to the inflaton field as (ia = — l^(0)(i0/V0, as usual. Then, the functional form 
of is determined as 

/(0) = e-2" = e'^i'^ (4.17) 

For the polynomial potential V oc 0", for example, we have / = e~ . The above case is, in 
a sense, a critical one. What we want to consider is super-critical cases. For simplicity, we 
parameterize /(</>) by 

/(0) = e"^i'^ (4.18) 

where c is a constant parameter. (The above may be compared with the analysis in subsection 
13.21 and in particular to f l3.40p .) 

Let us consider c > 1 case. Note that f l4.18p can be written as 

^ = 2c^. (4.19) 
Then, the condition c > 1 can be promoted to the condition 

^ > 1 . (4.20) 



2 fV 

Thus, any functional pairs / and V which satisfies f l4.20p in some range could produce the 
vector-hair during inflation. The equation for the inflaton becomes 



-Sad) — Va, 



2c 

ev V{ct>) 



(4.21) 



where we used (14.190 and the slow- roll parameter ey is defined as fl2.14p . In this case, 
if the vector field is initially small pv/V{(l)) -C ey/2c, then the conventional single field 
slow-roll inflation is realized. During this stage / oc e~^'^" and the vector field grows as 
oc e^^'^"^)". Therefore, the vector field eventually becomes relevant to the inflaton dynamics 
fl4.2ip . Nevertheless, the accelerating expansion of the Universe continues because pv/V{(t)) 
does not exceed ey jlc. In fact, if p„/'K(0) exceeds ey l'2.c, the inflaton field (f) does not roll 
down, which makes = p^/(0)~^e~'^"~^'^/2 decrease. Hence, p„ -C ^(0) always holds. In 
this way, there appears an attractor where the inflation continues even when the vector field 
affects the inflaton dynamics. We stress that even if new infrared perturbations become 
relevant, those may be renormalized to the background vector field and the anisotropy stays 
at the attractor value. 

Let us make the above statement more precise. Note that the inflaton dynamics is 
determined by solving the slow-roll equation: 

- 3d0 -V^ + vlr^Ue-^''-'" = . (4.22) 

Here, we see gauge kinetic function changes the effective potential. Using the slow-roll equa- 
tion (14.160 . this yields 

^ = ^ = + 2c^e-'"-^^-^^^^'^^ . (4.23) 
da a V Va, 
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This can be integrated by neglecting evolutions of V, and a 

' = f-T^ + l]e-^(-i)"+^'^] , (4.24) 



where f2 is a constant of integration. Substituting this back into the slow-roll equation fl4.23p . 
we obtain 

^ = _Y± + h + ^e-4(c-l)«4-4.1-l _ (4_25) 

da V c V 

Initially a — >■ — oo, the second term can be neglected. While, in the future a — >■ oo, the term 
containing Q disappears. This clearly shows a transition from the conventional single field 
slow-roll inflationary phase, where 

to what we refer to as the second inflationary phase, where the vector field is relevant to the 
inflaton dynamics and the inflaton gets 1/c times slower as 

^ = AYt U27) 
da cV ■ ^ 

In the second inflationary phase, we can use fl4.24p discarding Q term and rewrite the energy 
density of the vector field as 

P. = ye Jy,'' = -^evV{<P) , (4.28) 



which yields the anisotropy 

S Ic- 1 



ev . (4.29) 



H 3 

Moreover, from (14. 8 p and (14. 9 p the scale-factor slow-roll parameter (12. lip is 

a 1 / 1 2 \ 1 , , 

where we neglected the anisotropy and used (]4.16p and (]4.27p . Thus we have a remarkable 
result inn 

S Ic- 1 , , 

Jf = -—- (4.31) 

Therefore, for a broad class of potential and gauge kinetic functions, there exist anisotropic 
inflationary solutions, and the anisotropy is of order the slow-roll parameter e. 
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4.2.1 Example: chaotic inflation 

In order to make the statement more precise, we consider chaotic inflation with the potential 

V{<P) = ^m202 , (4.32) 

where m is mass of the inflaton. For this potential, the coupling function (I4.18P becomes 

= e'^^''^ . (4.33) 
It is instructive to see what happens by solving fl4.8p - fl4.lip numerically [91]. In Figj2l 

2e-006 


■2e-00e 
-4e-006 

dt 

-6e-006 
-8e-006 



-le-005 

-2 2 4 8 8 10 12 

Figure 2: Trajectory of the inflaton fleld (j) in phase space. This trajectory is drawn for 
parameters c = 2 and m = 10~^ and with initial conditions (pi = 12 and <pi = 0. There 
are two different slow-roll phases, namely isotropic and anisotropic inflation. The transition 
occurs around = 9. 

we have shown a trajectory in — space where, as our earlier analytical discussion also 
showed, we see two slow-roll phases. The flrst one is the conventional isotropic inflationary 
phase and the second one is the anisotropic inflationary phase. As usual, inflation ends 
with oscillation around the bottom of the potential. This tells us that isotropic inflation 
corresponds to a saddle point and ansiotropic inflation is an attractor in the slow roll phase. 
In FiglHl we have calculated the evolution of the degree of anisotropy S/if = cr/d for various 
parameters c under the initial conditions ^/c(j)i = 17. As expected, all of solutions show a 
rapid growth of anisotropy in the first slow-roll phase which corresponds to the conventional 
inflation. However, the growth of the anisotropy eventually stops at the order of the slow roll 
parameter. Notice that this attractor like behavior is not so sensitive to the parameter c. As 
one can see there is a sufficient amount of number of e-folds during an anisotropic inflation. 
We should emphasize that the anisotropy of the expansion disappears after the inflation. 
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Figure 3: The time evolutions of the anisotropy S/if for various c with respect to the 
number of e-folds. The anisotropic inflation phase can be sufficiently long. 



Hence, the anisotropy affects only the generation process of primordial fluctuations and the 
evolution of fluctuations during inflation. 

As we have proved for general cases, the anisotropy satisfies the inequality 



This inequality holds universally for any potential functions. This result exemplifies the 
extension of Wald's cosmic no-hair theorem in infiationary setup discussed in [8] and reviewed 
in Appendix [Dl It is noteworthy that deviations from a de Sitter space during a slow-roll 
infiation, the isotropic and anisotropic deviations alike, are labeled and characterized by the 
slow-roll parameter e. 

4.2.2 Anisotropic infiation in a variety of models 

There are many models which realize anisotropic infiation. We can generalize single field 
infiation to multi-field infiation models |98] or readily show that the small field or hybrid 
infiation models also admit anisotropic extension along the line we discussed here. A wider 
range of anisotropic infiationary models are discussed in [99]. We can extend the standard 
kinetic term to the Born-Infeld type |100j which may be useful in finding stringy realization 
of anisotropic infiation. We can introduce a mass term to the vector field, that is, a vector 
curvaton [101] . It is interesting to study cosmological consequences of the vector curvaton 
scenario in detail. It is also possible to extend the model to non-Abelian gauge fields [OSl 
llU2j . In this case, we have more complicated dynamics which would lead to interesting 
phenomenology. In a straightforward way the Bianchi type-I model we discussed here can 



E 



(4.34) 
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be extended to other Bianchi type models |103j . (For a discussion on Bianchi cosmologies 
seeO) 

Instead of a gauge field, which is a one-form, one may realize anisotropic inflation through 
higher form fields |104j . For example, we can consider the action for a 2- form field -B^j/ 



S = d'^Xy/^ 



(4.35) 



where H^„\ = dfj,B^x + d^Bx^ + dxE^j^y is the field strength of the 2-form. Moreover, the 
parity violating term 



(4.36) 



can also induce anisotropy in the expansion. Here, the axionic coupling function j(0) should 
be chosen appropriately. 



4.3 Exact power-law anisotropic inflation 

In the previous subsection, we have seen anisotropic inflation appears in many different 
forms and scenarios. To gain more profound understanding of the mechanism of anisotropy, 
it would be useful to have exact solutions which exhibit anisotropic expansion |105] . We 
know there exists exact isotropic power law inflation for the exponential potential 



y(0) = e^"^ . 

Hence, it is natural to consider the exponential kinetic function 

/(0) = 

in order to obtain exact anisotropic power law inflation. 



(4.37) 



(4.38) 



4.3.1 Exact solutions 

We start by recalling that isotropic power-law solutions can be obtained through the ansatz 

a = Clogt, = eiogt + </.o , (4.39) 

where 



<4 



corresponding to FLRW geometry 

ds^ = -dt^ + t^/>^' {dx^ + dy^ + dz') 



(4.40) 



(4.41) 
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Thus, for A ^ 1, we have power-law slow- roll inflation. 

We now seek exact anisotropic solutions. We start with the power-law ansatz 

a = Clogt, a = 7]\ogt, </. = eiogt + 0o . (4.42) 

From the hamiltonian constraint (14. Sp . we get two relations 

Ae = -2 , pe + 2C + 2r/ = 1 , (4.43) 

to have the same time dependence for each term. The latter relation is necessary only in the 
non-trivial vector case, pa 0. Then, for the amplitudes to be balanced, we need 

-e + v' + ]:e + l:u + ]:w = 0, (4.44) 
b 6 b 

where we have deflned variables 

u = Voe^''"' , w = plfo^e-^P*° . (4.45) 
The equations for the scale factor (14.91) and the anisotropy (I4.10p . under (I4.43P yields 

-C + SC-u-]:W = 0, (4.46) 

-ri + 3Cr]--w = . (4.47) 
3 

Finally, from the equation for the scalar (14. lip we obtain 

-^ + 3C^ + \u~ pw = . (4.48) 
Using (I4.43p . (I4.46P and (I4.47p . we can solve u and w as 

9 9 ^ 3p ^ 1 P / . .^N 

= -9C= + ^c + fc-^-f . (4.50) 

Substituting these results into (14.481) . we obtain 

(3C - 1) [6A (A + 2p) C - (A^ + 8pA + 12p2 + 8)] = . (4.51) 

li C = 1/3, we have u = w = which is not our desired anisotropic solution. Thus, we 
have to choose 

_ A^ + 8pA + 12p^ + 8 

^" 6A(A + 2p) • ^^-^^^ 

Substituting this result into (I4.46P and noting (14. 43 p . we obtain 

A^ + 2pA-4 
^ = 3A(A + 2p) ' 

e = (4.54) 
45 



Finally, (09|) and f g30|) reduce to 

_ (pA + + 2)(-A^ + 4pA + 12p^ + 8) _ (A^ + 2pA - 4)(-A^ + 4pA + 12p^ + 8)^ ^ 

2A2(A + 2p)2 ' 2A2(A + 2p)2 ^^'^^^ 

Note that fl4.44p is automatically satisfied. Thus, we have obtained exact anisotropic power- 
law solutions. 

Recalling the definition (14. 42 p . we see C ^ 1 or |A| ^ 1 is necessary for slow- roll infiation. 
From the solution f l4.52p . it turns out that this requirement can be achieved if A <^ p. For 
these cases, u is always positive and since w is by definition also positive, 

A^ + 2pA > 4 . (4.56) 

Hence, p must be much larger than one. 
The spacetime metric is 

The average expansion rate is determined by C, and the average slow roll parameter (12. lip is 
given by 

^ _^ ^ 6A(A + 2p) 

d2 A2 + 8pA + 12p2 + 8 ' ^ ^ ^ 

In the slow-roll limit A ^ 1 and p ^ 1, this reduces to e = A/p. 
The anisotropy is characterized by 

2(A^ + 2pA-4) 

H a A2 + 8pA + 12p2 + 8 ' ^ ' ' 

where if = d is the isotropic expansion rate. One may rewrite the above in the form (I4.3ip 

E c - 1 A^ + 2pA 



if 3c ' 4 



> 1 . (4.60) 



We see the anisotropy is positive and proportional to the slow roll parameter e, with a 
proportionality constant which is always less than 1/3. 

To summarize, we have constructed a model of anisotropic power-law infiation which 
can be exactly solved. Although the anisotropy is always small, it persists during infiation. 
Clearly these exact solutions give rise to counter examples to the cosmic no-hair conjec- 
ture [891190]. 



4.3.2 Stability of the anisotropic infiation 

We showed that a model with (I4.37P and (I4.38P admits both isotropic (I4.40p and anisotropic 
power-law solutions. Here, we investigate the phase space structure and stability of these 
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solutions. It is convenient to choose e-fold number as the time coordinate da = adt and use 
the dimensionless variables 



a 



Y = - 

a 



Z = /(0)e 



-a+2a 



a 



With these definitions, we can write the hamiltonian constraint equation as 

Since we are considering a positive potential, we have the inequality 

3{X^ - 1) + ^Y^ + ^Z^ < . 



(4.61^ 



(4.62) 



(4.63) 



Using the hamiltonian constraint f l4.62p we can eliminate (p from the equations of motion 
and write them in the autonomous form 



dX 
da 
dY 
da 
dZ 
da 



-Z\X +1)+X |3(X2 - 1) + -Y^ 
3 12 

(F + A) UiX^ - 1) + + -YZ^ 



3{X' 



1) + -Y^ - pY + 1-2X + -Z^ 
2 3 



(4.64) 
(4.65) 
(4.66) 



Therefore, we have a 3-dimensional space with a constraint (14.631) . A fixed point in this 
phase space is defined by dX/da = dY/da = dZ/da = 0. 

We first analyze the isotropic fixed point X = 0. From f l4.64p we learn Z = 0. The 
remaining equation fl4.65p yields Y = —A or Y'^ = 6. The latter solution does not satisfy the 
constraint ( I4.63p . Thus, the isotropic fixed point is 



(x,r,z) = (o,-A,o) . 



(4.67) 



This fixed point corresponds to the isotropic power-law solution (I4.4ip . Indeed, one can 
check that the solution fl4.4Up leads to the above fixed point @ 



Next, let us discuss an anisotropic fixed point. From (I4.64p and fl4.65p . we have 



Y 



Z' 



3p+^]X-X 



3X 



3{X^ - 1) + -Y^ 



X+1 

Using the above results in (I4.66P we obtain 

(X - 2) [(A^ + 8pA + 12p2 + 8) X - 2 (A^ + 2pA - 4)] = 



(4.68) 
(4.69) 

(4.70) 



^Note that the fixed curve Z = and + — 6 does not satisfy the constraint (|4.63p . 
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The solution X = 2 does not make sense because it implies Z'^ 
and fl4.69p . Thus, an anisotropic fixed point is expressed by 



-18 - 36p2 < by fOSjl 



X 
Y 
Z' 



2 (A^ + 2pX - 4) 
A2 + 8pA + 12p2 + 8 ' 
12 (A + 2p) 
A2 + 8pA + 12p2 + 8 ' 
18 (A2 + 2pA - 4) (-A2 + 4pA + 12p2 + 8) 
(A2 + 8pA + 12p2 + 8)^ 



(4.71) 
(4.72) 
(4.73) 



Note that from the last equation, we find that A^ + 2pA > 4 is required for this fixed point 
to exist under the condition A ^ p. This nicely dovetails with the discussions of previous 
subsection 14. 3 . 1 1 and this fixed point corresponds to the anisotropic power-law solution (I4.57p . 

Now, we examine the linear stability of the fixed points. Eqs. (14.641) . (14.651) . (14.661) in the 
linearized form are given by 



d6X 

da 
d6Y 
da 



d6Z 
da 



3 



-r^ - 3 + XY5Y + - (X + 1) Z5Z , 



(4.74) 



6X {Y + X)6X+ |3 (X2 - 1) + ly2 ^ r (y + A) + ^Z^j 6Y 



+ l^-Y + 2p + \j Z6Z , 

2(3X - 1)Z6X + {Y -p) Z6Y 
1 



3X' 



-Y' 



2X - pY - 2] 6Z . 



(4.75) 



(4.76) 



First, we consider the stability of isotropic fixed point (I4.67p . In this case the above equations 
reduce to 



dSX 
da 

d5Y 
da 

d5Z 
da 



h'-3]6X, 



6Y 



pX 



6Z. 



(4.77) 
(4.78) 
(4.79) 



We see that the coefficient in the right hand side of above equations becomes negative when 
A^ + 2pA < 4 during infiation A ^ 1, which means the isotropic fixed point is an attractor 
under these conditions and the isotropic fixed point becomes stable in this parameter region. 
In the opposite case, A^ + 2pA > 4, the isotropic fixed point becomes a saddle point and 
unstable. In the latter case, existence of a background gauge field, no matter how small, 
destabilizes isotropic infiation. Next, let us discuss the stability around the anisotropic fixed 
point. Since we are considering the infiationary universe A ^ 1, the condition A^ + 2pA > 4 
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Figure 4: Several trajectories in X-Y-Z space are shown for A = 0.1, p = 50. The trajectories 
converge to the attractor corresponding to anisotropic inflation. 



implies p ^ 1. Under these conditions the linearized equations are approximated as 

-3SX, 

-36Y + v/6(A2 + 2pA - 4)6Z , 



da 
d6Y 

da 
d6Z 



l-^/6{\^ + 2p\-4)6Y . 



da 2 
The stability can be analyzed by setting 

5X = e'^'^SX , 5Y = e'^'^SY 

Then we find the eigenvalues cu are given by 



5Z = e'^'^SZ 



id 



|± J3(A2 + 2pA-4)-^. 



(4.80) 
(4.81) 
(4.82) 

(4.83) 
(4.84) 



As the eigenvalues have negative real part, the anisotropic fixed point is stable. Thus, the 
end point of trajectories around the unstable isotropic power-law inflation are anisotropic 
power-law inflation. 

In FigJU we depicted the phase flow in X-Y-Z space for A = 0.1, p = 50. We see 
that the trajectories converge to the anisotropic fixed point indicated by yellow circle. The 
isotropic fixed point indicated by orange circle is a saddle point which is an attractor only on 
Z = plane. Thus, anisotropic power-law inflation is an attractor solution for parameters 
satisfying A^ + 2pA > 4 [TU5] . 
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4.4 Inflation with multi-vector hair 



In the previous subsection, we have considered a single Abelian vector gauge field. It turns 
out that there are qualitative differences in the dynamics when the number of Abelian vector 
fields becomes more than two |1U6] . For example, isotropic inflation with vector fleld hair 
can be realized dynamically (as an attractor solution) in the multi-vector fleld system with 
a uniform coupling between an inflaton and gauge flelds |106[ 1107] . 

We consider the exponential potential with parameters Vq and A and N copies of Abelian 
gauge fleld A^T"* with fleld strength Fji^^ = df^A^J^^ — d,yA^f^^ described by the action 



m=l 



(4.85) 



The N copies of Abelian gauge fleld A^j^^ are coupled to an inflaton </> with coupling constants 
gm- As before, the Greek letters denote the full space-time indices and the Latin letters are 
reserved for spatial parts. 

For the analysis in this subsection, it is convenient to use the tetrad formalism |108] . As 
used in Appendix |Dl let us introduce the unit normal n'^ of the homogeneous hypersurface in 
a generic Bianchi I spacetime. The evolution of the spatial slice is governed by the equation 



n 



H{g^,u + n^Uy) + a^^ , (4.86) 

where H is the averaged expansion rate and the trace-free tensor a^i, represents anisotropy 
of the universe. Similarly to (E25]), the energy-momentum tensor can be decomposed as 

T^u = pn^riy + 2g(^ni.) + Pig^^u + n^n^) + , (4.87) 

where p, p, and n^y are respectively energy density, pressure, energy flux and anisotropic 
pressure seen by the observer with 4-velocity n^'. The electric flelds seen by an observer with 
the four-velocity are deflned by 

F(-) = Fj™)n'^ (4.88) 

and we assume, for simplicity, that the magnetic parts vanish 

e^.a/3F^'")"^n'^ = , (4.89) 

where e^ya^ is the volume four-form (see Appendix |A] for our conventions). The Einstein 
equations then take the form 

1 ^ 

+ Voe^'^ + - XI e^'-^F^"^)^^'^)^ (4.90) 





m=l 



P 

q, = (4.91) 



P 



I, 



2 Voe^"^ + - V ef'-'^F^'")^^™)'^ (4.92) 



2^ " 6 ^ 

m=l 



Nr. 



m=l 



(4.93) 
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Here, an overdot denotes a derivative with respect to the proper time t associated with an 
observer with proper velocity n^. 

In an orthonormal frame, the Einstein equations take the form given in |109j . We assume 
that ^ and E^^ are non-vanishing and not parallel to each other and choose the spatial 
frame such that 

Ef^ = (E,0,0), Ef^ =Q. (4.94) 
This fixes the choice of directions in spatial sections and results in 

Vti = (723, ^^2 = -0-13, ^3 = 0-12, (4.95) 

where Qi represents rotational velocity of the frame with respect to the Fermi propagated one 
(for a discussion on Fermi transport see e.g. |110j ). Let us adopt the following dimensionless 
variables 

1 , , . „ CT,; 

= ] 

and 



S± = ^(^22±a33), ^«i = ^ for i^j, (4.96) 



where capital Latin superscript indices from 3 to N. Then, the evolution of geometry is 
governed by the following equations 

SV =(g - 2)S+ - (S?2 + S?3) + 2£2 + 2(£f ))2 - (£f )2 

TV 

+ ^(2(£S^))^-(£(^))^-(£(^))^), 

yl=3 

TV 



S'_ =(g - 2)S_ - S?, + + 2S^3 - 3(£f )^ - 3 $^((£f ^ - (£^ 

A=3 

N 

=(g - 2 + 3S+ + S_)Si2 + 2E13E23 - 6£f ^£?) - 6 ^ £[^^£^ 



(4.98) 



^3 

A=3 

N 

^2 ) 
A=3 

N 

=(g - 2 + 3S+ - S_)Si3 - 6 el^^Ss^^ , 

A=3 

N 

S'23 =(g - 2 - 2S_)S23 - 2S12S13 - 6 £2^^£3^^ ' 

A=3 

where a prime represents a derivative with respect to the "e-fold number" time coordinate 
r defined by 

dr = Hdt , (4.99) 

(note that, being in an anisotropic background, r is not the conformal time) and g = — 1 — 
H/H^ is the deceleration parameter 

3 N 

g = 2S^ + |(S^„ + SL + S?3 + Sy + V-n + £^ + (£f))^ + (£r)^ + 5^5^(£l^))^ (4.100) 

i=l A=3 
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(4.101) 



One may also work out equations of motion for the scalar and gauge fields 

3 TV 

n' = {q- 2)n - 3An + sg^e,' + 39,{{Ef^)' + (£f )n + 39A J2 

i=l A=3 

n' = (2g + 2 + \n)u , 

£'=(g-l-f^-2S+)£, 

(£f))' ={q-l-^Q- 2S+)£f ^ + 2Si2£f , (4.102) 

(£f)'=(g-l-|l] + S+ + S_)£f , 
(£S^))' =(q-l-3An- 2S+)£S^) + 2Si2£^^^ + 2Si3£i^\ 

(£(^))' =(g _ 1 _ + S+ + S_)£^^) + 2S23£i^\ (4.103) 

(ef^y -ffi + s^-s_)£r^ 

Moreover, there exists a constraint among the variables 



3 Af 



+ £2 + (g(2))2 ^ (£(2))2 ^ J2 $^(£rV = 1 . (4.104) 



1=1 A=3 



We are now in a position to study the fate of anisotropy in inflation with multiple vector 
fields hair. 



4.4.1 Universe tends to be isotropic 

Here, we consider the cases of uniform coupling constants g = g„i and show that isotropic 
inflation is an attractor in the phase space when the number of vector fields is greater than 
two. 

The equations (14.981) to (I4.103P form a dynamical system of dimension 3(A^ + !)• To 
understand the dynamics, it is useful to find out equilibrium points and their linear stability. 
Let us look for equilibrium points with non- vanishing £, £3^^ and £3'^''. For their time 
derivatives to vanish, we need 

g-l = |fi, S+ = S_ = 0. (4.105) 

From the evolution equations for £[^'' and £["^\ we also obtain 

S12 = Si3 = S23 = . (4.106) 
Defining the total energy density parameter of the electric fields 

£2 = £2 + (£S2))2 + (£(2))2 ^ ^(£(^))2 (4.107) 

A 
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and using the equilibrium conditions for Q and 11, we can deduce relations 

\-g „ 4 



n 



A(A + (?)-4 



(4.108) 



Apparently, we need the condition X{X + g) > A. In order to have an accelerating universe, 
we also require g < or, equivalently — g^ < 0. What remains is to determine the relative 
strengths and angles among the vectors through S'^^ = 0. There are six equations for the 
remaining 3(A^ — 1) variables 



3 ) ; 



For = 3, they lead to an orthogonal solution 

A(A + ^?)-4 



.(2)n2 



(2)n2 



£^ = (£^ 



g(2) ^ g(3) 



£ 



(3) 



. 



A(A + ^)-4 

(4.109) 



(4.110) 



3iX + g)^ ' 

For more general cases > 3, it is convenient to introduce three N — 2 dimensional vectors 



Then, the magnitude of £3 is given by 

|£3| 



o(4) 



A(A + (7)-4 



(4.111) 



3(A + ^?)2 ' 

and, the magnitudes of £2 and £1 are given through relations 

(er)^ + |£2 



(4.112) 



(2)^2 , I? ,2 _ A(A + ^) -4 



2 ' 



£^ + (£fV + |£i 



3(A + (7) 
A(A + (7)-4 



(4.113) 



3(A + (7)2 • 

They introduce three arbitrary constant parameters. We also know from 

£1 • £3 = £2 • £3 = 



(4.114) 
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that £3 is perpendicular to £1 and £2. The angle between £1 and £2 is fixed once we choose 
the three parameters controling their magnitudes by using 

£f^£f + £1 ■ £2 = 0. (4.115) 

Apart from these geometrical conditions, we can take arbitrary combinations of the com- 
ponenets for these three N — 2 dimensional vectors. Due to the 0{N — 2) invariance of the 
equilibrium values of E^c^\ they are a family of equilibrium points that spans a 3(A^ — 3) 
dimensional submanifold in the state space, for which the expansion rate becomes isotropic 
even though the vector fields have non-zero background values. 

To examine local stability of an equilibrium point, we need to examine linearized equa- 
tions around it. It is not so difficult to see that the problem can be reduced to the stability 
analysis of the orthogonal equilibrium point with three electric fields. Hence, we shall do the 
stability analysis for = 3. Since 11 can be eliminated by using the Hamiltonian constraint, 
there remain the twelve independent varibales which are grouped into six pairs: T,i2,8,\ , 
Si3,£f^; S23,£?^; S_,£^^^ - £f ; S+, 2£ - fif - £f and, Q, 8, + E^^ + 8,f. The sta- 
bility analysis gives the following results. The first five pairs, which represent anisotropic 
perturbations, share common eigenvalues given by 

g - 2 ± v/(g - 2)2 - 16£=^ 
c^s,£ = • (4.116) 

Their real parts are negative since — 1 < g < 2. Furthermore, the isotropic mode has 



un,. = ^-2±V(g-2)2-4£2(4 + 3,(A + ,)) ^ 

If we require q = < 0, it follows 4 + 3g{X + g) > 0, which means that both of the 
eigenvalues have again a negative real part. In other words, these isotropic equilibrium points 
are local sinks for a range of the parameters for which the universe undergoes accelerated 
expansion. 

One may confirm the above result beyond linearized stability level by solving the dynami- 
cal equations f l4.98p - f l4.103p numerically, for the case of three vector fields. We solved the ba- 
sic equations with the parameters X = 2, g = 5, and the initial conditions 'E± = S12 = S13 = 
S23 = 0,n = 0.6, £ = 0.1, £S^^ = 0.2 , £5^^ = 0.001 , £f^ = 0.1 , £f ^ = 0.01 , £f ^ = 0.0001. 
The initial condition for 11 is determined by the constraint equation. The parameters are 
chosen so that the universe is accelerating, namely g < 0. From FigJSl we see the anisotropy 
disappears after transient anisotropic infiationary phases. This confirms that the isotropic 
infiation is an attractor. It should be noted that the vector fields possess a nontrivial or- 
thogonal configuration. We have solved the basic equations for other sets of parameters and 
initial conditions. Although the transient behavior depends on the parameters and the ini- 
tial conditions, the system always asymptotically approaches the isotropic infiationary final 
state. The transient anisotropic phases correspond to the saddle points. However, since du- 
ration of anisotropic infiation is sufficiently long, the anisotropy at each of the saddle points 
would be relevant to CMB observations in the realistic cases where infiation ends with a 
finite duration. It is worth mentioning that non-Abelian gauge fields also exhibit a similar 
isotropization |lllj . 
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anisotropy 




Figure 5: For a uniform coupling g = 5, we plotted time evolutions of anisotropic expansion 
rate normalized by the Hubble parameter with respect to the e-fold number time r. The 
components S+, S_, S23, S13, and S12 correspond to thick, dashed, blue, red, and green 
lines, respectively. 



4.4.2 The cosmic minimum-hair conjecture 

Based on the results and experiences of anisotropic models presented so far, we put forward 
the "cosmic minimum-hair" conjecture stating that the hair in inflation tend to be minimized 
by the dynamics of the system in a quasi-de Sitter expanding Universe. We provide sup- 
portive evidence for this conjecture through studying more general anisotropic trajectories. 
A closely related model independent analysis is also presented in Appendix 

It is natural to ask whether isotropic inflation with non- vanishing vector fields is a unique 
attractor of the system. It is difficult to extract a clear conclusion for a system of this 
complexity. It is also a parameter dependent problem. However, we can argue that multiple 
vector fields are expected to repel each other and try to become isotropic in an expanding 
Universe. Let us first look at the stability of axisymmetric infiating solutions discovered in 
|105] . In the present model, they are located on lower dimensional boundaries of the full 
state space. For example, we have an equilibrium point 

_ 5A^ + 2Xg -3g^-8 2(A^ + - 4) 12{X + g) 

^ iX + 3g){X + g) + 8' + (X + 3g){X + g) + 8' (X + 3g){X + g) + 8' 

E_ = S12 = = S23 = eS'^ = = 

3((A-3^)(A + ^)-8)(^(A + ^) + 4) 3(A(A + ^) - 4)((A - 3^)(A + g) - 8) 

[(X + 3g){X + g) + 8]^ ' [^x + 3g){X + g) + 8]^ 

There are many others which represent physically the same spacetime, but lie on different 
boundaries. This was a stable attractor solution for the single-vector-field model. For this 
equilibrium state, we notice that E_|_ > by looking at the evolution equation for E_|_ and 
requiring £^ > 0. More generally, any electric field in 1-direction tends to support positive 
by generating a tension along that direction. It is also easy to see that for any electric 
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perturbation in 2- or 3-direction, the eigenvalue is given by 3S+ > 0. This axisymmetric 
solution is hence a saddle point in the general multi-field state space. These points together 
with the mathematical structure of the Maxwell's equations imply that the positive Z1+ 
created by the 1-component of a field comes with negative sign in the evolution equation in 
1-direction while it has plus signs in the other directions. Thus, if there are more than one 
vector fields and one of them is dominant, it creates such an anisotropy that destabilizes 
orthogonal components of the other fields. This instability does not show up for single field 
models since it would merely cause a rotation of that vector. In summary, we expect multiple 
vector fields to rearrange their orientations to minimize anisotropy of the space which would 
cause an instability in some of their components. 

In the present model the vectors cannot decouple from cosmic dynamics because they are 
kept excited through coupling to the dominant scalar inflaton field. However, they still re- 
distribute themselves to achieve as much isotropy as possible within the given circumstances. 
As an example, if = 2, the attractor solution contains two orthogonal electric fields with 
the same amplitude: 

_ 2A^ + 2Xg -3g'-2 _ A(A + (?) - 4 ^ _ 3A(A + (?) - 12 



{X + 3g){X + g) + 2' ^ 2{X + 3g){X + g) + 4' " 2{X + 3g){X + g) + 4' 
^ ^ 6{X + 2g) ^ 3igiX + g) + 2)((?(2A + 3g) + 2) 



(A + 3g){X + g)-2' ((A + 3g){X + g) + 2^ 



.(2).2_ 3(A(A + r?)-4)(^7(2A + 3(7) + 2) 



2((A + 3^7)(A + ^7) + 2)2 ' 

Si2 = Sl3 = S23 = £1^=0. 

Shall we included a third gauge field, it becomes unstable because of the eigenvalue for £3 

cj£3 = 6S+>0. (4.118) 

From the above argument, we can understand the global structure of the phase space for 
cases of uniform coupling constants. There are a bunch of saddle points consisting of those 
axisymmetric signle-field and orthogonal two-field equilibrium solutions on the boundaries, 
which are attractors when restricted in properly chosen invariant subsystems. Since these 
invariant sets have dimensions smaller than 3(A^ + 1)? the initial condition from which the 
orbits are attracted to any of the anisotropic attractors is of measure zero in the entire state 
space. Whenever an orbit starts from a point not included in these subsystems, it is attracted 
to the isotropic final state after a sufficiently long time. For example, in order for an orbit 
to be attracted to the two field attractor, it must satisfy £ -"^^ = for all A = 3, ■ ■ ■ , A^. Any 
small deviation from this subsystem would lead it to isotropy due to the linear instability 
demonstrated above. A typical orbit is first attracted towards a nearby saddle point. Then 
the instability explained above kicks in and another vector component rises. It might come 
across another instability and go to another saddle point. An orbit continues this routine 
until it finally settles down to one of the isotropic attractors. 

In order to see what happens for general couplings, we took couplings gi = 4.8, g2 = 
5.0, = 5.2 for the three gauge fields. As one can see in FigjGl the attractor is an anisotropic 
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anisotropy 




Figure 6: Time evolutions of anisotropic expansion rate normalized by the Hubble parameter 
versus the e-fold number time r. The plots are for the couplings gi = 4.8, g2 = 5.0, = 5.2. 
The components S+, S_, S23, S13, and S12 correspond to thick, dashed, blue, red, and green 
lines, respectively. 
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Figure 7: Time evolution of the total anisotropy S^'^S^,y/2 with respect to the e-fold number 
time r. 
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inflation. In this general case, the behavior is even more interesting. Since there are several 
saddle points, we have different anisotropic inflationary stages. During inflation the pre- 
ferred direction in the Universe is changing, which would lead to interesting observational 
consequences. For example, the anisotropy can be scale dependent. Even if we change initial 
conditions, the attractor is the same. Hence, the predictive power of the theory is main- 
tained. In FigJTl we plotted the square of the magnitude of anisotropy Yj^^Yj^^yjl. While 
the system is not close to any of the fixed points, the orbit follows a complicated dynamical 
trajectory. However, once it starts to be attracted towards one of the inflating solutions, the 
evolution is fairly regular. The result shows monotonic decrease of the total anisotropy after 
the orbit goes into this regular tracking regime despite the fact that each component of the 
anisotropy may exhibit a complicated behavior. It is inferred that when an orbit goes from 
a fixed point to another, the latter is always less anisotropic than the former. 

The above evidence supports the cosmic minimum-hair conjecture: the Universe organizes 
itself so that any anisotropic feature in the spacetime during inflation becomes minimum. 

4.5 Primordial magnetic fields revisited 

Noting that space-independent background gauge field configurations we discussed in this 
section and the backreaction of superhorizon gauge field configurations analyzed in section 
|3] have the same mathematical appearance in the equations of motion, we can return to 
the question of primordial magnetic fields of section |3] and address the backreaction issue 
discussed there using the elaborate computational tools of previous subsections. 

Whenever the produced primordial magnetic field has a scale invariant spectrum c = 3/2, 
f l3.27p implies that the electric fields have a red spectrum. Hence, the dominant contribution 
to the energy density comes from large scales. We can assume the coherent electric fields with 
a definite direction dominate the energy density of the electromagnetic fields. Contrary to 
the simple expectation that the energy of the electromagnetic field will eventually overcome 
that of the inflaton and terminate inflation [61], as we discussed, growth of the electric 
field initiates anisotropic inflation [Ml 1112] and one should revise and improve this simple 
picture by considering its backreaction effects. Explicitly, before the backreaction becomes 
important, we have the relation 



However, backreaction effects changes the dynamics of the inflaton such that (cf. fl4.24p ) we 
find the attractor behavior 



That is, the effective Ceg changes from Ceff = c to the critical value Ccfj = 1 due to backreaction. 
In other words, the scale invariant spectrum of a gauge field is an attractor [99l 1112] . Using 
fl3.16p . (I3.17P and (13.220 . before the backreaction becomes important we have 




(4.119) 




(4.120) 
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where H is the Hubble during inflation which is almost constant, and hence the transition 

I' 



point ah occurs at Pem ~ Pa ~ 10 "^p^ ~ 10 '^M^^H'^, i.e. 



-10-^^. (4.122) 

We can now calculate the power spectrum of magnetic flelds. First, we consider the 
modes which exit the horizon before ab- The superhorizon evolution of the mode function 
before is given by 

1 fa 



u^r) = --= - , (4.123) 

fkV2k \akj 

where we should note fk is deflned by ft = {ab/ afY'^~'^{ak/ af)~'^'^ from the continuity at a^. 
Since the evolution after a^ becomes Wk oc a^, we obtain the mode function after as 

.. = ^f*V''7^V. (4.124) 



From f l3.15p and that = H'^Pb, we obtain magnetic fields 

Tj2 / \ 2c-2 



B{X„rf) = ^{X,Hr~U^ , (4.125) 

where Tf denotes end of inflation and Xp = af/k is the physical wavelength of the mode at 
the end of inflation. 

Compared to the cases with no backreaction fl3.27p . the amplitude is reduced by the 
factor 

\ 2c-2 / N 2c-2 

-) -10-^^f^l . (4.126) 
afj H \afj 

For the flat spectrum c = 3/2, we can deduce magnetic flelds at the end of inflation as 

B(.,,.,).10-.§i(^)^. (4.127) 

Without backreaction, we anticipated 10~^^Gauss for the scale invariant case c = 3/2. How- 
ever, by taking into account the backreaction, we have a suppression factor ab/af in f l4.125p 
which is about 10~^^. Hence, we can expect at most 10~^^Gauss on Mpc scales at present. 

For modes which exit the horizon after the transition time ab, by setting c = 1 in fl4.123p . 
we obtain 

1 fa 
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Figure 8: The magnitude of magnetic fields is schematically depicted as a function of 
wavenumber k for the case c = 3/2. There exists a break at kb = abH in the spectrum 
due to the backreaction. As can be seen, the amplitude of magnetic fields on Mpc scales 
gets a suppression by 10~^^ due to this break. 



where fk is now defined by fk = {ak/aj) ^. Thus, we can calculate magnetic fields at the 
end of inflation as 

The resultant spectrum for c = 3/2 is schematically depicted in FigJHl As expected, we 
have a flat spectrum on large scales before the backreaction becomes relevant. However, 
once the backreaction becomes important, the spectrum becomes blue. Thus, after taking 
into account the backreaction, we realized that primordial magnetic fields on large scales 
from inflation cannot be expected [112j . 



4.6 Phenomenology in anisotropic inflation 

So far, we have shown the possibility of anisotropic hair in a variety of inflationary models. 
To explore observational consequences of this anisotropy we need to work out anisotropic 
cosmic perturbation theory. Since the background is anisotropically expanding, we cannot 
use the standard cosmological perturbation theory reviewed in section [2l2l This analysis has 
been carried out in litereature |113] . Note that, even in a fixed anisotropic background, the 
behavior of fields is non-trivial |114] . In this subsection, in a direct extension of standard 
analysis review in section [22] and motivated by the Bianchi type-I background isometries, we 
classify perturbations under the 2-dimensional rotational symmetry and obtain the quadratic 
actions for 2-dimensional scalar and vector sectors. 
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4.6.1 Classification of fiuctuations 



In order to gain some intuition let us start with the spatially homogeneous and isotropic 
FLRW Universe. For simplicity, we consider K = fiat space which in the conformal time 
coordinate r its metric is given as 

ds^ = a^(r) [—dr'^ + 5ijdx^dx^~\ . 

In the C = E = 0,Wi = gauge (c/. ( 12. 24^ )1^ the metric perturbations are parameterized 



as 



rfs^ = [-{1 + 2A)dT^ + 2(5,i + Vi)dTdx' + {Sij + hij)dx'dx^] , 



(4.130) 



where we imposed Vi^i = and hijj = ha = 0. As usual, dealing with linear perturbations 
it is convenient to work in the Fourier space. Since there exists 3-dimensional rotation 
symmetry, we can take the wave vector to be k = {k, 0, 0). Then, the perturbed metric has 
the following form 



( 



\ 



* 
* 




* 





-a'^h+ J 



* is symmetric part. 



(4.131) 



Here, we utilized the special choice k = (A;, 0, 0) to solve the constraints Vi^i = and hij j 



-h. 



ha = 0. With the same token only B .^ remains, and hyz 

To prepare the setup for studying the cases with only 2-dimensional rotation symmetry 
in y — 2;-plane we consider k = {k^, ky, 0). This case can be obtained from the previous case 
upon the rotation 




1 

k 



ky 





k,, 



kx 








)l 









(4.132) 



where k^ = k^ 



kl- Under this rotation, the perturbed metric becomes 



/ -2a^A 



^9, 



\ 



* 



* 



—a?h+ J 



(4.133) 



To simplify the perturbations, we use the remaining part of gauge transformations induced 
by spatial dependent diffeomorphisms 

u;fj, ; 



13 



If we ignore vector and tensor perturbations Vi , hij , this gauge is called flat slicing gauge. 
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where 



we obtain 



flU 



\ 



-2a^A 
* 
* 



a?hA 



* 

* 



aV3 + --- \ 



,2_fc 



y 





(4.134) 



where we have omitted some unimportant parts. It should be noted that the above gauge 
transformation does not change flat slicing gauge for which the 3-dimensional scalar curvature 
vanishes. 



In the anisotropic models with the background metric 

dsl = a{Tf{-dT^ + dx^) + h{Tf{dy^ 



dz' 



(4.135) 



that is, a = e" 6 = e""'"°', dr = dt/a, o the isometries are smaller than the previous case. 
Nonetheless one can show that, as we have done in fl4.134p . in appropriate gauge the most 



\ 



brou^ 


^ht to the form 






Sgox 


Sgoy 
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SQxx 





5gxz 
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* 


Sgyy 
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* 


* 


5gzz 



(4.136) 



where we can impose further conditions so that the perturbed metric goes back to (14.1341) 
in the isotropic limit. 

One can classify the perturbed metric using the rotational symmetry in y — 2;-plane. 
In 2-dimensional fiat space, an arbitrary vector m" where a = y,z can be decomposed 
into the scalar part m"^ ^ and the vector part m"^ = 0. In 2-dimensions the tensor 
can be constructed from 2-dimensional vector and scalar. Thus, the most general metric 
perturbations can be classified into the scalar sector and the divergence-free vector sector: 
There are three 2-dimensional vectors and seven 2-dimensional scalars. There are, however, 
three scalar and one vector gauge transformations and hence we deal with two gauge invariant 
vectors and four gauge invariant scalar combinations. Thanks to the symmetry in the y — z 
plane, without loss of generality, we can take the wavenumber vector to be k = {kx,ky,0). 
Hence, the vector sector in 2-dimensional classification can be identified as Sgoz, ^gxz in the 
above perturbed metric. The remaining components Sg^Q, Sgo^, Sgoy, ^gxx, ^gyy: ^gzz belong 
to the scalar sector. (Note that out of these six scalars, only four are independent.) 

In order to calculate the statistical properties of primordial fluctuations from anisotropic 
inflation |115[ I116[ 1117] , we need to reduce the action to the one for physical variables. 
Then, we can quantize the system and specify the vacuum state. We note that besides metric 



14 



Notice that the conformal time in anisotropic inflation is the conformal time in 2-dimensional part (r, x) . 
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perturbations, one has to consider gauge field(s) and inflaton scalar field perturbations too. 
As in the 3- dimensional case and due to rotation symmetry, the scalar and vector parts do 
not mix at linear level and we can hence analyze the vector and scalar sectors, separately. 

In what follows we start with the anisotropic inflationary model introduced and discussed 
in section 14.21 



2d vector perturbations. The metric perturbations that belong to 2d vector sector can 
be written as 

/ b^p, \ 
b'^T 

* 

* * 
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\ 



(4.137) 



where we have incorporated the anisotropy while keeping the spatial scalar curvature to be 
zero, and the 2d vector part of gauge field perturbations can be taken as 



^^vcctor 



(0,0, 0,D) 



(4.138) 



Note that we have no residual gauge transformation and, in particular, D is a gauge invariant 
under Abelian gauge transformations. And, as we have seen in fl4.134p . F corresponds to the 
cross-mode polarization of gravitational waves in the isotropic limit a = b. Using this gauge, 
we can calculate the quadratic action as 



S 



vector 



drdPx 



b'^ fc2 f2 n2 



■^/2 



Aa? 



2a2 

2^2 



P Pv'P 



(4.139) 



Since the perturbed shift function P^, does not have a time derivative, it is not dynamical 
and its equation of motion is a constraint. There are two propagating dynamical degrees 
of freedom F and D in this 2-dimensional vector sector. Eliminating the non-dynamical 
variable (S^ and defining canonically normalized variables 



b\ky[^ 

V2k 



D = fD 



(4.140) 



we obtain the reduced action for physical variables 



S 



vector 



, 1 ff 



4\D 



+ 



1 fv a ky 
y/2 a b k 



f 
T'D* 



f2'2 p 



fc2 



\D\ 



(4.141) 



F* D 



{k/b)' 
[k/b) 



(YD* + T*D) 
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where k is time dependent and given by 



+ (4-142) 

and becomes constant in the isotropic hmit a = b. In the isotropic hmit a = b, T and D 
represent the cross-mode of gravitational waves and vector waves, respectively. The third 
line in the action fl4.14ip describes how these modes interact with each other. 

Next, we use the slow roll approximation to simplify the action. To this end we start 
with f OT]) 

where 

/ = ^— ^ , (4.144) 

c 

and integrate the above assuming the slow-roll approximation e'/e <^ a /a. The resultant 
expressions are 

a = (-r)^l-^ b = {-ry^-'-^' . (4.145) 

Note that the range (1, oo) for c corresponds to (0, 1) for /. Recalling f l4.13p and fl4.15p we 
have 

f2 '2 

^ = 3(-r)-2/e. (4.146) 
From (14. 7p . the background equation for the vector can be found as 



, (4.147) 



and hence 



f' 

J = i-r)-' [-2-3e + r]- 2Ie] , (4.148) 
where ?7 is a slow-roll parameter (12. lip 

- = 2^ (2e -v) = 2(2e - r/)(-r)-^ . (4.149) 



Furthermore, we obtain 

f" 

- (-r)"2 [2 + 9e - 3r/ + 6Je] . (4.150) 



/ 

Note that all the above equalities are in the first order in slow-roll approximation. 



64 



Substituting these results into the action, we obtain the action in the slow roll approxi- 
mation |117] 



S 



vector 



drd'k ^|f'|2 + i [-P + (-r)-2{2 + 3e + 3/e + 3/esin2^}] |fp 
+ ^\D'\^ + ^ [-k^ + {-ry^ {2 + 9e - St] + eiesin^ e}] \D\ 



+ ^-— (-r)"^ sin^(rD* + T*'D) 



{-ry^ sin 9 {T D* + T*D) 



where we have defined 



sin 6' 



kyd 



(4.i5i: 



(4.152) 



This 6 represents the direction dependence. In the isotropic limit / = 0, the Lagrangian for 
r becomes the familiar one for gravitational waves in an FLRW Universe, cf. discussion of 
section 



2d scalar perturbations. In a similar way, we can derive the quadratic action for physical 
variables in the 2-dimensional scalar sector for the model discussed in section 14.21 The four 
gauge invariant metric perturbations may be parameterized as 



scalar 



( 



\ 



-2a2$ a/3i 
* 2a^G 



a/32 


2\?G 
* 



\ 




-2\?G ) 



(4.153) 



where we have kept the spatial scalar curvature vanishing. The perturbations of the scalar 
(infiaton field 0) will be represented by 50. The variable G and 50 are the gauge invariant 
variables that correspond to the plus mode of gravitational waves and the scalar pertur- 
bations, respectively, in the isotropic limit a = b. The gauge field will also have two 
2-dimensional scalar perturbations 



SA 



scalar 



{SAo,0,J ,0) 



(4.154) 



where we have fixed the Abelian gauge by putting the longitudinal component to be zero. 
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S' 



Prom these ansatz, we can calculate the quadratic action as 



scalar 



d xdr 



52 11 u2 1 

Qj ~r ^ Gi Qj ^ 

fa' h'\ fa' h'\ P 1 



1 P2 ,2 , ,2^/2 .2^2 ,2^2 , 1 ^2^^/2 _ ^^^2 _ b 2 _ 1 2j,2 

2 ^ 2 2 '"^ 2 



2a2 



(/0 + fU<p) ^'"^<P' - a'b^V^' + ^f^"G^ - 2a%'V^G - 2bb'<!>'G 



(4.155) 



^scalar gQvems the dynamics of seven scalar perturbations $, /32, G, 5Ao, 50 and J. Among 
these, $, /?!, /32 and (JAq are non-dynamical and can be eliminated. In the slow roll approxi- 
mation, we obtain the following reduced action [117] 



scalar 



where diagonal parts are given by 



= ^|GT + ^[-A;2 + (-r)-2{2 + 3e + 3/e + 3/esin2^}]|Gp 



{2 + 9e-37/-6/esin2^}] |Jp 
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(4.157) 
(4.158) 

m\ (4.159) 
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^ (G50* + G*5(/>) 



and the non-diagonal "interaction" parts reads 



Here, we defined canonical variables 



G = V2bG , J 



50* J + 50 J* j - (-r)-2 (50V + 
-r)"^ (G*'J+G'J*) - (-r)-' (GV + GJ*) 



■ J , 50 = 650 . 



(4.160) 
(4.161) 
(4.162) 

(4.163) 



66 



Note that G , J and dcf) represent the gravitational waves, the vector waves, and the scalar 
perturbations, respectively. The above action shows there exist the interaction among these 
variables. We notice the scalar part (14.159^ contains anisotropic factor / ( 14.143^ without 
suppression by a slow-roll parameter e. Therefore, to obtain the quasi-scale invariant spec- 
trum of curvature perturbation, I itself has to be small. Note also that the interaction terms 
all vanish in the isotropic 7 = limit. 



4.6.2 Anisotropic coupling vs. anisotropic expansion 

We are now in a position to calculate corrections to power spectrum of various variables due 
to the anisotropy |115l 11161 1117j . From the actions fl4.15ip and (14.1561) . we see there are 
two sources of statistical anisotropy for fluctuations. The first comes from the anisotropic 
expansion itself and is encoded in f l4.157p . f l4.158p . and f l4.159p while the second comes from 
the couplings f l4.160p . f l4.16ip and f l4.162p due to the background vector field. As we discuss 
below, the latter dominate over the former source. 

The first source of statistical anisotropy in fluctuations can be understood intuitively, 
recalling that the amplitude of canonical perturbations at horizon crossing is equal to the 
effective Hawking temperature i7eff/27r, where ifeff is the effective expansion rate. In the 
anisotropic case this i^cff depends on direction and hence we expect to see a direction de- 
pendent power spectrum. The essential structure of the statistical anisotropy of fluctuations 
due to the couplings can also be understood without complicated calculations, recalling the 
gauge field-inflaton coupling term 

and the order of magnitude of background quantities 

Pv'^ . U V 1 

For example, to obtain the J-G coupling, one of the -F^i^'s has to be replaced by the back- 
ground quantity v'. Hence, the coefficients in the J-G coupling is proportional to fv' which 
is of the order of \/Te. This explains the strength of the coupling in f l4.162p . Similarly, J-5(p 
coupling is proportional to f^v' because we have to take the variation with respect to 0. 
Hence, we can estimate its magnitude to be a/7. This explains the interaction term f l4.16ip . 
Finally, the coupling G-6(j) has a magnitude of the order of /^f'^ which is proportional to 
I^/£. This shows a good agreement with the coupling fl4.160p . Thus, we can understand why 
there is a hierarchy among the couplings of the gravitational waves, the vector waves and 
the scalar field. 

To proceed further, we should quantize this system by promoting canonical variables 
to operators which satisfy the canonical commutation relations, with appropriately chosen 
vacuum state. In the deep subhorizon limit kr -> — oo, the actions fl4.15ip and f l4.156p reduce 
to those of independent harmonic oscillators where we choose the Bunch-Davis vacuum state 
|0) by imposing the conditions aa^k|0) = at an initial time Tj. Here, a^^k is an annihilation 
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operator whose commutation relations are given by 



We are interested in the power spectrum of the scalar perturbations 

(O|50k(r)50p(r)|O) = P5^(k)5(k + p) , 

and the power spectrum of the cross and plus mode of gravitational waves 

(0|fk(r)fp(r)|0) = Pr(k)(5(k + p) , 
(0 Gk(r)Gp(r)|0) = PG(k)5(k + p) . 



(4.164) 



(4.165) 



(4.166) 
(4.167) 



We can also calculate the cross correlation between the plus mode of gravitational waves and 
the scalar perturbations 



(O|50k(r)Gp(r)|O) = PsM^)5i^ + P) ■ 



(4.168) 



As is mentioned in the previous subsection, / has to be small. And, the anisotropy in the 
expansion rate is much more small. Hence, we can treat the anisotropy perturbatively and 
estimate its magnitude by using perturbation in the interaction picture. In the interaction 
picture, the expectation value for a physical quantity 0(r) is given by 



(m|0(r)|zn) 







Texp [ i Hj{t )dT 



0(r) 



Texp 



Hi{T)dT 







(4.169) 

where \in) is an in vacuum in the interaction picture, T and T denote a time-ordered and an 
anti-time-ordered product and Hi denotes the interaction Hamiltonian. This is equivalent 
to the following 



{in \0{t) \ in) 



■N 



N-i ■ ■ ■ 



T2 



dri 



X (0 |[i/,(n), [Hi{r^), ■ ■ ■ [Hj{rM). 0(r)] ■ ■ ■]]! 0) . 



(4.170) 



In our analysis, we assume the noninteracting part of Hamiltonian to be that of free fields 
in de Sitter spacetime, which for each mode k it is 



E 



where 



Qn,\,{r) 



^i(i")a„„k + M(r)*ajj _k, 



u{t) = 



-ikr 



2k 



1 - 



kr 



(4.171) 

(4.172) 
(4.173) 
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where Qn represent the physical variables Z), F, G, J, 50. The rest of the Lagrangian (14.1510 - 
(14.1621) may then be regarded as the interaction part —Hi = L*^^) — Lq. To see the leading 
effect on the anisotropy in the scalar perturbation, which is of the order of /, we evaluate 
the correction due to the interaction given by 



<t>j 
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sin 6 



(-r)-i (5yj + 5'(p'.p) + (-r)-2 f50^J + 50jtjl .(4.174) 



From (14.1591) and in the analogy with the slow-roll parameter in the ordinary slow-roll 
inflation, we find that the mass term for 6(j) is proportional to 



/ sin^ e6(l)6(j)^ 
and is expected to give the anisotropy 

6 {in |50k50p| in) 



sin^ eiN{k) 



(4.175) 



(4.176) 



where N{k) is the number of e-folds from the horizon exit of fluctuations with wavenumber 
k to the end of the inflation. We will show below that the corrections coming from the 
interaction term Uf'^ is proportional to N{kY and hence dominate over the correction due 
to the mass term. Thus, the leading correction is given by 



5{in |50k(T)50p(T) I in) 



dri (0 







(4.177) 



Using (I4.172P and commutation relations for the creation and annihilation operators, we 
obtain the anisotropy expressed as follows 



5 {in 
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sin^ 6 I dv 
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UiT 



-Im 



-T2) U [T2)U [T) 



-T2) U[T2)U [T) 



U[Ti)U [T2) 



(4.178) 



where Im denotes the imaginary part. Substituting the function form of u (I4.173P and 
introducing time variables x = kr, xi = kri and X2 = kT2, we have 



6 {in |50k50p| in) 
(O|50k50p|O) 
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sm tl 
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(4.179) 



The contribution to the integral from the subhorizon — xi 3> 1 is neghgible. In the hmit of 
superhorizon — Xi <^ 1, we also have — X2 Ij^X 1- Hence, the integrand in fl4.179p 
approximately becomes 8/xiX2- Thus, the anisotropy can be evaluated as |117j 



5 {in 




in) 




50k50p|O) 
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ix) = 7 sin e dx2 dxi 



X rX2 



1 - ^ J-i i-i X1X2 

24/ 

= ——sin^9N\k), (4.180) 

where N{k) = — ln(— /cr) is the number of e-folds after the horizon exist (to the end of 
inflation). 

Similarly one can compute anisotropy in both polarizations of gravitational waves |117] 

6Iesin^ e N^{k) , (4.181) 
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where we used the interaction term in the action fl4.15ip for f and that in f l4.162p for G. It 
is interesting to calculate the cross correlation. The leading contribution comes from Hj'^ 
and Hf"^. The result is as follows |117] : 



{in 


l^'/'kGp 


in) 




50k50p|O) 



-24/^^ — jN\k) . (4.182) 

As we will soon show, this might be detectable. We remind the reader that in these results 
we have ignored the anisotropy of the background expansion while considered the anisotropic 
coupling due to the existence of the background gauge field. 



4.6.3 Statistical anisotropy from anisotropic inflation 

We found anisotropic inflation is an attractor in supergravity motivated inflaton-gauge field 
coupling given in ( 14. ip with a wide range of gauge kinetic functions. The metric during 
inflation approximately reads 

ds^ = -de + e^^* [e-^'^'dx^ + e^^* [dy^ + dz^)] , (4.183) 

where H and S describe the average expansion rate and the anisotropic expansion rate, 
respectively and their values are specified by f l4.143p . Remarkably, the degree of anisotropy 
is at most of the order of the slow-roll parameter e. The point is the existence of anisotropic 
coupling during inflation. 

We can now discuss cosmological implication of an anisotropic inflationary scenario. 
There are many interesting phenomenology in anisotropic inflation. We start with the 
anisotropy in the power spectrum we is parameterized as 

P(k) = P{k) [I + sin^ 6] , (4.184) 

where will be replaced by Qs for (scalar) curvature perturbations and gt for tensor pertur- 
bations. Then, we can predict the following 
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• There exists statistical anisotropy in curvature perturbations of the order 

Qs = 2AIN^{k) . (4.185) 

In [116j . it is pointed out that the sign of predicted by our models is different from the 
observed one [ID] . Note that our parametrization in fl4.184p is different from that used 
in |15] . The positive sign in our definition means the anisotropy is oblate, on the other 
hand, the negative sign means the anisotropy is prolate. In |10], they found prolate 
anisotropy. However, it is possible to modify the model so that the sign of Qs is flipped. 
As reviewed in subsection 14.4.11 the dynamics of vector fields tends to minimize the 
anisotropy in the expansion of the universe and leading to the orthogonal dyad |1U6] . 
Then, the orthogonal direction to the plane determined by two vectors becomes a 
preferred direction. In this case, the sign of Qs becomes opposite. The reason is that 
the vector field tend to make the expansion slow along the direction of the vector. 
Hence, for a single vector, we had oblate anisotropy. While, the same mechanism gives 
rise to prolate anisotropy for the dyad. We can also utilize anti-symmetric tensor fields 
|104] to achieve the same goal. 

• There exists statistical anisotropy in gravitational waves of the order 

= 6IeN^{k) . (4.186) 

• These exists the cross correlation between scalar perturbations and gravitational waves 
of the order of — 24/y^iV^(/c). Using the definition of curvature perturbations "Jlc = 
d(j)/\/2e, one can translate the cross correlation (14.1821) between the scalar perturba- 
tions and gravitational waves to that between the curvature perturbations and gravi- 
tational waves normalized by the power spectrum of curvature perturbations: 



{in 

Tr. = 



, , ■ ;\ ^ , \ = -24V2IN\k)e . 4.187 



Here, we should notice that there appears an enhancement factor N'^{k) in the above formula. 
This is because the interaction on superhorizon scales persists after horizon crossing during 
inflation. Because of this enhancement, even when the anisotropy of the spacetime is quite 
small, say S/if ~ 10"'' in our example, the statistical anisotropy imprinted in primordial 
fluctuations can not be negligible in precision cosmology. In other words, we have to assume 
/ ^ 1. In fact, from the observational upper bound Qs < 0.3, there is a cosmological 
constraint 

where we used the result in [16] . Since / is derived from the gauge kinetic function (I4.144p and 
the e-fold number N{k) can be determined once reheating process is clarified, the constraint 
on Qs implies the constraint on the gauge kinetic function. 

After taking into account the observational constraint, the anisotropy in the background 
metric of anisotropic inflation is negligibly small. Nonetheless, as we discussed the anisotropic 
effects may arise through perturbations which can have observational prospects. 
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Figure 9: The off diagonal TB and EB spectra C/^^i, Cig^i induced by the cross correlation. 
As a reference, the conventional diagonal BB spectrum induced by isotropic part of the tensor 
perturbations is plotted with a dotted line. The parameters gg = 0.3, r = 0.3 are taken. 



It is useful to notice that there exist consistency relations between observables 

Agt = egs, r, = -V2e g, , r = 16e . (4.189) 

The consistency relations allows us to test anisotropic inflation in a model independent way. 
Let us explain how to use it. It is known that the current observational limit of the statistical 
anisotropy for the curvature perturbations is given by gg < 0.3 [IB]. Note that, according 
to [IS], a signal as small as 2% can be detected with the PLANCK. Now, suppose that we 
detected gg = 0.3. We also assumed the tensor-to-scalar ratio to be r = 0.3. Note that 
this particular number is not important, as is explained later, as long as we can detect 
gravitational waves. Then, the consistency relations would give us predictions. Namely, 
anisotropic inflation implies the anisotropy in the gravitational waves 

gt - 10~2 (4.190) 

and the cross correlation 

= -V2gse ~ -4 x 10"^ , (4.191) 

where we used gg ~ 0.3 and e ~ 10"^. If these predictions are confirmed by the CMB 
observations, that must be a strong evidence of anisotropic inflation. 

We showed in the above how the consistency relations are used to make predictions. So, 
the next task is to clarify how we see these features in the CMB. The answer is that the 
anisotropy induces off-diagonal TB,EB spectra C^^^+i. Here, we define the angular power 
spectrum C^^/ = {ae^rnCie',m') with coefficients aim of the spherical harmonic expansion. In 
Fig. IHl we have depicted examples of TB and EB correlations Cj^g^^, Cf^^^ |118] with the 
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parameters r = 0.3 and Qs = 0.3. As a reference, the conventional BB diagonal spectrum 
induced by the isotropic part of tensor perturbations is also plotted with a dotted line. Note 
that unlike parity violating cases for which odd parity correlations Cj^, exist [TTl [841 
l85l 186] , our model predicts even parity correlations such as C^i^^ as the result of parity 
symmetry of the system. The ratio of TB correlation induced by cross correlation to the 
isotropic BB correlation is not dependent on e (or r) for a fixed value of Qs in our anisotropic 
inflation model. For the optimistic value of Qs ~ 0.3, both amplitudes become comparable. 
This simple order estimation implies that the TB signal could be comparable to that of B 
mode correlation induced by primordial gravitational waves. Hence, anisotropic inflation 
can be a potential source of off-diagonal TB correlation. Since the current constraints on 
the TB/TE ratio is of the order of 10~^ [20], we need to improve the accuracy by one more 
order, which might be achieved by the PLANCK. 

It should be stressed that the statistical anisotropy is a natural consequence of the fact 
that gauge fields are ubiquitous and there remains many other phenomenological aspects 
inflation models with gauge field hair to be explored. 



4.6.4 Non-Gaussianity in anisotropic inflation 

As we explained in section[3], we know the vector field can induce a non-trivial non-Gaussianity 
in the presence of a gauge kinetic function. Since the anisotropy in the expansion is quite 
small in realistic cases and the dominant effect is a coupling between scalar and vector, one 
can calculate non-Gaussianity in the isotropic universe by neglecting the metric backreac- 
tion. Concrete calculations in this regards has been carried out in |120j with the statistical 
anisotropy in the power spectrum 

I^.ICMB > 0.1 — . (4.192) 

Moreover, they have evaluated the non-Gaussianity as 

/■ nr \9s\cMB , . 

/nl ^ 26 ^ ^ . (4.193) 

This is another consistency relation which anisotropic inflation predicts. Namely, there exists 
a large non-Gaussianity if the statistical anisotropy is also large l^f^l ~ 0.1. 

The authors in |120] pointed out quantum fluctuations of the electro-magnetic field are 
comparable to the classical background when c— 1 is severely constrained by the observations. 
For these cases, quantum diffusion process would be important and those affect the classical 
background in the long wavelength limit. Nevertheless, they succeeded in making a statistical 
prediction for the non-gaussianity. Namely, the consistency relations are robust and the 
models are testable. 

^^There are other mechanism to produce the statistical anisotropy |119j . However, the consistency relations 
can be used to discriminate anisotropic inflation from others. 
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5 Gauge-flation: Inflation driven by non-Abelian gauge fields 

As mentioned in the previous sections the energy density which drives the slow-roll infla- 
tion dynamics is generically provided by one or more scalar fields and that scalar fields 
are the natural choices if we are interested in homogeneous and isotropic FLRW models. 
Nonetheless, as we will review in this section, it is possible to drive isotropic infiation by 
non-scalar fields and in particular vector gauge fields. This mo del/ scenario will hence be 
called gauge-fiation \121\ 1122] . In this section we spell out how it is possible to get quasi- 
de Sitter expansion from a gauge field theory minimally coupled to Einstein gravity and 
study the stability of the infiationary trajectories against classical isotropic and anisotropic 
perturbations |123] . We then study gauge-fiation cosmic perturbation theory and compute 
the information which could be extracted from our second order (quadratic) perturbation 
analysis and finally compare these results with the cosmological observation data. 



5.1 The theoretical setup, how to get an isotropic model 

Our starting point in gauge-fiation is a non-Abelian gauge theory which is minimally coupled 
to Einstein gravity in four dimensions. The gauge group can be any non-Abelian compact 
group G which could be simple or a generic product of simple groups. Generators of algebra 
of G will be denoted by T^, A = 1, 2, ■ ■ ■ dimG. Lagrangian of the models of our interest is 
of the generic form 

L = L{g^,, F^,) = \r + ^g{F^u) , (5.1) 
F^u is the field strength of the gauge fields A^: 

F^u = dX - - af'^'^A^^A^^ , (5.2) 

where g is the gauge coupling and f'^^'" are the structure constants of the gauge group G, 
\T^,T^] = if^^'^T'^. We choose the standard normalization Tr(T^T^) = ^S^^. Under a 
generic local gauge transformation U E G 

A. A'=-iU-^dM + U-^AM, 

^ r 1 (5-3) 

F s. F = U F U 

and £jG{.F^y) is a generic gauge invariant action made out of powers of F^y whose Lorentz 
indices are properly summed over using only metric g^j_i, and/or the Levi-Civita tensor e"^/^'^ 
and the gauge indices are summed over by taking trace. That is, JZg is a minimally coupled 
diff-|-gauge invariant action. Moreover, with our choice [Lq being only a function of F^jjy) 
gauge field equations of motion will be at most of second order in time derivatives. In 
particular, we stress that may include Yang-Mills term —^TrF^^ but is not limited to 
that. 

We would like to get an infiationary background driven by gauge fields as infiaton. If 
we ignore the gauge indices for the moment, i.e. considering an Abelian gauge field, turning 
on a vector gauge field in the background will break rotational symmetry of the fiat FLRW 
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geometry, unless only the temporal component is turned on and homogeneity implies that this 
component should be only time dependent. This gauge field will hence have a vanishing field 
strength, it is a pure gauge. Now, let us consider non-Abelian gauge fields. Considering gauge 
indices will change the situation: Gauge fields are defined up to gauge transformations and 
two gauge fields which are related by gauge transformations (15.31) are physically equivalent. 
To be precise, let us fix the background metric coordinates to the standard spatially fiat 
[K = 0) FLRW one (12.21) . In this coordinate system we may turn on a generic vector 
gauge field with spatial components Ai and homogeneity implies Ai = Ai{t). Note that Aq 
may always be set to zero in a particular gauge, the temporal gauge; i.e. one can always 
find U = U{t) for which A'^ = 0. This choice, therefore, fixes the gauge transformations 
up to time independent gauge transformations]^ We have hence fixed the gauge freedom 
up to space-time independent global gauge transformations. Therefore, up to global gauge 
transformations, our background gauge field is (note that for the ease of notation we have 
suppressed the explicit gauge indices) 

Ao = 0, A = A{t). (5.4) 

As discussed, turning on a vector gauge field A^ will break rotational invariance. We 
use the remaining global gauge transformations to remedy this: Two fields related by gauge 
transformations are physically equivalent, i.e. 

(Ai)G = U'^AiU =G Ai, with constant U eG . (5.5) 

The above gauge transformation acts upon the gauge indices which are suppressed in the 
above equation. On the other hand, upon global spatial rotations 

Ai {Ai)R = R,jAj . (5.6) 

If the background configuration Ai is chosen appropriately such that Ar = Ac then our 
physical gauge configuration will preserve rotational symmetry. In other words, rotational 
non-invariance caused by turning on vector fields in the background may be compensated 
by the global gauge transformation. To see this, we first note that Ai is in vector (triplet) 
representation of the rotation group S0{3)r. On the other hand, any non-Abelian gauge 
group has an SU (2) subgroup. This SU (2) may be identified with S0{3)r for some particular 
gauge field configurations. As far as our current discussion is concerned, without loss of 
generality, we can choose the gauge group G to be SU{2) or 5*0(3) and choose the T'^'s 
to be SU{2) generators in the triplet (adjoint) representation icr", a = 1,2,3 where a" are 
Pauli matrices. Let us choose the background to be 

A^ = a{t)m5-, 

^^We may set dimG number of relations between gauge fields and their first derivatives through gauge 
choice, and in the temporal gauge Aq = we have used all dimG choices. Note also that one may consider a 
generic space-dependent gauge transformations U = U{xi), will not change Aq = 0, however, will introduce 
x' dependence into the spatial components Ai and hence we will not consider these gauge transformations. 
Therefore, going to temporal gauge completely fixes the local gauge freedom. 
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where a{t) is the scale factor and ilj{t) is a scalar under rotations. This latter may be seen 
immediately, noting that one can rewrite the above as the following form 

At = me%, (5.7) 

where are the triads of the spatial hyper- surface. It is readily seen that for a gauge field 
of the form (15.71) . the effects of any generic rotation by angle 9, R = e^°"'^, can be viewed 
as gauge transformation U = e^"'^ with A = —9. To summarize, we have shown that the 
ansatz 

" { a{t)i^{t)5'i ^ = z ^^-^^ 
leads to a rotationally invariant, homogeneous background. 

The "homogeneous-isotropic" non-Abelian gauge field configurations in the context of 
FLRW cosmology (and not necessarily inflationary models) has been discussed in the liter- 
ature before, e.g. see |124[ 1125] ,^ Here we discussed the non-Abelian gauge fields which 
as we will see are suitable for building inflationary models. From our discussions above, 
however, one can readily see that it is possible to maintain homogeneity and isotropy with 
non-Abelian global symmetry (i.e. having several Abelian U{\) gauge fields rotating among 
each other by a global symmetry G, which has an S0{?>) subgroup). Such Abelian gauge 
field configurations in the context of cosmology has been called "cosmic triads" |126l [5l [91]. 

In the above, and hereafter, we will only consider SU{2) gauge theory, however, our 
arguments can be directly generalized to an SU (2) subgroup of a generic non-Abelian gauge 
group G. With the above choice fl5.8l) all gauge freedom, local and global, has been employed. 
Note also that, although one can always set = by a gauge choice, taking as we have, 
is our choice, ansatz, for the gauge field. That is, we have chosen to work with a very 
particular configuration of the gauge fields. We hence need to study two issues: 1) The 
choice ( 15. Sp is compatible with the equations of motion of the action ( 15. ip with an FLRW 
ansatz for the metric |121l 1122] and, 2) this choice is stable under both classical |123j and/or 
quantum fluctuations or perturbations. In what follows we will study these two questions. 
The first question will be dealt with in 15.11 and the second in 15.51 and in 15.61 



Consistency of reduction 

Here we show that if we start with a gauge field of the form (15. 8p at, say t = tj, the dynamics 
of the system does not take us out of this "homogeneous and isotropic sector" (where the 
gauge field is of the form (15.81) and metric is FLRW). To this end we should analyze the 
equations of motion of the theory: 

G^^ - - -^7^ I5.9aj 

^.|^ = 0, (5.9b) 



^^For the closed Universe case with K ^ 1 (|2.2I) . there are other choices for the homogeneous isotropic 
gauge field configuration |125) . 
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where is the covariant derivative and is defined as 



for a generic tensor X^ in the adjoint representation of the gauge group G and is the 
diffeomorphism-covariant derivative computed for the FRLW metric (12.21) . Restricting to 
the G = SU{2) case, 



We are now ready to insert the ansatz (15.81) into the above equations. We first note that 
F^, = a{t){^ + Hm , F^j = -9a{ty^\%j , (5.11) 
where dot is derivative w.r.t. comoving time, and that 

T,, = 2^n, + (7^,£G, (5.12) 

in the above Lq = ^g{F^u^ g^^) (the minimal couphng condition) has been used@ Inserting 
fim]) into flHTT^ one finds 

T't = diag{-p,P,P,P), (5.13) 

where 
and 

(j) = a{t)^{t) . (5.15) 

(Note that 0, unhke ■0, is not a scalar.) C^red. is the reduced Lagrangian density, which is 
obtained from calculating L{F^^ \ g^^) for field strengths F^^^^ given in (15. lip and FLRW metric 
(12. 2p |122] . Since the energy momentum tensor is of the form of a homogeneous-isotropic 
perfect fiuid, the metric ansatz will remain of the FLRW form when we evolve in time using 

dSHD. 

To show that the gauge field ansatz (15.81) maintains its form in time we should consider 
(]5.9bp . Plugging the ansatz into the equation of motion (I5.10p one can readily see that the 
1/ = components of the equation of motion are trivially satisfied and the v = i components 
become proportional to (5f , leading to the single equation 

a-^dt d(j) (J(f) 

where Lred.i.4'^ 0! '^(0) is the reduced Lagrangian introduced above. 

Technically, we have shown that there exists a consistent truncation/reduction of the 
gauge field theory to the homogeneous-isotropic sector specified by the scalar field ip (or 
0). Moreover, the Lagrangian governing the dynamics in this sector is simply obtained from 
computing the complete Lagrangian over the reduction ansatz. 



^®For Lq which is a polynomial in F^j^^, one can show that the first term in (|5.12p is symmetric in 
indices. 



77 



5.2 Choosing the gauge theory action, the gauge-flation model 



As discussed regardless of the form of the action one can always reduce the theory to 
a isotropic-homogeneous sector. The question we would like to tackle now is to find the 
appropriate action(s) which can lead to successful slow-roll inflation, and for this purpose 
one may only focus on the reduced Lagrangian L^-ed- To gain some intuition about the form 
of reduced Lagrangian it is useful to work out explicit form of some terms appearing in the 
Lagrangian. The space-time indices may be summed over using metric, e.g. 

pa^pt^u b ^ 25'^^[-(^ + H^f + g^^^] , (5.17) 
or using the Levi-Civita tensor 

^ap^upa^pb^ = -Sgd^'^tp + HiIj)iIj'^ . (5.18) 

As we see Lorentz invariance (summing over all free space-time indices) implies that ip is 
always accompanied by a factor of Htp and, if we use metric for summing over the indices 
ip+Htp term appears together with the same power of ip"^. Therefore, the reduced Lagrangian 
is neither of form of a kinetic term plus a potential, it is nor of a simple K-infiation type 
where the action is only a function of X = {d^ipY or powers of ip] our gauge-fiation model, 
even in the homogeneous-isotropic sector, is hence not a special K-inflation model. Due to 
this particular features, one should work out the slow-roll conditions from the first principles, 
e.g. to have accelerated expansion p + 3P < and 

As a warmup example, let us consider the Yang-Mills action, for which 
1„ ^„.,. 3, 



PYM = \[{i^ + Hijf + g'ij'] , Pym = Ipym = l[ii^ + H^py + g'^'] 



(5.20) 



As expected for the Yang-Mills theory T^^^, is traceless, yielding p = 3P. For this case, 
p + 3P = 2p > and the theory does not lead to accelerated expansion. Noting the features 
discussed above, it turns out that satisfying p + 3P < and p > conditions is not trivially 
obtained for actions which only functions of Fp^^F^^'^ ^ (such terms in the context of cosmology 
has been considered e.g. in [124^ 1125] 1127] ) and one may consider terms involvin g e" ''^'^'^. The 
simplest choice which is gauge invariant and is not a total derivative is the terno 

{FAFf^^ie-f^'^-'F^^^F^.f. 

The dependence of this metric on metric g^^ is only through det g and is of the form 
1/V— det g. Therefore, the contribution to the energy momentum tensor from the above 



^^Note that A F is a total derivative and a topological term, i.e. it does not contribute to equations of 
motion and the energy-momentum tensor. 
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F'^ term will have P = —p, making it particularly interesting for inflationary model building. 
We therefore choose the gauge-fiation action to be of the form |121l I122j 



S = / d^Xy/^ 



PI f.- 

pa p fJ.u , / fluXcT -pa pa \2 



{5.21] 



with the reduced gauge field Lagrangian 

^red = l[{^ + Hijf - + t^g^i^ + , (5.22) 

and ^ 

P = pYM + Pk, P = ^Pym - pK , (5.23) 

where pyu is given in fl5.20l) and 

p^ = %,g\^^E^f^\ (5.24) 

In the above action k is a dimensionful parameter, of mass dimension —4. The Lagrangian 
f l5.2ip has then two parameters, a dimensionless Yang-Mills coupling g and a dimensionful 
coupling K. We also note that k always appears in combination Kg^ . Demanding for any 
field configuration the energy density is positive definite (weak energy condition) we take 
K > 0. 

Before entering into the details of inflationary trajectory analysis we would like to make 
some comments about the above choice of the gauge theory action. One may ask if it 
is possible to write the reduced Lagrangian in as a diff-invariant Lagrangian for a scalar 
field 'il'{t,Xi), such that when we restrict to ip = ■ip{t) we obtain Lred- Let us consider the 
Yang-Mills term first: 

Sred-YM = \ j dta^i^ + H ^jf - g' = ^-j dta'[i^^-g^'-{2H^+H)ij'] + ^- j d{a'Hij') . 

(5.25) 

Recalling that the Ricci scalar R for the FLRW geometry is i? = Q{2H^ + H), dropping the 
total derivative term the above action may be written in a diff-invariant form 

S^M-red = d'x[-{d,ij)' - \rH^ - g^'\ , (5.26) 

which is a A(/)^ theory with coupling g^ and a conformal mass term. This is expected, 
because the Yang-Mills theory is at classical tree level a scale invariant theory. Finding the 
diff-invariant form of K-term is more involved: 

Sred^^ = hg' j dta\^ + H^f^' = hg' j dt{a'[^'^' + ^j' - ^^'^^^] + tot.der.) 

(5.27) 

The first two terms have an obvious diff-invariant extension of the form {d^ipYil)^ and ip^ 
form, the last term which is proportional to 3if ^ + H, does not have a simple form in terms 
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of curvature invariants!^ Although it is possible to write the reduced gauge-flation action 
in a diff-invariant form for scalar field ip, this action is not in a minimally coupled form (e.g. 
note the conformal mass term), and the K-term takes a contrived, practically useless form. 
We should finally point out that the equivalence of the diff-invariant reduced action and 
gauge-flation action (15.211) on the homogeneous-isotropic trajectories is only true at classical 
level; quantum mechanically and once we consider quantum fluctuations this equivalence 
does not persist. Taking note of this point will be crucial in studying gauge-flation cosmic 
perturbation theory, to which we will return in section 15.61 



5.2.1 Theoretical motivation for the gauge-flation 

As discussed the gauge-flation action ( ]5.2ip . and in particular the K-term, was primarily 
chosen in search for a model with quasi-de Sitter expansion behavior, i.e. P ^ —p. However, 
one may question the presence and naturalness of K-term from the gauge theory viewpoint. 
This term is a specific term. and higher power of F generically appear in the gauge 
theory (Wilsonian) effective action in one or hi ghe r loops and their effects in the cosmological 
contexts have been discussed in [1251 [l27l [128] O These terms typically appear in the gauge 
theory effective action below the charged fermion mass scale, once we integrated out massive 
fermions. Such terms for QED has been extensively discussed in the literature, in one and 
two loops levels. For example, for QED below the electron mass scale, at one loop level the 
effective action is 

^QED-ane loop = -\f^ + Oi{F^f + Z?^' • (5-28) 

(At level there are only two ways to contract Lorentz indices of photon field strength 
field F^y using only metric, [F'^Y ^^^d P^- Therefore, by gauge and Lorentz invariance only 
these two terms are expected.) The coefficients a and /3 may be computed using standard 
field theory techniques and their values for QED may be found in |129j . What is important 
in our discussion here, is however, the parametric dependence of a, /3; a, /3 terms are one 
loop effects and hence should be proportional to Og^^ = (f;;)^ where e is the QED coupling. 
Moreover, they are of mass dimension —4 and hence they should both be proportional to 
m~^, where me is the electron mass scale. That is, 

(47r)^mg 

One would expect to be able to repeat the same argument for higher loops, to obtain 

type terms at loop / level with coefficient ai where «/ ~ aK The above effective loop 

expansion is of course valid for small couplings and for the low energy photons {Ephoton ^ f^ie)- 



"^^FoT example, one may rewrite 3H^ + H = ■^R + ^Rq. To contract the Rq indices we need to have 
square of time-derivative of ip, which we do not. The closest diff-invariant action that one can write will be 
something of the form ^ip^R which misses terms proportional to H. 

^^In our setup we only work with Einstein gravity. Extensions of our discussions to cases of modified 
gravity, e.g. f{R) gravity, or when there are terms in the action of the form R TrF^ has been discussed in 

[Eg. 
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After the above brief review of gauge theory loop effects, let us return to the K-term. 
Given its Lorentz index structure, this term cannot arise from fermion loops. Nonetheless, 
the axion-type coupling xF ^ F involves the right type of indices. Explicitly, let us consider 
a Yang- Mills theory coupled to massive axions |130j : 

^a.ion-YM = '^FtF^ " \{d,X? " \m\^ + A^Tr (c'^^^'^F^^jF^,) (5.30) 



where x is a massive pseudoscalar with mass M, / is the axion scalcj, taken to be much 
larger than M and A is the dimensionless axion coupling. For energies below axion mass one 
may integrate out the axion, using standard field theory techniques (for a detailed discussion 
see [13H 1132] ). This leads to and action of the form (I5.2ip with 

^ = -\^ / = fM^ (5.31) 

where we have presented the above in terms of mass scale yU which is the cutoff of the axion 
theory. 

Such "massive axions" are commonplace in beyond standard model and string theory 
motivated particle physics models and hence our gauge-flation action is well motivated with 
particle physics models. Of course there remains two important questions 1) whether one 
can ignore all the other loop effects while having a significant contribution from certain loops 
and, 2) whether the range or value of parameters or scales natural to these particle physics 
models allow for a successful model of inflation. Here we will discuss the former and will 
return to the latter in section [61 To make sure that at one loop level the K-term dominates 
over the fermionic loop effects, a, (3 terms, it is enough to check or demand that a^P, 
i.e. 

jjL \ A 



where Mf is the typical fermionic (or generic charged matter) masses in the systemic This 
condition can be met for small enough gauge couplings. Moreover, small gauge coupling will 
also be needed to suppress the higher fermionic loop effects. 

The condition f l5.32p guarantees that the K-term dominates over all the other dimension 
eight or higher contributions coming from gauge field or fermionic loops. However, (slow-roll) 
inflation, as was implicit in our earlier discussions, can take place if the K-term can dominate 
over the Yang-Mills term. The above discussion clarifies this issue too: The /t-term, unlike 
all the other gauge theory loop corrections, comes from integrating out an axion field, which 
is nothing but eliminating the massive axion field by evaluating the action on its classical 
trajectory |131] . Therefore, the /t-term is in fact representing the potential of the axion field 
which can naturally dominate over Yang-Mills term and this all can consistently happen 
within perturbation theory. More discussion on this will be presented in section 16.21 



^^The axion field x takes values in the [0, tt/] range. 

^•^The above equations are obtained for a standard quantum field theory loop analysis on flat background. 
For inflationary background, one should note that the Compton wavelength of fermions 1/Mf which con- 
tribute to the loops should be larger that the Hubble horizon size H~^. 
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For completeness of the discussions, we would like to also mention that it is possible to 
obtain gauge-flation action, naturally, within low energy string theory setup. Explicitly, as 
discussed in |133] , consider compactification of ten dimensional heterotic supergravity theory 
(which arises as the low energy effective theory of heterotic string theory) to four dimensions. 
Consistency of the theory implies presence of certain ol corrections, which include a term 
like the K-term. 



5.3 Gauge-fiation inflationary trajectories, analytic treatment 



The equations of motion f l5.9ap and (]5.9bp . for the flat FLRW metric and gauge field ansatz 



5.8p . in terms of = a{t)ip{t) field, takes the form 



+ ^ + (5.33a) 
.^-(|.^), (5.33b) 
1 + + + «:^)^^ + (1 - 3fi:^)— ^ = . (5.33c) 

We are interested in slow-roll dynamics specified by e, 77 <^ 1. Using the Friedmann equations 
dESSaj), fl5:33bl) and (l539|) we have 



(5.34) 



Pym + Pk 



That is, to have slow-roll the K-term contribution should dominate over the Yang-Mills 
contributions Pym^ Pk ^ Pym during slow-roll period. As we will see the time evolution 
of the system increases p^j^, with respect to pushing inflation to its end. The end of 
accelerated expansion is marked by e = 1 which happens when p^^^ = p^. 

For having slow-roll inflation, however, it is not enough to make sure e ^ 1. For the 
latter, time-variations of e and all the other physical dynamical variables of the problem, 
like rj and the if) field, must also remain small over a reasonably large period in time (to 
result in enough number of e- folds). As discussed in section [2] for simple slow- roll models, 
?7 through equations of motion also measures rolling velocity of the inflaton (I2.13p . In our 
case, however, this relation is modified. It is useful to define 

..-A. (5^35) 

5 is not an independent slow-roll parameter and is related to e and 77 through equations of 
motion: 

e =2 - k/^^(1 - 5f, (5.36a) 
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To have a "standard slow- roll" , e ~ He^ and rj ^ e,we should demand that 5 ~ e^. Explicitly, 
the equations of motion fl5.33ap . (]5.33bp and fl5.33cp admit the solution 

e~ ^2(7 + 1), (5.37a) 

Kg^i;^c^2-e, (5.37c) 
where ~ means equality to first order in slow-roll parameter e and 

7 = -[TT- > or equivalently H ~ = — — — . (5.38) 

i/^ e — 7(7-1-1) 

In the above 7 is a positive parameter which is slowly varying during slow-roll inflation. 

Recalling (I5.35P and that 5 ~ e^, (I5.38P implies that •jH'^ remains almost a constant 
during the slow-roll inflation and hence 

- Ir^tl^ (5.39) 



li + l li 



where e^, 7^ and Hi are the values of these parameters at the beginning of inflation. As 
discussed the (slow- roll) inflation ends when e = 1, where 

Using the above one can compute the number of e-folds Ne 



f 



dH 7, + l 7, + l 



N,= Hdt = - — ~ In . (5.41) 



Before moving to numerical analysis some comments are in order: 



• ?7 <^ 1 implies that our field values ip is sub-Planckian (as 77 ^ ip"^). 

• If 7j is a parameter with order 1-10 values, then (15.380 implies that ~ g^e. 

• Number of e-folds, as in simple single field models, is proportional to l/e^. 

^''^Based on our numerical analysis, there exists a range of initial values and parameters for which starting 
with S/{HS) ~ 0(1), while ipf ^ e 1, after a very short time (less than an e-fold) 6 becomes very small, 
of order e^. That is, the dynamics of our system suppresses ip very fast ii e,T] <^ 1. Therefore, for all the 
inflationary period we may confidently use S ~ 6(7+!)^^ ~ (^ ^ 7])e/6. See section [5.41 for a more detailed 
discussion. 

^^Note that all the dimensionful parameters, i.e. H,ip and k, are measured in units of Mpi; H, ijj have 
dimension of energy while k has dimension of one-over-energy density. 
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• Further analysis of the equations of motion, also confirmed by the numerical analysis, 
reveals regardless of initial values of 6, for the region of parameter space which leads 
to long-enough slow-roll period, 6/H6 remains small and hence (I5.37bll leads to 

±r.llL (5 42) 

The above equation may be integrated to obtain ip^'^ — ipj"^ ~ — i-j/j^"^, i.e. ipj = 
That is '0 field displacement during inflation is of order the field itself and is sub- 
Planckian, and this is independent of the initial value of 7 parameter. 

• The above conclusion (field roaming being of the same order of the field itself) may also 
be reached noting fl5.37cl) ng'^ip^ll — 5)^ = 2 — e. During slow-roll we can drop 5-term, 
then the ratio of ipf (which is computed for = 1) to ipi is obtained as ^ 
This estimation matches with the above up to percent level (-\/3/4 vs. ^1/2). This 
is a good evidence for the validity of slow-roll approximation for the whole period of 
inflation (till e becomes one). In other words, the above indicates that end of inflation 
and exit from slow-roll happens very fast in the last couple of e-folds in the end of 
inflation. This is the expectation which is confirmed by our numerical analysis of the 
next subsection. 

• The gauge-flation model has other trajectories which could be relevant to cosmology 
and are not necessarily within the slow-roll trajectories discussed here. An analytical 
and numerical discussion of these trajectories has been discussed in [134] . 

5.4 Gauge-flation inflationary trajectories, numeric analysis 

In this part we make a more thorough parameter space analysis of the gauge-flation model 
through a numerical analysis. Gauge-flation model has two parameters, the gauge coupling 
g and the coefficient of the term k, the former is dimensionless and the latter is of 
mass dimension —4. The degrees of freedom of the model in the homogeneous-isotropic 
sector are the scalar field ip and the scale factor a{t). Inflationary trajectories are hence 
specified by four initial values of these parameters and their time derivatives. These were 
parameterized by Hi, ipi and 5i (or ifji). (The initial value of the scale factor a(tj) is not a 
physical observable in flat FLRW model, as it can be absorbed into the redefinition of spatial 
coordinates Xi.) The Friedmann equations, however, provide some relations between these 
parameters; assuming slow-roll dynamics these relations are fl5.37ap - fl5.37cl) . As a result each 



inflationary trajectory may be specified by the values of four parameters, (V'i, V"*; 9i i^)- 
what follows we explore the parameter space of our model more thoroughly and present the 
results of the numerical analysis of the equations of motion fl5.33ap . fl5.33bl) and fl5.33cp . for 
three typical sets of values for {ipi,ipi; g, k). 

5.4.1 Generic features of the slow-roll trajectories 

Our numeric analysis shows that it is possible to get slow-roll trajectories with enough 
(A'e > 60) or even arbitrarily large TVe, for a large region of the parameter space. The 
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Figure 10: The classical trajectory for tpi = 0.035, = -10"^°; g = 2.5 x 10-^ k = 1.733 x 10 
These values correspond to a slow-roll trajectory with Hi = 3.4 x 10~^, ji = 6.62, = 
9.3 X 10~^,5j = 8.4 X 10~^. These are the values very close to the range for which the 
gauge-flation is compatible with the current cosmological and CMB data (c/. discussions of 
section ISTTj) . Note that k, Hi and ipi are given in the units of Mpi. 
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slow-roll region in the parameter space has the following features: 

• Small gauge coupling g, g ^ 10~^ — 10~^. This latter is important because we are 
dealing with a non-Abelian gauge theory with a specific F'^ term, the K-term. As was 
discussed in section 15.2.11 smallness of gauge coupling g is crucial for justifying the 
theoretical stability of gauge-flation action f l5.2ip . 

• Naturalness of k. A typical value of n for slow-roll trajectories is k ~ 10^'^ — 10^^ (in 
Planck units). The n parameter may be related to the cutoff scale fi of an axion model 
f l5.3ip . leading to /i < lO^^GeV, which is a natural scale within grand unified theories. 

• Dependence on initial field values. It appears that with the above values for g and k, 
typical value for the field ip and its roaming is of order 10~^ — lO^^Mpi, in accord with 
e ~ ?7 ~ 10~^. 

• Slow-roll trajectories are not so sensitive to the initial value of ip (or 6 parameter). 
That is, getting slow-roll inflation mainly depends on the values of g, k, ip and to a less 
extent on 6. This feature is also shared with the ordinary chaotic models. This latter 
will become more clear in the analysis of section 16.31 

Moreover, within the above range of parameter space, and regardless of their precise 
values, the slow-roll trajectories show the following generic features. These may also be seen 
from our figures and time evolution diagrams. 

• If we start with a small S, S ^ 10~^, the field ip remains almost constant during inflation 
and at the end of inflation falls off very fast and starts oscillating. However, if we 
start with a large 5, 6 ~ 0(1), ip has a very fast single oscillation, reducing 6 to around 
10~^ and then it follows a small-5 trajectory mentioned above. 

• After inflation ends the dominant term in dynamics of the system is the Yang-Mills 
term, the energy momentum of which behaves like a radiation gas@ Therefore, it 
is expected that: 1) the slow-roll e parameter should asymptote to e = 2 (as is also 
indicated by f l5.34p ) and, 2) the system should essentially behave as a theory. (This 
is expected as the Yang-Mills part of the reduced Lagrangian behaves like a g'^ip^-) That 
is |135j . ip and hence other dynamical variables, should oscillate (following a Jacobi- 
cosine function) with a period of the order gipf ~ Hi fl5.22p . Moreover, the amplitude 
of oscillations should decrease in time like 

• The phase diagram of "effective inflaton field ip" during slow-roll is hence like a straight 
line during slow-roll, while becomes very similar to that of a chaotic model [2]. 

^^Recall that for a radiation dominated Universe with P ~ pji^ a{t) ^ t^/^, H = l/(2t) and e = 2 
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Figure 11: The classical trajectory for ^pt = 0.025, = -10"^°; g = 2.507 x 10~^ k = 1.3 x 10 
These values correspond to a slow-roll trajectory with Hi = 3.63 x 10~^, 7^ = 2.98, = 
2.5 X 10"'^, 6i = 1.1 X 10~^. These figures show that it is possible to get arbitrarily large 
numbers of e-folds within the slow-roll phase of our gauge-flation model. 
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5.4.2 More detailed discussion for various parameter choices 

Here we present evolution of the system for three different sets of the four parameter values 
{g , k; ip , ip) . There are four graphs for each set. Three of them depict evolution of effective 
infiaton field ip, number of e-folds N^. and slow- roll parameter e vs. comoving time t (in 
Hubble units H~^). The last one is the phase diagram of the effective infiaton field, for 
presentation purposes we have drawn <p/a{t) vs. 4>/a{t) = ip. These numeric analysis match 
very well with our slow-roll analytical equations. 



5.5 Stability of the isotropic background 

We showed in previous subsections that gauge-fiation can lead to a successful isotropic and 
homogeneous inflationary dynamics. We demonstrated that in spite of having a non-zero 
gauge field value on the background, the isotropy of the background is preserved due to 
the time independent part of gauge transformations (c/. discussions of section [5TT!) . How- 
ever, one may ask if the isotropic gauge-fiation model is a stable setup against the initial 
anisotropics. In other words, due to the gauge- vector nature of infiaton field in gauge-fiation, 
a question that may arise naturally is the classical stability of gauge-fiation against the initial 
anisotropics and choice of initial conditions. In this section, we will investigate this issue 
and in order to study the generality of the isotropic FLRW gauge-fiation, here, we will study 
gauge-fiation in a homogeneous but anisotropic background. This subsection is mainly a 
review of |123] . Here, for practical reasons we consider gauge-fiation in an axially symmetric 
Bianchi type-I setup but we argue that our results is expected to be valid for more general 
anisotropic cases. (For a quick review on Bianchi cosmologies seeO) Similar question for a 
general infiationary setup was posed and analyzed in [8] and reviewed in Appendix [Dl 

Consider Bianchi type-I axially symmetric metric (14. 6 p 

where e'^^*-' represents the anisotropy and e"*-*^ is the isotropic scale factor. Due to the 
symmetries of the metric, as before, the temporal gauge is an appropriate choose for the 
gauge fields A"-q = 0, and we consider the following modification to the isotropic ansatz (15.81) 

A", = e%iJi = c/m^(e°W-2'^(*Vi, e"W+"(*V2, e"W+"(*V2), (5.43) 

where ipi act as two scalar fields and 6""^ are the triads associated with spatial metric. Without 
loss of generality and for some technical reasons which will be clear shortly, we parameterize 
tpi and ip2 in terms of and A as 



^1 = ^2 = ^^ = A^. (5.44) 



Furthermore, as in the isotropic case (I5.15p . the equations take a simpler form once written 
in terms of cf) and a 

<p{t) = a{t)ip{t), a{t) = e". (5.45) 
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Figure 12: The classical trajectory for 51 = 4.004 x 10"^ , k = 4.73 x 10^^ ; tpi = 8.0 x 
10~'^, ipi = — 10~^. These values correspond to a non-slow-roll trajectory with 5 ~ 2, 
Hi = 6.25 X 10~^, 64 = 6.4 x 10"'^. As is also seen from the phase diagram, the bottom-left 
figure, we start far from the slow-roll regime for which 5 ~ <C 1. Despite of this, as we see 
from the top-left figure, after an abrupt oscillation the field i/j loses its momentum and falls 
into the standard slow-roll trajectory. As shown in the bottom-right figure, for this case we 
get a large number of e-folds which seems to be a fairly robust result not depending much 
on the initial value of 6. 
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Plugging the ansatz f l5.43p and the metric (14.60 into T^j/ for the gauge fields, we obtain 
a diagonal homogenous tensor 

T"^ = diag(-p, P-2P,P + P,P + P). (5.46) 

Following the decomposition (15 .23 p . we can write the energy density p and the (isotropic) 
pressure density P as 

P = Pk + Pym and P = -p^ + ^p^,,, (5.47) 
where as before, and p^^ are respectively the contributions of n and Yang-Mills terms 



3 ng^^^ 



= o— (5-48) 



= 2|^3A^(a-'^"+A)a) ^a + ^'-A^a^ + ^A^^ j ' ^'-''^ 

while P, which represents the anisotropic part of the pressure density, is given by 

P = i(l - A') f ^(^ - 2(a + ^jtf _ _ + (2^ - (. + ^(5.50) 

3 \A*a Aa y j A \a \ a I a 

It is readily seen that in the isotropic case A^ = 1 (d = and A = 0), P vanishes and p^j^j- 
and reduce to their expected values ( I5.20p and (15.241) . Note that p^ is only a function 
of (j) and not A which makes parametrization ( I5.44p perfect for studying the quasi-de Sitter 
solutions. 

The independent gravitational field equations are 

(5.51) 
(5.52) 

p 2 

^ - = -^Pym- (5-53) 

From (I5.52P we learn that the evolution of & depends on the anisotropic part of the pressure 
P, i.e. P is the source for the anisotropy a. Thus in the absence of P, the anisotropy & is 
exponentially damped, with time scale a~^. Note that the independent Einstein equations 
are three independent equations for four unknowns, so we need an extra equation which can 
be obtained from the action. 

After substituting the axi-symmetric Bianchi metric and the gauge field ansatz into the 
total Lagrangian ( I5.2ip . the total reduced (effective) Lagrangian is obtained as 

\ ^ a'^ ' 2A^ 

JA6-1) A00 (2 + Ae) A^0^ ,,, (2 + A^) ^V 3/.^V0^\ 
+ ^~)^A^ + ~)^A^^ ^A^ ^ + 2^^^^' ^^-^^^ 
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-^2 = 


P 


3' 




- 3CKT = 




a 4 


-3&' = 


P 



where dot denotes derivative with respect to the time coordinate t. As we see, o" is a cychc 
variable in the above action, therefore, it's conjugate momentum is a constant of motion 
which gives 

^ = : I2T ■ (5.55) 

2a2(3 + A-4(2 + A6)g) 

Here we used this fact that the anisotropy & should vanish for the isotropic gauge field, i.e. 
when = 1. 



5.5.1 Analysis in quasi-de Sitter regime 

Up to this point T^^, the reduced action and & have been determined without imposing the 
slow-roll inflation condition, i.e. the K-term dominance. Hereafter, we simplify and analyze 
the equations assuming quasi-de Sitter inflation, in the sense that the parameter e 

a 

is small during inflation. Instead of d, it is more useful to work with its Hubble-normalized 
quantity T,/H, which is defined as 

§ - -• (5.56) 
H a 

In fact, one should only impose the slow- roll condition on the isotropic sector of the dynamics 
and, as the equations imply, the rest of variables are not enforced to be the slow-roll. 

Combining f l5.5ip and fl5.53p we obtain e in terms of p^, Pyj^j and &: 

2py„ + 



Pk + Pym + 3(^2 



(5.57) 



which demanding a very small e, implies that PyJ^^ and should be much smaller than p^. 
Moreover, noting that P is a term from the contribution of Yang-Mills part, and that the 
energy momentum tensor of the Yang-Mills terms T satisfies dominant energy condition 
(c/. Appendix [D]) , we learn that the absolute value of all the elements of T^^^ are less than 
its energy density p^. Thus, from (15.461) we always have 

P < \pym, (5.58) 
which combining with (I5.52p and (I5.57p . gives 

- < e and e ~ (5.59) 
H p« 

where, as mentioned before, ~ means equality to the leading order of e. 

Using the above relation, (I5.48P and (15. 49 p . and ignoring the S terms in p^^^^ we find 



3a2A4 



9 \2 ^\ 

(2 + A^)(A^ + 4t^) - ^(1 + A^) + (1 + 2A^ 



(5.60) 
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where 7 = {cf. (15.381) ). As we will show through analytical calculations, ^ is at most an 
order one quantity. Therefore, in order to have a successful quasi-de Sitter inflation (e ^ 1), 
we should have 

a " aX 



< 1 and rvi < 1 • (5-61^ 



In other words, recalling f l5.43p . as in the isotropic inflation case, our field values ipi^s should 
have (physically reasonable) sub-Planckian values during the quasi-de Sitter inflation. 

So far we have hence shown that the quasi-de Sitter inflation still means "p^ dominance" . 
Precisely, we have shown that e ~ ^^^^ and / H"^ = |^ < . To proceed further we need 
to analyze dynamical field equations. Using f l5.5ip . f l5.57p and ignoring a terms, we obtain 
the following relations similar to the isotropic case 

^ Pk , 4> — Oicf). (5.62) 

Then, neglecting a terms, from fl5.54p . we can deduce the field equation corresponding to A 

as 

(2 + A^)(AA + aXX + 2^AA) - Q}A ^ + X\X^ - 1)(^ + ^ + ^) = 0, (5.63) 

which, using fl5.62p and keeping the leading orders, can be simplified to 

(2 + A*^)(AA + 3dAA) - QX? + X^{X^ - 1)(2 + X^-i)a^ ~ 0. (5.64) 

One can see that A = is a singular point of the above equation and that the dynamics 
do not mix A > and A < regions. That is, if A is initially positive (negative), it always 
remains positive (negative) during inflation. 

Since fl5.64p is a nonlinear second order differential equation which has no explicit time 
dependence, its solution will be of the form A = A(A) and hence 

X^^A. ,5.65) 

In terms of derivatives with respect to dN = adt, and denoting ^ by a prime, we obtain 

dX' 1 ,6A' X(X^-l), ,,,, , , 

which implies that A'(A) is an odd function of A, A'(— A) = — A'(A). We note that A — )■ —A 
(together with — > —0) is in fact a symmetry of our theory; this symmetry is nothing but 
the charge conjugation symmetry of the original gauge theory we start with. 

Using Mathematica, the above equation can be studied by the phase diagram method, and 
in Figure Uni we have presented the behavior of the solutions in the A' — A plane. Apparently, 
all of trajectories approach to the isotropic fixed point A^ = 1. Next, we give the asymptotic 
analysis to confirm that the isotropic inflation is an attractor in the phase space. 
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Figure 13: The phase diagram in the A' — A plane, the vertical axis is A' and the horizontal is 
A. Both figures show existence of attractor at |A| = 1, corresponding to the isotropic FLRW 
background. The left figure shows the phase diagram over a large range of values for A, while 
the right figure shows the phase diagram for A in the vicinity of the attractor solution A = 1. 
The left figure explicitly exhibits the A'(A) = — A'(— A) symmetry. 

Asymptotic analysis 

In order to investigate the system analytically, it is convenient to rewrite fl5.64p as 

(2 + A^)(^ + 3^) - 6^ + (A^ - 1)(2 + A^) ^ , (5.67) 

Moreover, using f l5.62p in (15. 55 p . we can write S as 

|. + 4, (,«-l))_g,. (5.68) 

From (I5.67p . one can distinguish three different regions for the value of A: A in the vicinity 
of one, A close to zero and large A values. In order to have a better understanding of the 
system, we will solve (I5.67P in these three limits and can then determine S using (I5.68p . 

I) If A is initially in the vicinity of one, we can approximate A as 

A = ±1 + -5A, |(5A|<1, (5.69) 
6 

and equation (I5.67P as 

(5A" + 3(5A' + 2(2 + 7)(5A~0. 

yielding 




(5.70) 
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As we see, the dampinEf ratio is equal to — , ^ , which is less than one and hence SX is a 
^ ^ ^ 2^/2(2+^' 

damped oscillator, damped within one or two e-folds. Therefore, A = ±1 + 6X exponentially 
approaches the attractors at A = ±1. 

Inserting (15. 69^ in (15.680 . we can determine S 

which implies that S and 6X have opposite signs. For A^ in the vicinity of one, a has the 
following behavior 

I — |<e-t"*, (5.72) 

here the subscript denotes an initial value. Thus, |ct| is exponentially damped, with a time 
scale 2/ (3d), as is also indicated by the cosmic no-hair theorems (see Appendix iDj) . 



II) In the limit of very small A values (|A| ^ 1) and considering the leading orders, 
equation (15.671) has the following form 



which can be simplified as 



A" \' 

- + 3j-3^-lc.O, (5.73) 



(i)" + 3(^)' + 2lc.O, (574) 



Solving the above equation, we obtain 

1 



^2 - Aie-^"* + A2e-°*, (5.75) 

which represents an exponential increase in |A| value with time scale of the order d"^. Thus, 
in the limit of initially very small A^ values, |A| is growing very rapidly and escaping quickly 
from the vicinity of zero. As a result, the above approximate solution is only applicable in 
first few e-folds where A^ is far from one. Despite the monotonic increase of A, as seen from, 
interestingly, this is not necessarily the case for a. To see this we evaluate S for two different 
initial conditions in which (i) Ai = and (ii) A2 = 0. 

(i) Putting Ai = in ([5775]), we have A = Aqcs"*, and (ESS]) yields 

S/if ~ -^A^ ~ -^e"', (5.76) 

i.e. |S| is exponentially increasing in time. (Recall that = - is a constant in the 
leading order of e.) However, as mentioned above, this exponential growth of E can 
only be sustained for the first few e-folds, after that A gets close to one and |S| is 
exponentially damped. That is, our gauge-fiation does not strictly follow cosmic no- 
hair theorem [7], as our model does not satisfy the theorem's assumptions. 
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(ii) Putting ^2 = 0, we obtain A = Aoe"*, leading to 



which is quickly damped. 



Note that in this limit ^ 1 and S has always a negative sign. 

Ill) In the limit of large A values (|A| ^ 1) and recalling (15.611) . ^ ^ 1, we have A^7 <^ 1. 
As a result, up to the leading orders, we obtain the following approximation for (I5.67P 

A" + 3A' + 2A ~ 0, (5.78) 

which is identical to (I5.74p that governs the evolution of in the limit of |A| ^ 1. Thus, 
the behavior of A in the limit of |A| ^ 1 is identical to the behavior of ^ in the limit of 
|A| <^ 1. More detailed analysis of this case may be found in |123] . 

To summarize, assuming a system which undergoes quasi-de Sitter inflation in the sense 
that e is very small, we determined A and S. We see that regardless of the initial A values, 
all solutions converge to A^ = 1, within the few first e- folds. Note that, A^ = 1 corresponds 
to two values A = ±1, which are the isotropic solutions. As we saw before, A cannot pass 
through zero during its evolution, so its sign does not change in time. As a result, as we have 
shown analytically and will be demonstrated numerically in the next subsection, system's 
trajectory eventually meets its attractor solution 

A ^ 1 if Ao > 0, 

A ^ -1 if Ao < 0. 

Furthermore, comparing (I5.74p and (I5.78P which are the approximate forms of (I5.67P in the 
limits of very small and very large A^ values, we find out that the behavior of A in the limit 
of |A| 1 is identical to the behavior of in the limit of |A| <^ 1. These results may be 
traced numerically too, which is depicted in FigJTSl 

Although A^ evolves toward the FLRW isotropic solutions, it is shown that in some 
solutions, grows rapidly at the first few e-folds saturating our upper bound of S ~ e for 
a short time. However, this growth stops fast (within a couple of e-folds) and S is damped 
for the rest of the quasi-de Sitter inflation. 



5.5.2 Numerical Analysis 

The degrees of freedom in our system consist of two scalar fields if) and A, the isotropic 
expansion rate a and the anisotropic expansion rate a. Thus, our solutions are specified 
by eight initial values for these parameters and their time derivatives. The gravitational 
equations, however, provide some relations between these parameters. Altogether, each 
inflationary trajectory may be specified by the values of six parameters, ((?, k; t/^q, ipo, Aq, Aq), 
here subscript indicates the initial value. 
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In what follows we present the results of the numerical analysis of the equations of motion 
f l5.5ip . fl5.55p and the and A field equations, for two sets of parameters corresponding to 
two different positive initial A values (Aq = 0.1 and Aq = 10). Similar diagrams for two more 
values of Ao, |Ao — 1| ~ 0.01 has been presented in |123] . The top left figures in Figures 
[H] and [15] show classical trajectories of the field ip with respect to dot, while the top right 
figures indicate dynamics of e. As we see there is a period of quasi-de Sitter inflation, where 
i(j remains almost constant and e is almost constant and very small, which is essentially 
the same behavior we saw in our gauge-fiation isotropic analysis of previous subsections. 
Bottom left and right figures respectively show evolutions of our two dimensionless variables 




Figure 14: The classical trajectory for k = 3.77 x 10^^ g = 10'^, ^jq = 0.6 x 10-^ ^jq = 10"^°, 
Aq = 10, Aq = —3.6. Here A = —a\ which corresponds to = case in (15. 78 p . As expected 
from our analytical calculations there is a short period in which S/if is positive and rapidly 
increasing, but this lasts very short and quickly (in a couple of e-folds) decreases to 
become almost zero. These values lead to a trajectory with do = 4 x 10~^, eo = 0.24 and 
{a/H)o = —5 X 10^^. The initial value of e is rather large, but with in one number of e-folds 
it decreases and reaches 10"^. Note that value of e at the point of maximum T,/H is equal 
to 0.05 (S/i7 ^ ie), almost saturating our upper bound for anisotropy S/if. 
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A and S/if during the first several e- folds. In the left bottom figure of Fig. [TU we see 
that A which started from Aq = 10 quickly decreases and gets close to one. The right 
bottom figure indicates that S/iJ, which is initially equal to 1.2 x 10~^, shows a phase of 
rapid growth and saturates our upper bound S/if ~ e. More precisely, the peak value of 
Ti/H is T,/H\t^^^^ = 0.016 which is about |e|ip^„j.. After its sharp peak, S decreases quickly 
and within few e-folds becomes negligible. At that point, system mimics the behavior of 
isotropic inflation. The left bottom figure of Fig. [15] shows A, initially equal to Aq = 0.1, 
quickly evolves towards one. The right bottom figure indicates T,/H, which is initially equal 
to —2 X 10~^, is exponentially damped and becomes negligible. As a result, after a few 
e-folds, the system undergoes an essentially isotropic quasi-de Sitter inflation. 




Figure 15: The classical trajectory for k = 1.99 x 10^°, g = 10"^ Vo = 0.099, ipo = 10"^°, 
Aq = 0.1, Aq = 1.288 X 10"'^. Here A = dA which corresponds to the A2 = case in fl5.74p . 
As has been predicted by the analytical calculations, S is negative and monotonically, ex- 
ponentially damped. These values give a trajectory with do = 4 x 10~^, eo = 0.16 which is 
a rather large initial e value. As we learn from the top right figure, e decreases within a few 
number of e-folds to 0.5 x 10~^. For these values do = —0.8 x 10~^. 
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5.6 Gauge-flation cosmic perturbation theory 

We showed that the gauge-flation model can lead to a fairly standard slow-roll inflating 
Universe with enough number of e-folds and examined the stability of isotropic FLRW back- 
ground against initial anisotropies. However, as reviewed in section [21 we now have a wealth 
of CMB and other cosmological data which can be used to test inflationary models beyond 
the background dynamics and any analysis of inflationary models is incomplete without 
discussing the cosmic perturbation theory. In this section we work out the gauge-flation 
cosmic perturbation theory. Dealing with non-Abelian gauge flelds brings many new fea- 
tures and complications compared to the standard scalar-driven models reviewed in section 
El Although, as discussed in previous subsections, at the level of the classical background 
the gauge-flation effectively resembles a single-scalar driven model, our consistent reduction 
arguments does not extend to quantum level; in principle, all the gauge fleld modes can 
contribute to the quantum fluctuations around the background. Due to isotropy of the back- 
ground one can still use the scalar, vector and tensor decomposition for the perturbations. 
In what follows we flrst classify all these perturbations and as we will see gauge invariant 
perturbations consist of flve scalars, three vectors and two tensor modes. We study their dy- 
namics, at classical and quantum levels and flnally compute the power spectra for curvature 
perturbations and the tensor modes and the corresponding spectral tilts. 

5.6.1 Classification of gauge invariant perturbations and field equations 

The perturbed metric can be parameterized as in the standard case (I2.24p : 

ds^ = -(1 + 2A)dt^ + 2a{d,B + Vi)dx'dt 

+ ((1 - 2C)5ij + 2dijE + 2d(^iWj) + hij) dx'dx^ , 

where di denotes partial derivative respect to and A, B, C and E are scalar perturbations, 
Vi, Wi parameterize vector perturbations (these are divergence-free three- vectors) and /ijj, 
which is symmetric, traceless and divergence-free, is the tensor mode. The 12 components 
of the gauge fleld fluctuations may be decomposed as 

5A\ = 5%{Y + <PE) + 5i{u^ + <PW,) , (5.79) 

6A'', = 5-(g-|:C) + 5"^'9,fc(Z + 0E) + (70e">9fc(Z-Z) 

£1 

+ 6id,{vj + <j)W,) + <f)eywj + 6''^<j)Chij + ^h,,), (5.80) 

where following our discussions in 15.11 we have identifled the gauge indices with the local 
Lorentz indices and the perturbation is done around the background ansatz (15. 8p . In (I5.79P 
and (I5.80p . we have parameterized the gauge fleld perturbations by four scalars Y, Q, Z, Z, 
three divergence- free vector modes Uj,Vj,Wj and one symmetric, traceless divergence-free 
tensor mode hij, adding up to 4 -|- 3 x 2 -|- 2 = 12. However, as we see explicitly in f l5.79p and 
f l5.80p . due to the gauge (vector) fleld nature of A'^^, it turns out to be more convenient to 
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parameterize the gauge field perturbations in a way which also involves metric perturbations. 
To be more precise, in fl5.80p we recalled (15. 7p and hence 



A\=ije% ^ 5A% = ij5e% + 5^^.A%, (5.81) 

where is the perturbation of the spatial triad and induces the metric perturbations on 
the gauge field, while 5 ^.^^j represents the spatial fluctuations of the gauge field itself. Thus, 
in the 5^4 "j, all the genuine gauge field fluctuations are coming from the and the metric 

fluctuations are from the (5e"j term. Similarly, for the 5A"-q in (15.791) . Y and Uj are chosen in 
such a way that they are invariant under any inflnitesimal space-time transformations, that 
is 

A\ = Q ^ 5A\ = ct^5i{d,E + W,) + 5^^A\, (5.82) 

where 5gfA"-Q is invariant under any inflnitesimal space-time transformations. 

Because of the gauge symmetries of the problem, not all 25 metric and gauge fleld per- 
turbations are physical. The gauge transformations in our case, in addition to the standard 
ones discussed in section [2] which are generated by inflnitesimal space-time transformations 
(I2.27P we also have the inflnitesimal internal gauge field transformations 

5A\ ^ 5A\ - U,X^ - e\^X'A% , (5.83) 

where the gauge parameters A" can be decomposed into a scalar and a divergence-free vector 

A-^ = 5'''d,\ + 51X^ . (5.84) 

In what follows to distinguish the above two gauge freedoms we call the space-time gauge 
transformations (I2.27P "x'^-gauge" and the gauge fleld internal gauge transformations (15.831) 
"A'^-gauge" where 



• x^-gauge act on the perturbed metric 5gfj,v, hence A, B, C, Wi and Vi are transformed 
under the space-time gauge transformations, while they are invariant under the internal 
gauge transformations (I5.83p . 

• A°--gauge act on the the genuine gauge fleld perturbations, S^^A"-^, thus Q, Y, Z, Z, Ui, Vi 
and Wi are transformed under internal gauge transformations, but they are invariant 
under the inflnitesimal space-time gauge transformations. 

• The two tensor modes hij and hij are invariant under both type of the gauge transfor- 
mations, and hence are physical quantities. 



Field equations. After discussing the decomposition of metric and gauge fleld pertur- 
bations we study their dynamics governed by perturbed fleld equations. As in the general 
multi-scalar fleld case reviewed in section |2] these equations are nothing but the perturbed 
Einstein equations and perturbed gauge fleld equations 

5G,, = ST,,, 5(^k^^)=0, (5.85) 
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where by 6 in the above we mean first order in field perturbations. Nonetheless, as in the 
scalar-driven case the integrability condition of perturbed Einstein equations 5(V^G^jjy) = 
implies some of the gauge field equations of motion, but of course not all. 

To see how many gauge field equations are needed to be added to perturbed Einstein 
equations to obtain the full set of equations, let us recall the counting of these equations. 
Among the ten perturbed equations, there are four scalars, two vectors and one tensor {of. 
discussions of section [2]). So, as far as scalar and vector modes are concerned, we need 
one more scalar and one more vector equation. These extra equations may be provided by 
the n = component of perturbed gauge field equations. The equation of motion for A°-q 
component is a constraint enforcing the gauge invariance of the action@ Therefore, this 
extra constraint is independent of the Einstein equations. The full set of equations needed 
to deal with scalar and vector perturbations are hence provided through^ 

6G,. = 6T,^, ^( ^^g^^ )=0. (5.86) 

We stress that to deal with the tensor perturbations, we need to consider another equation 
among perturbed gauge field equations ^(^^g^^) = 0. 

Instead of working with perturbed equations of motion at first order, one may alterna- 
tively consider action expanded up to second order in field perturbations; the equations of 
motion for second order action should obviously reproduce fl5.85p . Working out second order 
action is also needed for the process of quantization of perturbations, to read the canonical 
momentum conjugate to the dynamical fields. Therefore, in what follows we also present the 
part of (but not full) second order action which is require for our analysis. In terms of the 
second order action, the gauge field equation in fl5.86p may be written as 

«4(^^^)^0, Sf{?M^)^o, (5.87) 

where 62 stands for second order in perturbations. The equation of motion for the tensor 
mode hij will also be obtained from the corresponding part of the second order action. 

To work with perturbed Einstein equations we use the same decomposition developed in 
section [2] and decompose perturbed energy- momentum as in (12.251) : 

ST,j =P5gi, + (^5,y(5P - ^V'tt^) + 9,,-7r^ + dii^] + 9,-7rf + vrj 

5r,o =P%o - {dM + ^ql) , 
5Tqq = - p6goo + 6p , 



^^Dcaling with a gauge invariant action, A^q docs not appear in the Lagrangian density L , the momentum 
conjugate to A"q, given as is identically zero. 

^^Note that for a single scalar case the perturbed Einstein equations are enough and they contain the 
equation of motion for the perturbations of inflaton field [2], whereas for the multi-scalar cases with n 
number of inflatons we need to consider n — 1 perturbed equations for inflatons. Similarly, for the models 
with one vector gauge field, again the equation of motion for gauge field perturbation is included in perturbed 
Einstein equations (this may be seen from our discussions of section [3]), while if we have more than one gauge 
field we need to consider gauge field equations. 



100 



where p and P are the background quantities and , nY , vr^ represent the anisotropic 
inertia and characterize departures from the perfect fluid form of the energy-momentum 
tensor; tt^ and 6q\ are divergence-free 3-vectors and a tensor mode, tt^^ is a symmetric, 
divergence-free traceless tensor. 

We summarize the above discussion in the following table. In the left hand side of the 
table we have the fields d.o.f, while in the right hand side we summarized the number of 
independent equations governing the dynamics of each part of the system: 







5,fAl 


-gauge 


^"-gauge 


Gauge- 
invariant 


Einstein 
Eqs 


Gauge 
Field Eqs 


Total 
Eqs 


Scalar 


4 


4 


-2 


-1 


5 


4 


1 


5 


Vector 


2 


3 


-1 


-1 


3 


2 


1 


3 


Tensor 


1 


1 








2 


1 


1 


2 


Total 
d.o.f 


10 


12 


-4 


-3 


15 





Table I: Gauge-fiation perturbation modes 



In the table SgfA^^ represents the genuine gauge field fluctuations, "x^-gauge" denotes 
the space-time gauge transformations and the "A^-gauge" represents the internal gauge field 
transformations. Note also that, in the gauge field Eqs column we have only counted the 
equations which are independent and not included in Einstein or other gauge field equations. 

5.6.2 Treating the equations of motion, classical solutions and quantization 

After determining the field equations, we are now ready to eliminate the gauge degrees of 
freedom and solve the equations. In the following, first by constructing the gauge-invariant 
combinations we remove the gauge d.o.f, then we solve the equations for scalar, vector and 
tensor perturbations. 

• Scalar modes 

In the scalar sector of the perturbations, A, 5, C, E are coming from the perturbed metric 
and, Q, y, Z and Z from the perturbations of the gauge field. Under the action of the 
transformation f l2.27p . the scalar fluctuations of the metric transform as 

A-5t, C ->C + H6t, 

St . (5.88) 

B ^ B ~\ a6x , E ^ E - 6x . 

a 

while the genuine gauge scalars {Q, Y, Z and Z) are invariant under the coordinate trans- 
formations. On the other hand, under the action of the internal gauge field transformation 
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of the form fl5.83p . the gauge field perturbations transform as 

Q^Q, Y^Y--\, 

f (5.89) 

z ^ z , z ^ z - -\. 

9 

We note that Q and Z are invariant under both internal and space-time gauge-transformations. 
Then, we can construct five independent gauge-invariant combinations. One such choice is 
as follows. The standard Bardeen potentials fl2.28p 

m = C + a^HiE ), (5.90) 

a 

d ( 2.-^ B 



from the metric perturbations, in addition to 

Q =Q, (5.92a) 

M =^Z, (5.92b) 

M =^{Z - Y) , (5.92c) 

which are the three gauge invariant combinations coming from the genuine gauge field fluc- 
tuations. 

The scalar part of the first-order perturbations of the gauge field strength is given as 
= St{Q-^C) + 6^^d,,{Z-Y + ^E)+geyd,{{<pZ)-<f>{Z-Y)-<pZ + <f>'E^ 

5F% = 25f.a,](g -j^C- g^\Z - Z)) + 2g<Pe%d,^k{Z + <PE) - 2g<\>e\^{Q - ^C) , 

which is made of gauge-invariant combinations (Q, M and M) as well as Z which is not 
gauge- invariant. However, note that are not gauge-invariant and under an arbitrary 
gauge transformation, they transform as 

F% ^ - eVA^F;,. (5.93) 

It is straightforward to show that in the gauge field strengths ^F"^^ the gauge freedom part 
can be used to remove Z terms and therefore, in the equations govern the dynamics of the 
system Z terms cancel out. 

Since we are working with gauge-invariant combinations, we should write the Einstein 
equations in a gauge-invariant form. In order this we note that (5T^j/ has four gauge-invariant 
scalar parts bpg, 5Pg, 6qg f l2.29p 

Sp = Sp- pa^E - —) , SP, = 6P-Pa^(E-—), (5.94) 

a a 

Sqg = 6q + {p + P)a\E - -) , (5.95) 

a 
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and a^n^ [2]. Perturbing the energy-momentum tensor to the first order, we obtain 

a^TT^ = 2(M-M), (5.96) 

6q, = -2(M + 3H6 X M - HM - ^J^M +-(- - ^^)), (5.97) 

a a an 



— — ~ 77* + 6 1 + —)— 3 1 + — ^ — ' 



+ 3.^$-(1 + ^^)^;a#-2(1 + ^:^)^;m, (5.98) 

6P, = (l-3^)4(g-^vl/) + 2(l-3^)^^-(l-3^)^$ 

- (4,3«T)^*-(i-^£!:^)^!M-2(i-^f)^!M. (5.99) 

a* 3 3 a"^ 

The perturbed Einstein equations then take the form 

a'^dijTi' = dij{^ , (5.100) 
di{5qg + 2{i! + H^))={), (5.101) 

- ?,H5qg + 2—^ = 0, (5. 102) 

(X 

2 k'^ 

dPg + + 3H5qg + 2eif2<|, _ - $) = . (5.103) 

Although it is not independent of the Einstein equations, here for later convenience we also 
write the equation of energy conservation 

/ 2 / 2 / 2 

6p„ - 3H6q„ + 3eH^6q„ - 6H—^ + H—(^ - $) + 2— (^ + H<^) = 0. (5.104) 

As mentioned before, the four scalar modes of the Einstein equations do not suffice to deal 
with five gauge-invariant scalar degrees of freedom and we need one more equation. This 
last equation is provided with the scalar part of f l5.87p . or by the equation of motion for M 
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field coming from tlie second order action for the scalar perturbations 



a a? a?(j)laj CL \ cjp' a? J 

+ I _ u-^h + 2(3 - 2^4)%^ + (12«!|!(1 - 4) + 2^!)^* - 6«!f (1 + 4)U 

a a'^ a'^cp/a a-^ a'^ a-^ a a-^ a 

Kg^'k^ M U^eM\Q ^g^\-. ,g'<\>\.^ ^9^^\,,2 3W^2 



2^— 2— : — — — ^—!—'^^ - 3£--:_^<|) — — ^-i—^ 



a2 a2 a2 g^(f)^/a^ a? a? ^ \ a? ?, a? j a? g^(l)^/a^ 



, g^cf)^ 1 , lKg^(P\k^\k'^ (, Kg^(l)\,(t)- , , g^(j)\ M 

+ ^ + :7>r7T(l + o^)-I -iM' + (1 + ^)(-^ + $) + (2 - e^)vl/ - 2^ 



02 20Va2' 3 a4 ^a2y a2 V > 

M 2ft02A;2 M \A;2 ~ / k02 • ^ k02^ ^ k02^ \ P . 

This equation after using (15.1021) takes the following form0 

ei/^^ + 5x + H^) - -^e^ + -^(^ + $) = 0. (5.105) 



With the above constraint we have enough number of equations for the gauge-invariant scalar 
perturbations. Note that the rest of field equations does not lead to any new equation. In 
particular, as expected, the field equation of $ is equal to (15.1021) . while the field equations 
of M and M reduce to (I5.105p . On the other hand, using the constraint (I5.105p . the field 
equation of Q is equal to the energy conservation equation (I5.104P while field equation of \1/ 
is identical to (I5.103P plus the field equation of Q. 

From the combination of (I5.10ip - (l5.103p and (I5.96p . we obtain 

g204 p 1 Q 1 g203g 12 

^ + i7$ - 2^$ + — 2M - ^ - 2- - - — ^ - 4^- + 2--e^ = 0. 5.106 

a2 a a Ha a-^ a 

Furthermore, one can write (I5.102P as 

(6^. - 3^)2 + il2H^ - 6^ )g - 1^(1 + ^)(v^ + $) + ^(3 + ^)v^ 



k\ Kg^d)^ K(l)\ g^6^ o204^ 2^4 

— ^ 3 + + 2-^ M + 3^$ + 3^— + 3e^^ = 0. 5.107 



Also using (I5.96p . we can omit M in (I5.10ip and obtain 

^ + i/(l + -)$ - M + H{1 + -i)M - (- - '^)Hm - ~ 0. (5.108) 
2 2 a a 



^^The constraint (|5.105p is equal to the gauge field eonstraint equation ( ^^.f ) = 0. 
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Equations flS.lOSp . fl5.106p . fl5.107p and f l5.108p provide enough number of equations for 
Q, M, \1/ and $. Note that the first three equations are exact in slow-roll parameters, while 
the last one has been written in first order in slow- roll. Using ( I5.96P , one can then determine 
M in terms of the rest of variables. 

In order to solve the equations and determine the dynamics, we first write equations in 
the two asymptotic limits of subhorizon if) and superhorizon scales <ti H) and 

then combining them together, we derive the closed form differential equations governing 
the dynamics of our dynamical variable Q. For convenience, we will rewrite the equations 
in conformal time r (r = f ^dt). 



Asymptotic past limit: In the asymptotic past limit kr ^ 1, the equations flS.lOSp and 
f lS.lOSp take the following forms respectively 

A;2(^ + $) = and M'-*' + -g = 0, (5.109) 

a 

where prime represents a derivative respect to conformal time. Using the above constraints 
in (15.1041) and fl5.106p , we can omit \l/ and M in terms of Q and $ which leads to the following 
set of coupled equations for $ and Q respectively 

V + ^^(^^Q + |-$')=o, (5.110) 

a ^7 a 07 

(2^Q + $')" + A;2(2^g + $') =0. (5.111) 
a a 

The above equations, then imply that Q oc in the kr —00 limit. Moreover, fl5.109p 
indicates that we have a non-zero scalar anisotropic inertia a^n^: 

a^Ti^ = -2$. (5.112) 

Note that regardless of the details, for all the scalar inflationary models in the context of 
GR the anisotropic stress a^ir^ is identically zero. 

In the set of equations fl5.110p - fl5.llll) . while the first one has a complicated form, the 

second one is simply a wave equation for 2^(5 + with a sound speed equal to one. Then, 
multiplying the former by a factor of 7 + 1 and subtracting the result from the latter we 
obtain 

(7 - l)^Q" - m" + {^)k'{{i - l)^Q - $0 = 0, (5.113) 
a 37 a 

which is a wave equation for the variable (7 — 1)^Q — with a sound speed square equal 
to (-3^)- In fact, in terms of 

Xi = 2^g + $' and X2 = (7-l)-g-$', (5.114) 
a a 

the set of equations fl5.110p - fl5.llll) is diagonalized into two wave equations for Xi and X2 
with sound speeds = 1 and = (-^) respectively. 

Correspondingly, in the asymptotic past limit: 
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One can decompose Q a.s Q = Qi + Q2 where Qi^2 satisfy 



Qi + k^Qi = , g'2' + ^^k^Q2 = 0. (5.115) 

37 



• Then, we can decompose $ as $ = $1 + $2 such that 

$" + A;2$i = , $2 + ^^A;2$2 = 0. (5.116) 

37 

Besides that, we also have the following two constraints 

$; = (7 - l)H^Q^ , $'2 = -2HijQ2, (5.117) 
which couple $ and Q fields in /cr — — 00 limit. 

We note that, although the sound speed of Q2 is negative for 7 e (0, 2), after fitting gauge- 
flation with the cosmic data, we will see that this region is not physically interesting and the 
data requires 7 > 2 values. 

Up to this point, we worked out the equations governing the dynamics of system in the 
asymptotic past limit and the only quantity which is left to be determined at this limit is 
the the canonical normalized field. In order to read the canonical normalized field, here we 
determine the form of the 2nd order action at the asymptotic past limit which after using 
constraints and in terms of Qi and Q2 can be written as 



SoS,^, ~ / drd X 



(7 - 2) 37 



. (5.11^ 



Thus, our canonically normalized field is given as 

= V2(l + 7)Qi. (5.119) 

The superhorizon limit: In the superhorizon limit /cr <^ 1, fl5.105p takes the following 
form 

On the other hand, up to the leading terms in slow-roll parameters, (15.1071) leads to 

| + ^(e-r/)(<|. + |)^0. (5.121) 

Recalling the slow-roll background relation ( I5.37bl) (5 ~ |(e — r])e), and comparing f l5.120p 
and (15.1211) . we obtain 

^~-i(e-7/)$ and ^ ~ ^. (5.122) 
(p 6 
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Use of the above result in (15.1061) leads to the following equation for $ 

$" - 2(e - r/)J{2$ ~ 0. (5.123) 

Moreover, form (I5.104p . we have the following equation for Q 

g"_ j{2(2 + 8e + 6(e-r/))Q~0, (5.124) 

which indicates that the superhorizon behavior of Q is similar to the Sasaki- Mukhanov 
variable v (12.551) . As we see, in this limit, we have only one equation for both of Qi and 
Q21 similarly both of $1 and $2 are described by the same equation. Moreover, the Bardeen 
potentials $ and \1/ are both constant on super-Hubble scales {kr <^ 1), similar to all the 
other adiabatic perturbations. 

Eq. (I5.122P reveals that at the superhorizon limit, the scalar anisotropic stress a^vr"^ is 
non- vanishing and is given by 

a^n^ ~ (5.125) 

Before this, (I5.112p showed that gauge-flation has a non-zero a?Ti^ at the asymptotic past 
limit. The above relation indicates that this quantity has a non- zero value also at the 
superhorizon, thus is an observable quantity which we will come back to it later. This is 
a unique and specific feature of the non-Abelain gauge field inflation, not shared by any 
scalar-driven inflationary model. 

Working out the field equations of Q and $ in asymptotic past and the superhorizon 
limits, in the following we combine them and read the closed form differential equations 
corresponding to each field and study the system. As we see in (I5.119p . Q is our quantum 
field, while \1/ is the classical field. 

► Quantization of the scalar perturbations 

Combining (I5.115P and (I5.124p . field equations of (5i,2 take the form 

+ )Qi - 0, (5.126) 

z 

Q2 + i-^k' - -)Q2 ^ 0, (5.127) 

with the following algebraic constraint at superhorizon scales (15.1220 

+ Q2 = 0{e)(j)^ at kT<^l. (5.128) 

The effective mass term is given as 

z" 

— - (2 + 8e + 6(e - 7/)):kI 

On the other hand, up to the leading orders in slow- roll, we have !K ^ — (1 -(- e)/r, which 
makes it possible to write ^ as 

z" ly^ -- 3 

— = ^ ^ ^ where z/q ~ - + 2(3e - r/). (5.129) 
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The general solution to the equation (15.126P is a linear combination of Hankel functions 



Qi{k,T) ^ ^^e^(^+2^'5)-/^(giijW(fc|r|) + gii7(J(fc|r|)). (5.130) 

On the other hand, the general solution of fl5.127p is 

^(q2K,J.f^k\T\) + q^I Lf^k\T\)), 7 - 2 < 

'^(^g2<(x/^fc|r|) + g,<)(/Mfc|r|)), 7 - 2 > 0, 



which as we see in case that 7 — 2 < 0, it is expressed as a linear combination of modified 
Bessel functions, otherwise it is expressed in terms of Hankel functions. Note that in fl5.13ip . 
the coefficients are chosen such that in both cases, Q2 has the same superhorizon value. Later 
fitting the gauge-fiation with cosmic data, we will see that 7 > 2. Hence the region 7 G (0, 2) 
in which the square sound speed of Q2 becomes negative, is physically uninteresting and is 
excluded by the cosmic data. 

The undetermined coefficients is fixed by the initial conditions. Note that, despite of 
dealing with two second order differential equations fl5.126p and f l5.127p . which in general 
needs four initial conditions to be fully determined, the algebraic constraint f l5.128p relates 
Qi and Q2 such that we need only two initial conditions, as in the single scalar field models. 
Imposing the usual Minkowski (Bunch-Davis) vacuum state in the asymptotic past limit 
{kr — )■ —00) of the canonical normalized field (I5.119p . we obtain 

Q^or^ = VWr+^)Qi ^ (5-132) 

which corresponds to 

qi = and q^ = q2 = 0. (5.133) 
V2(7 + 1) 

From (15.1281) and after using the asymptotic forms of Hj)^\z) and K^{z) in the z <C 1 limil0, 
one can read q2 as 

:.-YI"-^lV:, , ' fl^-^'V (5.134) 

V 37 y ^2(1 + 7) V 37 y ^ ' 

After obtaining gi 2 and qi^2, the quantum field Q is fully determined and we turn to find 
the classical field $ sourced by Q. 



■^"Thc Hi^\z) and K^{z) functions have the foUowing asymptotic forms in the limit of z ^ 1: 

ffW(z) -l^(.)(|)-^ K^{z) ^ lmC-)-\ 
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► Classical scalar perturbations 



To determine $ = $i + $2? we note that $12 are governed by 

$i + (/c'-^)<l>i ~ 0, (5.135) 

^'2 + {^^k^-j)^2 ^ 0, (5.136) 

dynamical equations the solutions of which are subject to the constraint equation 

$; ~ (7 - l)ijHQi, at $'2 ~ -2i>HQ2 , (5.137) 
in the asymptotic past limit. Here ^ = 2[K^(e — 77), which can be written as 

^ = ^^^^, where i^r ^ + 2{e - rj). (5.138) 
Similar to Qi,2, the general solution of the above equations are given as 



<^>i(^,r)^^(6i<)(A;|r|) + 6ii/(5(A:|r|)), (5.139) 



^(z62i/S(V^fc|r|) + h,E^^l{J^k\r\))^ 7 - 2 > 0. 



and 

$2(fc,r)~<^ X^'^^ --V^' '-"^V ^ ' (5.140) 

l^(^^.^S(v^^M) + 62i/£Hv^ ■ 

In ( I5.140p . the coefficients are chosen in such a way that for both cases $2 satisfies ( I5.137P 
with the same value. Putting f l5.133p and fl5.134p into fl5.137p and after using the asymptotic 
form of Bessel function^ in the /cr — t- — 00, we have 61 = 62 = 0, while 

— \\ 2i / h — 2\\^ 
hi^ ^ ^ Hih and ^2 ^ , — -] Hip. (5.141) 

Havin g th e coefficients above and using the background slow-roll relation e ^ (1 + 7)V^^, we 
obtain [f^ 

$ = $i + $2^^i/(^)^~'", A:|r|«l. (5.142) 



•^^Thc Hi^\z) and Ki,{z) functions have the foUowing asymptotic forms in the limit of z ^ 1: 

\ ZTT \ ZTT 

^^In the asymptotic past hmit $2 is a mode which can have negative for 7 < 2. Nonetheless, our 
analysis above shows explicitly that this does not render our perturbation theory analysis unstable, because 
what is physical is the total $ after imposing the constraint equations on the superhorizon scales. In other 
words, $2 niode in the asymptotic past is fixed by the constraints on the dynamical equations and not an 
independent mode. Note also that what matters in the CMB data, is the superhorizon field values. 
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We are now ready to compute the power spectrum of $ and curvature perturbations. 
The power spectrum for the Bardeen potential $ is given by 

A|=(^l*P. (5.143) 
which after using f l5.142p and f l5.138p . has the following form on the large scales {k ^ aH) 

Recalling f l2.35p and using f l5.122p . we learn that the comoving curvature perturbation is 
given by 3i ^ — in the superhorizon scales. The power spectrum of the comoving curvature 
perturbation !Ji is then given as 



A? ^ 



1 /H'' 



TT 



(5.145) 

k=aH 



which becomes constant on super-Hubble scales. Note that the scalar power spectrum in our 
model is exactly equal to the power spectrum of the comoving curvature perturbation in the 
standard single scalar field model fl2.6ip . The spectral index of the curvature perturbations, 
77.;; — 1 = 3 — 21^11, to the leading order in the slow-roll parameters is 

Us-l^ -2(e - r]) . (5.146) 

We note that the spectral tilt fl5.146p is always negative in our model. 

In addition to the power spectrum of the scalar and its spectral tilt, as one of the specific 
feature of the non-Abelian gauge field inflation, we have a non-zero scalar anisotropic stress 
value with the following power spectrum 



(5.147) 

k=aH 



which becomes constant on super-Hubble scales and is hence a physical observable. This 
is in contrast with all kinds of scalar inflationary models in the general relativity, in which 
a^ir^ is identically zero. 

A non-zero a^s is a feature present in the infiationary/cosmological models with modified 
gravity [136] and causes difference between the two Bardeen potentials $ and "if, this provides 
a way to trace its observable effects by affecting the structure formation analysis, for a more 
detailed discussion on this point see |136j and papers referring to it. 



• Vector modes 

In the vector sector, we have Vi, Wi, Ui, Vi and Wi which under the action of an infinitesimal 
"vector" coordinate transformation (12.360 

= x' + 5x^ , {di5x[ = 0) , 
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transform as 



Vi^Vi- a5x\, , Wi^Wi- 5x\, , (5.148) 

while Ui and Vi as our genuine gauge field fluctuations remain invariant. On the other hand, 
under the vector part of infinitesimal gauge transformation fl5.83p . 

Ui^Ui--Xy, Vi^Vi--Xy, Wi^Wi + Xy, (5.149) 

9 9 

while Vi and Wi obviously remain invariant. 

We can construct three gauge invariant divergence-free vector perturbations, one from 
the metric fluctuation 

Z^ = aWi-Vi, (5.150) 
and two from our genuine gauge field perturbations 

1 1 
Ui = -Wi + Ui , and Vi = -Wi + Vi. (5.151) 

9 9 

The vector sector of the first-order gauge field strength perturbations are 

6Fl^ =6id,i% - U, + <PW,) + eV(^70U, + g<l)'W, + <j>w,) , 
6F% =2e\ ''d,] {g<pVk + 9<P'Wk) + 2g<f)'6^ ^,w,] . 

As mentioned before, are not gauge-invariant, but transform as F^j, — )■ F^^, — e^^^A^F^^. 
Thus, Wi terms in the field strength tensor above are pure gauge terms such that all the 
physical quantities are independent of Wi. 

In order to investigate the dynamics of this sector we first work out vector parts of the 
perturbed energy-momentum tensor, SqY and nY ■ At the linear order in perturbation theory 
only the vector perturbations contribute to the vector part of energy-momentum tensor 
perturbations, thus we have 

S,v ^ _2£V(k, + *Z,V4(vx(^-U)) -Sj^fv.v) , (5.153) 

arrY = + ^(U, - V,) . (5.154) 

As mentioned before, both of them are gauge-invariant, hence can be written in terms of 
gauge invariant variables. 

The perturbed Einstein equations involves two vector equations, one constraint and one 
dynamical equation, given as fl2.38p 

di (^2a^7rJ - ^(a^^j)^ = , (5.155) 
2adqY + V^Zi = 0. (5.156) 
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On the other hand, deahng with three unknowns, the two Einstein equations are not enough 
to fully determine the system and we need one more equation, which as discussed in section 
15. 6. H is provided by the vector part of the second order action is the field equation of SA'^q. 
Explicitly, once we write down the second-order action for the gauge field perturbations, the 
momentum conjugate to Ui is vanishing (the vector part of f l5.87l) ). yielding 

_ 2i^(U. + ^Z.) H- 4(V X (f + :*2)). - ^(V X V). - ±V\U. - V.) = 0. (5.157) 
Using f l5.156p . the above equations leads to the following equation 

^ (V X (U + ^Z))^ - 4 V'(U, - V,) + ^V'Z, = . (5.158) 

This completes the set of equations we need for solving vector perturbations. Then, the 
combination of fl5.155p -f l5.156p and fl5.158p indicates that Z is damping exponentially during 
the inflation. From the combination of fl5.153p and fl5.155p . we then learn that Zi vanishes 
after horizon crossing. 

The above is the usual result of the scalar-driven inflationary models that the vector 
modes are diluted away by the (exponential) accelerated expansion of the Universe during 
inflation. In our model, despite of having vector gauge fields as inflaton, the power spectrum 
of the vector modes are unimportant in inflationary cosmology. 



• Tensor modes 

The tensor perturbations hij and h^j, as symmetric, traceless and divergence-free tensors, 
are both gauge invariant with two degrees of freedom. The contribution of these modes to 
the linear order perturbed gauge field strength is 



(5.159) 



The perturbed Einstein equations involve one equation for hij, sourced by the contribu- 
tion of hij to the energy-momentum tensor, reads as (12.401) 

+ 3i7/iy- + —hi, = 27Tf. . (5.160) 
Computing the linear order energy-momentum tensor, we obtain nfj 

= (2(^ - - ^(2^.- + h.,) + ^5,(2e-,/.,, + e-,/^,,)) . (5.161) 



112 



This equation and the equation for the tensor perturbation of the gauge field hih is provided 
by the second order action of the tensor modes which after using the slow-roll approximation 
{(f) ~ Hcf)), is given aJ^ 

54^) ^ \ [ d'xdta' (hi + ^)ihl - i\d,K,y) + ( - - ^{2h, + k,))h., - ^-^hl 
+ —{hij + hij)hij - 2(2— H —)hij jdk{hij + hijjdkhij + -(( —) + 2 — )— e ^ hkidihji 

Qj Qj Qj Qi ^ Q Q Q 

a"* a a'^ 2 a-^ J 

Having two degrees of freedom, each of hij and hij can be expressed in terms of the 
standard plus and cross polarizations which are the eigenvalues of (V^e"*^!^-^ = —k'^e^^^). 
However, they are not eigenvalues of parity- violating operator ( e^^^dk) and in order to 
decouple two degrees of freedoms in each of hij and hij, we need to write them in terms of 
the right and left circular polarizations. We can use the following parametrization for hij 



~ 2a 



( h^ + h^ _,(/i^_/,J o\ 

<h^-hj -{h^ + hj 
\ 0/ 



(5.163) 



where working with Fourier modes, we chose = (0, 0, k) and imposed the transversality 
condition. In a similar way, one can parameterize hij and Tifj in terms of right and left 
circular polarizations h^ ^ and tt^ ^ . 

In terms of h^ ^ and h^^ ^ , (15.1601) reads 

+ [k' - (2 - e):K2)/.,,, ^ 2a27rJ^, (5.164) 

where 

^ ^'(-2^^k. + 27^K^^,,, T 2k^VliKL + lh,J - 'Kh'^^ + 'K'h,,). (5.165) 



Here r is the conformal time dt = adr, ip = (p/a is the background effective inflaton field, 

2 / 2 

"K = d, prime denotes derivative with respect to the conformal time, and 7 = ^jp-- Equation 
(15.1641) implies that in the superhorizon limit, we have /i^^ oc a. 



On the other hand, the field equations of /z^ ^ are obtained from the second order action 
( 15.162P which after using ( ]5.33cp leads to 

+ (k' + (2(1 + 7) + e) J{2 ^ 2k^^^^^^ V, ^ nK.L - ^KJ ± k:K^-^h^^,. (5.166) 
A similar result and feature is also obtained in the chromo- natural model [137) . 
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Figure 16: undergoes a tachyonic growth phase in f2{'j) > —kr > 1 [cf. f l5.167p and 
fl5.168p ). In this figure, we have depicted f2 vs. 7. The minimum is f = 5 which is at 
7 ~ 0.6. 



Note that using the superhorizon scale behavior of /i^^ (^bz, a), from f l5.166p . we learn 
that while /i^^ behaves like a plane- wave at subhorizon scales, it is exponentially damped 
like hj^^ (X a"'-^"'"'^-' at superhorizon scales. However, the parity violating term plays an 
important role just before the horizon crossing, leading to tachyonic growth of around 
the horizon crossing. 

Before the superhorizon scales, we can neglect h^ ^ terms in the RHS of fl5.166p . which 
leads to the foUwoing wave equation for Kr^l 

dfhR,L + culdr, l)hR,L ^ 0, (5.167) 
where f = —kr and ujRj;^{f,'j) is given as 



2 

R,L 



While u1 is always positive, there is an interval f G (Ti,'f2) in which uj^ becomes negative. 
This short interval which has a negative uj^, leads to the tachyonic growth of hji (Fig. [T7|) . 
In Fig. [161 "we presented f2 vs. 7 and fi is almost one. 

The source term vr^^ (the RHS of (I5.164P ) vanishes at superhorizon scales kr 1. 

Nonetheless, due to the tachyonic growth of hn at the vicinity of the horizon-crossing, vr^ 
has the behavior of an impulse function in that region (see Fig. ITTl) . inducing the growth 
in and enhancing its superhorizon value. On the other hand, ttJ is small at the horizon 
crossing and has negligible effect on the superhorizon value of /i^. 

In the leading slow-roll approximation the standard Minkowski (Bunch-Davis) vacuum 
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Figure 17: This figure presents the tensor modes solution for = 5 x 10 ^, 7 = 10 and 

Hq = 10~®. In the top-left panel, we have the tensor field values ^ and ^ versus —kr, 
where Re and Im denote read and imaginary parts of the corresponding quantity. The small 

box presented the superhorizon behavior of the fields. The top-right panel shows In 
the bottom panels we presented the left-handed polarizations. 
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Figure 18: In the left panel we have depicted the sum of Left and Right gravity wave power 
spectra Pj^ + P^. In the standard scalar-driven inflationary models = = 1. The right 
panel the parity violating factor p^_^p^ versus 7 for ip = 10~^ and = 0.12 is shown. The 
power spectra have been calculated at kr = —0.01, long enough after modes have crossed 
the horizon and behave quite classically. As we see in the right panel, for very small and 
very large 7 values P^ ^ P^. 



normalization for the canonically normalized fields (c/. fl5.162p ) takes the form 



-ikr 



-ikr 



and hj^^ — )■ 



kr — )■ —00 . 



(5.169) 



The power spectra for the Left and Right moving gravitational wave modes is then obtained 
to be 



A. 



P 



k=aH 



k=aH 



(5.170) 

(5.171) 



where P^, P^ are functions of the parameters of the slow-roll background, in particular 7, ip. 
Numerical analysis reveals that P^ is a function very close to one (ranging from 1.0 at low 
7 to 1.25 at 7 = 10) while P^ varies quite considerably in its range. Since their explicit 
analytic expression is not illuminating we have only presented their ratio as function of 7 
has been depicted in Fig. [181 In our model, we hence expect to see birefringent gravity 
waves. 



Moreover, the spectral index of tensor perturbations, ut is given by 

ut — — 2e , 



(5.172) 
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which as expected are equal to their corresponding quantities in the standard scalar infla- 
tionary models, given in f l2.65p and fl2.66p . Note that due its exponential suppression on the 
superhorizon scales, the h mode does not contribute to the tensor power spectrum. For our 
model tensor-to-scalar ratio r is 

T = %{P^ + P^)e, (5.173) 

which can be written as r = — 4(Pjj + P^)nT- That is, our model respects a modifled version 
of the Lyth consistency relation [31]. This feature is also shared by the chromo-natural 
model [137] to be reviewed in the next section. 

5.6.3 Summary of gauge-fiation cosmic perturbation theory 

Let us summarize the main observational predictions of gauge-flation. Although it can be 
described by an effective scalar fleld at the level of the background, at the perturbation level 
the gauge-flation model deals with many flelds and has speciflc and interesting features: 

• Scalar modes: we have 5 gauge-invariant scalar modes, two of which are the standard 
Bardeen potentials and three of them are coming from the gauge fleld perturbations. As 
has been explicitly shown here and in |122j . due to the different constraint structure 
that we are dealing with, our scalar perturbations are adiabatic. In other words, 
similar to the scalar-driven inflationary models, all of the 5 physical scalar flelds gets 
constant at the superhorizon scales. On the other hand, being a non-Abelian gauge fleld 
inflationary model, gauge-flation predicts a non-zero scalar anisotropic stress o?tt^ ^ 0, 
in contrast with the scalar-driven models. 

• Vector modes: there are 3 gauge-invariant vector modes which are exponentially damp- 
ing, similar to the scalar-driven inflationary models. 

• Tensor modes: we have two tensor modes, one of which is the usual metric tensorial 
part, while the other is coming from the contribution of the gauge fleld. The former 
gets constant at the superhorizon scales, however, the later is diluted away at the 
superhorizon limit, thus unimportant in inflationary cosmology. Due to the existence 
of parity violating terms in the tensor sector,the right-handed and left-handed circular 
polarizations are not equal, but the right-handed mode is enhanced by the large value 
of its corresponding tensor anisotropic inertia tt'^ at the horizon crossing. 

• It is instructive to state the above in a different wording: In our gauge-flation action we 
have metric and SU (2) gauge flelds, which in the standard flat background terminology 
these yield one divergence-free traceless two-tensor (the graviton) and three transverse 
gauge flelds W^'s, since each gauge fleld has two polarizations, altogether we have 
eight propagating degrees of freedom. On the FLRW inflationary background and in 
our decomposition, where we have identified the gauge (internal) and spacial indices, 
the graviton appears as before, while the three gauge flelds will further decompose into 
a propagating transverse vector, two propagating scalars and one propagating tensor. 
Explicitly, let us focus on VF^. The temporal components may just be set to zero, and 
hence we deal with VT", such that VWf = 0. The propagating vector is the part of 
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Wj", Wi, where Wf" = 5^°'djWi. The remaining W"" is divergence-free on both of i and 
a indices and may hence be decomposed as 

where Wi,W2 are the two propagating scalars and h'} is the divergence- free traceless 
tensor. One can directly identify these propagating degrees of freedom in terms of our 
earher decompositions: The two scalars Wi,W2 correspond to the two modes Qi,Q2 
introduced in fl5.114p . The propagating vector is a linear combination of Uj and Vj 
which is governed by f l5.155p - fl5.156p and f l5.158p . (Note that these equations only 
involve first time derivatives of Hi and Vj, nonetheless they are mixing the two and 
once we eliminate one of the variables we end up with a single mode with second-order 
time derivative equation.) Finally, the tensor mode is exactly the mode governed by 

( Kim . 

In the following table we have recollected all the results of this subsection and compared 
them with their corresponding values in the standard single scalar model which has been 
reviewed in section O 





Gauge-fation 
model 


Single-scalar 
model 


A? 






Us -1 


-2(e-r/) 


-4e + 2r] 


r 




16e 


riT 


-2e 


-2e 




87r2 


identically zero 



Table II: Gauge-fiation summary results 

As we see, these two modes predict the same scalar power spectra, but different scalar spectral 
tilts and tensor power spectra footnote Note that is written in terms of the Hubble slow- 
roll parameter 77 fl2.1ip .. Moreover, as unique and specific features of the non-Abelian gauge 
field inflation, not shared by any scalar-driven inflationary model, gauge-flation predicts non- 
zero scalar and tensor anisotropic inertias, a'^ir^ and a'^ir^ . Regardless of the details, the 
anisotropic inertia a^n^ is identically zero in all scalar inflationary models in the context of 
general relativity [2] and a non-zero anisotropic inertia power spectrum in our mode 

is a consequence of having gauge fields in the system. As we see for our model 




(5.174) 
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As we see, the anisotropic inertia power spectrum, unlike what is usually perceived, is not 
necessarily attributed to modified gravity and in particular Galileon models |138] and these 
models can have different predictions for A^s/A^, than the gauge-flation. It is desirable to 
study the anisotropic stress and its cosmological observable consequences more thoroughly. 

The existence of a non-zero tensor anisotropic inertia tt^^ with different values for vr^ 
and vrj at the horizon crossing leads to the interesting effect of parity violating gravity wave 
power spectra, i.e. the right-handed polarization of the tensor modes gets enhanced while 
the left-handed mode is almost equal to the tensor modes in the standard scalar driven 
infiationary models. The observational consequences of such parity violation on CMB has 
been discussed in [TTJ [H] . 



5.7 Confronting with CMB observations 

Having worked out the basic observables of the gauge-flation model, i.e. power spectrum of 
curvature perturbations, spectral tilt, power spectrum of gravity waves and its tilt, we are 
now ready to confront the model with the data and restrict the parameter space of the model 
using the WMAP7 results [20], reviewed and summarized in section O First, we note that 
in order for inflation to solve the flatness and horizon problems it should have lasted for a 
minimum number of e-folds iVg. This amount of course depends on the scale of inflation and 
somewhat to the details of physics after inflation ends [2] . However, for a large inflationary 
scale, like H ~ 10~^ — 10~^Mpi, it is usually demanded that ~ 60. As for the CMB data, 
one may use the best-flt values given in fl2.69p - fl2.7ip . As is seen from Fig. [T], however the 
range for in the best-flt values, is associated with r = case. However, to perform a more 
precise analysis, we consider the allowed parameter space in — r plane depicted in Fig. [H 
together with the COBE normalization for A^. 

To perform the analysis, we use the results obtained in the slow-roll regime. We then note 
that the gauge-flation data are summarized in Table II are functions of e, t] and H. Number 
of e-folds A'^e f l5.4ip may also be written in terms of these three recalling that 7 ~ (e — t])/?]. 
Note also that H appears only through the COBE normalization and its value may be flxed 
once we obtain the bound on e and 77 using the other data. So for the moment we focus on 
rig, Nf,, r and choose to work with e and 7 as variables to perform analysis. Recalling our 
previous results we have 

Us 



Since < f{-f) < 1, > 1 
example, as depicted in [T9| 

ni^°)> 0.98, (5.177) 

and of course we have a less red spectrum (i.e. lowest value for Ug) if we demand larger A^e- 



= 1-^, /(7) = 7ln(^)<l (5.175) 

A^e 7 

— That is, in gauge-flation spectrum cannot be very red. For 
for' A", > 50, 
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Q\ , , , i , , , i , , , i , , ^1 , , , ^ 

0.92 0.94 0.96 0.98 1.00 1.02 



Figure 19: The shaded region exhibits the region which leads to A^e > 50. As we see this 
may happen for any value of 7 parameter. This also shows that our spectral tilt is always 
in 0.98 < ris < 1 range. Moreover, our model allows for arbitrary large Ng. 



Similarly, as depicted in the left panel of Fig. [201 we see that our model predicts a 
minimum value for r 

0.02<r< 0.28. (5.178) 

This is a very specific prediction of our model and gauge-flation may be falsified by the 
upcoming Planck satellite results. In the allowed region 

=^ (0.01 -0.1)Mpi, (5.179) 

where the max and min possible values of r respectively correspond to [ip = 0.01,7 = 
5, + = 6.3) and {ip = 0.01, 7 = 8, + = 77). Moreover , in the right panel of Fig. 
[2Ul we present the allowed region in terms of e and 7 which indicates that 

£ = (10-^-2x10-2), 7 = (0.1 -8). (5.180) 

Using the above and the COBE normalization, we can read the value of H: 



H \ r tt"^ 



-A? = 2 X 10"^e, 



and hence 

= (0.45 - 6.3) X 10-^ Mpi . (5.181) 

To determine the two parameters of our model, n and g, we can read them in terms of 
the parameter 7 

2/2 2 

7 = ^ |^=27rA27(7 + l) = 1-5x10-^(7 + 1), (5.182a) 

«:^7V^2 ^c^^^i±2)l = io7 y,t±ll, (5.182b) 
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0.92 0.94 0.96 0.98 1.00 1.02 



Figure 20: The left panel shows la and 2a contour bounds of 7-year WMAP+BAO+HO. 
The yellow area (region with lighter color) represents the gauge-flation predictions for il) G 
(0.01, 0.12) range. As depicted in Fig. [191 the region with enough number of e-folds restricts 
us to Us > 0.98 region, that is on the right-side of the = 50 line. Therefore, the allowed 
region is the highlighted region between = 50 and rig = 1 lines. The shaded region in 
right panel shows the allowed values for e and ijj, given in fl5.179p and fl5.180p . 

As we see from the left panel of figures [20} the value of 7 is restricted as fl5.180p . which gives 

^ ~ (0.15-3.7) X 10~^ A ~ (10"^ - 10-^)Mpi, k = A'K (5.183) 

As we see from (15. 182b ). smallness of the gauge coupling is directly related to the smallness 
of the power spectrum of CMB curvature fluctuations. This is a notable and interesting 
feature of our model. 

Nonetheless, if we restrict ourselves to la contour in the left panel of Fig. [201 we find 
more stringent bounds on r, Ug and hence H: 

0.98 < < 0.99, 0.05 <r< 0.15, ~ (3.4 - 5.4) x lO'^Mpi, (5.184) 

which leads to the following bounds for 7 and ip 

0.04 <^ <0.1 and 0.6 x 10"^ < e < 1.5 x 10"^ 

As we see the field ip is sub-Planckian and of order 10^^ GeV. Moreover, recalling the 
connection to massive axion model (c/. discussions of section [5. 2. ip and in particular f l5.3ip . 
one can see that ~ 5 x lO^^GeV and A ~ 2.5 x 10^. We will discuss in the next section, 
these values are all in a very natural range for energy scales and coupling of axions in beyond 
standard models. 
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6 Chromo-natural inflation 



One of the main theoretical issues concerning most of inflationary models discussed in the lit- 
erature is the "naturalness" issue. For example, in the simplest m^0^ or chaotic models, 
the field values are super-Planckian, ~ 10 — 20 Mpi, while the m ~ 10~^Mpi or A ~ 10~^^. 
(In these models the Hubble during inflation H is saturating the current observation bound, 
H ~ 10~^Mpi.) The naturalness issue, besides its theoretical unappealing feature, brings 
puzzles and problems in embedding these models within other high energy physics models, 
like GUT's or string theory. One of the issues regarding the naturalness of inflationary po- 
tentials V is coming from the slow-roll conditions, that V ~ (10"'^ — 10~^ ^pi)^, while its 
derivatives in Hubble units should be small, of percent level or lower. Without embedding 
in a UV completed theory one cannot guarantee protection of the inflaton potential against 
quantum loop corrections. Therefore, there has been many efforts in building inflationary 
models which are natural, speciflcally within a given consistent high energy theory setting, 
e.g. see [HE]. 

One of the main theoretical obstacles in building "natural" inflationary models is to 
flnd a setting which avoids fine-tuning in the potential (coming by demanding the slow- 
roll conditions) by relating it to an approximate symmetry in the potential. The simplest 
such approximate symmetry which does the job is the translation in the field space. Moti- 
vated by this idea the "natural inflation" model was proposed [139] . This model proposes 
using non-Abelian gauge theory axions as inflaton fleld. The axion field x classically cou- 
ples to -F A F term of non-Abelian gauge theory and is hence invariant under x ^ X + Xo 
for arbitrary Xo shifts. This shift symmetry is, however, broken to a discrete subgroup 
considering quantum effects (instanton contributions), which induce a cosine-type poten- 
tial V{x) ~ /^"^(l + cosx//), where and / are two scales [130] . The axion potential is 
pertubatively exact. 

The natural inflation model of |139j . although ameliorates the naturalness problem by 
seeking shift symmetry and that the axion potential does not receive perturbative quantum 
corrections, does not fully resolve it. As successful slow- roll inflation model is obtained for 
super-Planckian / parameter (while //^ ~ H'^M^^) which is not a natural scale within particle 
physics models. It is, however, noted that considering non-Abelian gauge flelds and their 
coupling to axions can resolve this problem |140j . In this section we review the chromo- 
natural inflation model of |140l I132j and study its connection with the gauge-flation model 
of previous section. 

6.1 The basic setup 

Axions, although gauge singlets, couple to the topological F A F term of the non-Abelian 
gauge theories and have a nonperturbatively induced potential. This system describes the 
chromo-natural inflation model with the Lagrangian 



c.n. 




Fa'" + l{d,xf + Al + cosh 



X 

8f 



A V- ^ ap-^ III 



(6.1) 
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where we have set SnG = M~^'^ = 1, ^ and / are parameters of dimension of energy, A is an 
order one dimensionless couphng. In our conventions the axion field x take values in [0, irf] 
range. 

As discussed in |140] . slow-roll inflationary dynamics for the above action can happen 
for yU, / -C Mpi, only if we also turn on a "rotationally invariant gauge field" of the form 
[ml [177] 

' a{t)ilj{t)Sf , fi = t 



in the temporal gauge. This is of course of the same form as discussed in the previous 
section in the gauge-fiation setup. Our discussions of section 5.1 regarding the preservation 
of isotropy and homogeneity with the above gauge field configuration almost trivially extends 
to the chromo- natural model too. In a similar manner one can show that reduction to the 
"isotropic and homogeneous sector" described by ip and x fields is a consistent reduction 
and the reduced Lagrangian may be simply obtained by plugging the gauge field ansatz (16.21) 
into the action (16. ip . to obtain 



'^reduced '^gr ~l~ 



^(^ + - ^/^^ + - /i^(l + cos j) - ^xi^M: + Hi:) 



(6.3) 

where dot denotes the derivative w.r.t. the comoving time t. One may then work out the 
equations of motion for the above system: 

^ + 3Hjc~jsmj = -3^^2(^ + Hij) , (6.4a) 

+ SHi; + {H^ + H)ij + 2/V'' = y V^':^ , (6.4b) 



2 

H = - 



(^ + Hi;)' + g^r] + + ^(1 + cos §) , (6.4c) 
J 6 3 / 

(^ + Hij)'' + /V'] - ■ (6.4d) 



We note that the axion-gauge field coupling term, the A-term, is metric independent and 
hence does not contribute to the energy momentum tensor of the the system. The effects of 
this term appears in the dynamics of the system through the terms in the right hand side of 
(16.4k ). (16.4b ). The A-term plays the important role of taming and flattening the steepness 
of the axion potential, without requiring the parameter / to be super-Planckian. In other 
words, considering the axion-gauge field coupling and turning on the gauge fields makes the 
model technically natural, justifying the name "chromo-natural" inflation model. This setup 
is particularly interesting because the axion potential and coupling terms are perturbatively 
and non-perturbatively exact. 

We will show in this section that equations (16. 4p indeed allow for standard slow-roll tra- 
jectories for which e = —H/H' and x/{HMpi), ip/Hip remain small during the inflationary 
period. Moreover, interestingly this happens for a large region of parameter space and ini- 
tial values which are technically natural. In our analysis we will first assume possibility of 
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slow-roll dynamics and then check the validity of assumption, both analytically and numer- 
ically. We can now look for slowly rolling inflationary solutions of this system of equations 
assuming x, ip, and if ~ and studying the resulting equations. In this limit (16. 4k ) reduces 
to a simple equation 

sin j ^ -^.e-'^'F^pF^^, ^ ^^^^^ ' ^^"^^ 

where, as before, ~ stands for equality up to the first order in slow-roll parameters e, rj, and 
for axion field values which satisfy (16. 5p we will drop the tilde and use and to obtain the 
second equality we have dropped against Htp. 

Alternatively, one may recall the potential for the axion field x'- 



?7(x,H) = /i^(^l + cos^ + S^J , (6.6) 

where S = —-^^°'^^'^F^pF^y This potential has an extremum at x = which to first 
order in slow-roll is given by (16. 5p . The potential has extrema if S < 1. This latter puts an 
upper bound on ijj (for given g and ii)o Eq. (l6.5p has two solutions, xl f ^-nd t\ — xU ^^-nd 
our slow-roll trajectories are such that the axion field rolls down along the minimum of the 
potential and inflation ends when x is close to vr/. It is easily seen that the solution with 
X > 7r//2 corresponds to the minimum of potential (16.60 and the one with x < 7r//2 to 
maximum. This suggests two "natural" initial conditions for axion field xo to start slow-roll: 
a "small axion model" where we start close to the top of the potential, and "large axion 
model" one with Xo/f close to vr. The former is the model discussed in [140] while the latter 
is nothing but the gauge-flation model discussed in the previous section (c/. discussions of 
section [5.2. ip . However, as has been shown in [142j and we will review in section [631 it 
is possible to get a successful slow-roll inflation for any arbitrary value of Xo- In addition, 
although both tp and x fields are slowly rolling, their rolling is "adiabatic" in the sense that 
during slow-roll inflation x and ip fields vary such that the minimizing condition (16. 5p holds. 

As may be seen from (16.4b ) or from (16. 5p during the slow-roll the energy is generically 
taken from the axion sector and injected into the Yang-Mills gauge field sector. In this 
process, which is the new ingredient of the chromo-natural model compared to the natural 
inflation model of [139j . the speed of axion rolling is reduced through the coupling to a speciflc 
combination of gauge fleld condensate. This condensate is proportional to the instanton 
number density F A F (though in the Minkowski signature) and it the cosmological setups 
(for K = 1 FLRW cosmologies) this has been dubbed as "cosmological sphaleron" [143[ 
I127j . Nonvanishing cosmological sphaleron conflgurations in Abelian and non-Abelian gauge 
theories can be a basis for lepto/baryo genesis model building. The former has been discussed 
in [144] and the latter in [145] . 



^^We note that in the discussions of |140j (|6.5p has been viewed as the condition minimizing effective 
potential for the tp field. 
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6.2 Chromo-natural vs gauge-flation 



If tlie value of x at minimum, x, is very close to vr/ one may expand the axion potential (16. 6p 
and to second order in x we obtain the potential for a massive axion discussed in (I5.30p and 
hence the chormo-natural model will reduce to the gauge-fiation, if the initial value of the 
X field Xo is close to tt/. That is, for the "large axion" region, the chromo-natural model is 
equivalent to gauge-flation. This can be made in more robust technical setting by integrating 
out axion field. To this end, one may expand the axion theory around the minimum of its 
potential at x = X (for X > ^/2) given in (16. 5 p and, integrate out the fluctuations [146] . 
This may be done by writing x = X + ^X^ expanding the action up to the second order in 
6x and computing the one loop effective potential [131] f§ 

Ueffix) = U{x) + ^ In det(n + ^ cos , (6.7) 

where U{x) is the potential (16. 6 p computed at the minimum x = X and □ is the Laplacian 
computed on the (almost-de Sitter) inflationary background. 

To compute the Indet term explicitly, as in standard fleld theory treatment, one should 

4 

integrate over the spectrum of the cos j operator. However, since we are doing this 

loop analysis on a non-flat space, only suhhorizon modes, modes with physical momenta 
larger than H, which are quantum modes will contribute. For these subhorizon modes one 
may approximate the □ term with its flat space value. We also note that in the region of 
parameter space we will be interested in the gauge-flation (c/. discussions of section WM 
fi"^ ^ H^f^ and H fi < f Mp\. The integral over the momenta can be cut off at a UV 
cutoff A which may be taken equal to /. Using equations in |146j . and in the MS scheme, 
we flnd 

U,ff{x) = U{x)+5U 
where sin ^ = H = -^e'"^'"' F^^F^^ and 



U{x) = -/i^ 1- Vl-H2-HarcsinH , (6.8) 

Note that to obtain the expression for U and 6U, we have evaluated for the x > 
(corresponding to the minimum of potential (16. 6p ) solution of (16. 5p . for which 

cos y = — Vl — . 

U in (16. 8 p is the minimum value of the axion potential (16. 6p . while 6U is the one-loop 
quantum corrections to this potential. The validity of one-loop analysis and integrating out 
axion around the minimum of its potential is quantifled by the 6U /U ratio; one can see that 



■^^We would like to comment that here we use the Wilsonian effective action language in which field 
renormalizations are reabsorbed into the (re)definition of fields. Field renormalizations should be considered 
if we used "renormalized perturbation theory" language |146| . 
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5U /U ~ few X /i^//'^. Therefore, if /i// < 1 one can safely approximate the potential by its 
classical value U{x)- To summarize, chromo-natural model, upon integrating out the axion 
field during the slow-roll evolution occurring around the minimum of the potential, reduces to 
a gauge theory whose action consists of the Yang-Mills term and U (S). This theory is not 
exactly our gauge-fiation model, as the potential U{E) is not just the term. Nonetheless, 
for |H| < 1, one may expand U to obtain gauge- flation K-term (c/. discussions of section 
I5.2.ip . Finally we comment that, our analysis of next subsection conforms that standard 
"slow-roll intuition" holds for our case and for slow-roll trajectories inflation is mainly driven 
by the i.e. Si/^M^ ~ ?7(S). 

In the above we focused on the field configurations and regions of parameter space for 
which the axion field is massive while the gauge field (at least some of its components, 
specifically the effective infiaton field if)) is light. One may consider another regime in 
parameter space for which the opposite situation happens. That is, gauge fields become 
massive (with effective mass of il' field to be around Hubble H) while the axions are light. 
This happens when 7 = / H"^ ^ 1. In this regime it is more appropriate to integrate out 
the gauge fields and obtain a theory which only involves axions. This has been carried out in 
|147] and leads to a natural inflation model |139] with (9x)^-type additions. Explicitly, for 
7 1, chromo-natural model below the cutoff scale A ~ /a/^/A reduces to a single scalar 
K-infiation model [H] with the action |147j 

S^ = j d'x^i^ - lidxr + ^,{dxr - /i^(l + cos(^))) + ■ • ■ , (6.10) 

here • • • denotes corrections suppressed by 7^^ and A. 



6.3 Slow-roll trajectories 

The slow-roll parameters (12. lip should remain small, e, 77 <^ 1, during the slow-roll. This 
together with (16.41) gives the following approximation for during inflation 

^1/(1 + cos ^), (6.11) 

where ~ means equality up to the first order in e. Dividing (16.4b ) and (I6.4ti ) gives e 



where 5 = —j^ and 7 = are defined as in (15.351) and (15.380 . It turns out useful to 
separate the gauge field and axion field contributions to e, z.e, 

e = e^ + e^, (6.13) 

where 
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Since both of and are positive, slow-roll condition demands that each of them should 
be very small during inflation. Also combing (12. lip and (I6.13p . one can determine t] 



V = Vx + Vth, (6.15) 
where again rj-^, rj^ are respectively x contributions to rj, 

A successful slow-roll inflation hence requires e^, e^, i]^ 1 and S < e. We will show below 
that depending on the initial value of x and xo, ^ should be a quantity of the order to e. 

We now consider the x ^'^^ equations of motion (I6.4b ..b). which in the leading slow-roll 
approximation leads to (16. 5p . One may then use (16. lip to replace H for x to find during 
the inflation 

\ 3/2 

A J ^^-^^^ 

where 

T3.^. (6.19) 

The above equation implies that 6 oc cos^. On the other hand, (I6.17P indicates that to 
have a slow-roll inflation we need S ^ e. Thus, in cases with small j, that is the "small axion 
model" discussed in [140^ 1132] . 6 should be of the order e, while for case with y = vr + 0(e), 
the case the model effectively reduces to the gauge-flation model |131j . we have 5 ~ ( Fig. 

m>- 

Within the slow-roll approximation the i/) equation of motion (16.4b ) yields 

(6.20) 

On the other hand, during the slow-roll inflation (16.4b ) ^ ^ fji'^'^^ which leads to 

ex<e^, (6.21) 

and consequently 

e ~ . (6.22) 

In other words, in the course of slow-roll inflation the energy budget of the Universe which 
is dominated by the axion potential (16. lip , is slowly, but predominantly converted into 
Yang-Mills gauge fleld energy rather than the kinetic energy of axion itselfO 



^^There are, however, regions in the parameter space where x can have a sudden change during inflation 
leading to interesting features in the dynamics of the system [142) . 
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Xo 



Figure 21: In the above we have set A^'e = 60 in the inequahty fl6.24p and plotted ^ versus A. 
The shaded area indicates parts of the A-xo plane which is excluded for successful slow-roll 
inflation. 



We can now compute number of e- folds using (16. 4b ) and (16.221) : 



AT 



1^1 + cos x)^/^(sin X 



ll/3 



^ 7x0 T2(l + COS x)4/3 + (2 sin xfl^ 



(6.23) 



where "0" subscript denotes an initial value and x = j. Before giving the complete analytic 
result for the above integral it is useful to note that, upon calling 



X = T(l+COSy)2/3 



the integrand takes the form x^+v^ ' which is always less than 1/2 and hence 



(6.24) 



We note that the above relation depends only on A and not T. Thus, setting A'^p ~ 60 in 
the inequality (16.241) . for the hilltop-type small Xo case |140] in which vr — y ~ tt, we should 
have A ~ 100. For the gauge- fiation case |121] in which n — ^ = 0(e), we need larger A 
values, A ~ 10^ |131] . Having the very simple relation (I6.24p . is very useful in restricting the 
parameter space of the model. Regardless of the values of other parameters, the shaded area 
in Figj2T] does not lead to a successful slow-roll inflation with enough number of e-folds. 
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6.3.1 Exploring the parameter space 

Up to this point, we demonstrated that the chromo- natural model can lead to slow-roll 
inflation in the entire region of Xo ^ (0,vr). Furthermore, having a preliminary relation 
between the value of Xo? ^ ^'^^ ^e, we could restrict the parameter space. Here, flrst 
working out the explicit form of the number of e- folds Nf,, we determine the parameter space 
in which successful slow-roll inflation {Ng > 60) is possible. Next, studying the system 
numerically, we show that the slow-roll analysis presented before is in excellent agreement 
with the dynamics of the system during inflation. 



The integral of number of e-folds f l6.23p may be computed analytically. To this end we 

2/' 



introduce y = sin^/'^(^), as 



Ne--TX I — ^— dy, yo = sm(^). (6.25) 



'yo 

The denominator always has a real and positive root, y 

= f--, (6.26) 

y 

since > 0, then y > 1. Furthermore, the other roots y± are given as below 



?/± = -1(1^^1-^), (6.27) 



which depending on the value of T, can be either both negative or both complex numbers: 

• T < ~ 0.73, {y > 45), y± are two real valued negative roots. Thus, the value of 
the integral is 



.r 3A 

TV ~ — 



y+ln{y-y+) y^ln{y-y^)\ 1 ^ yln{y-y) 



{y-y+) (y-y-) J {y+~y-) {y -y+){y -y-)\i 

(6.28) 



sin 1^ 



T > \/ ^ ^ 0.73, (y < 43), y^ are complex valued, y± = —y^ ± iyj, with yR = \ and 



I// = I w ^ — 1 . We then have 



.r 3A 4y2 

iV ^ 

' 2T 5^3 + 16 



\ ln((t/ + yiif + y]) + \n{y - y) - (^ + — ) arctan( ^^^^ ; 
.2 Ayi y yi 



2/ 

(6^29) 



The above expressions determine A^e as a function of T, A and Xo// which make it possible 
to explore the inflationary parameter space analytically. Note that in deriving the integral 
f l6.23p . we used the slow-roll attractor solution of ip (16. 5p . which determines it as a function 
of X- Thus, in the expression of A^'e, only xo appreases explicitly. 



129 



1000 




0.01 



50 100 



500 1000 



5000 1 X Ifft 




1000 1500 2000 3000 



5000 7000 10000 15 000 



Figure 22: In the above plots we have plotted the parameter space, A vs. T for A'^e = 60. 
The curves on the Left figure correspond to various initial values of axion field Xo] from left 
to right corresponding to y- — 10~^vr, 7r/3, 7r/2, 0.97r and O.QQtt. The Right plot shows 
the region on the parameter space which chromo- natural model is effectively equal to the 
gauge-flation, and the curves from left to right corresponds to y = O.Ott, 0.957r , O.QQtt and 
0.9927r. The "outer region of the curves correspond to N^. < 60 and are hence excluded. Note 
that A is typically (e.g. the minimum value of A) of order 10^ for the "small axion values" 
(in the left figure), while it is of order 10^ for the large axion values, in the right figure. The 



plots imply that log( 



> (log(^^^))^ while we have A^ 



lO^T 



In order for inflation to solve the flatness and horizon problems it should have lasted for 
a minimum number of e-folds which depends on the scale of inflation and somewhat to 
the details of physics after inflation ends [2]. However, it is usually demanded that N^. ^ 60. 
Using the explicit form of in fl6.28p and 06.291) . and after considering > 60, one can 
determine the accessible parameter space corresponding to each initial value of axion field, 
Xo/ f ^ (0,vr). In Figl22] setting A^'e = 60, we presented the parameter space, A versus T 
for several Xo/f values between 10~^, which is a "hilltop-type model" [140j . to 0.997r that 
is effectively reduced to "gauge-flation" [131j . The outer region of the curves correspond 
to A'^E < 60 and are hence excluded. As discussed before, A is typically of order 10^ for 
the small axion values which is the minimum value of A, while it is of order 10^ for the 
large axion values. Fig. [22] indicates that as the initial value of the axion field increases, 

A(^) T (^0 ) 

the value of A and T increase such that log( ^ (mJ — (log( "^Tk^i ))^ while we have 

A . (^) ~ lO^T . (m. 

After working out the parameter space which leads to sufficient inflation, at this point, 
we present the numerical result of studying the equations (16.41) for a set of parameters with 
different xo// values: lO'^, 7r/2, 37r/4 (Figj23]). As we see in the left panel of Fig. [231 
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Figure 23: In this figure we presented the classical trajectories with g=10~^, A=400, 
fi=7xl0~^, f = 0.009, started from different axion initial values, y. In both panels, the 
solid (orange) lines, the dashed (red) lines and the dotted (brown) lines correspond to y 
values equal to 37r/4, 7r/2 and 10~^, respectively. In the left panel, we see that axion field 
slowly increases during the inflation, then it gets equal to fir and inflation ends. As we see 
in the right panel, gradually increases during the inflation and after the inflation ends, it 
starts oscillating just like the gauge-flation model in Figs. dOlandlTTJ Note that as the initial 
value of the axion field increases, 6 gets smaller, in agreement with our slow-roll analysis. 
We also note that, as our analytical calculations show as we decrease initial value of x for 
the given {g, A, /i, /) values we get a larger number of e-folds. In particular, as we see, the set 
we have chosen gives 60 number of e-folds for Xo/f = 37r/4, while it gives about A'^g = 100 
for xo/f = 10-2. 

starting from Xo^ the axion field increases slowly during the slow-roll inflation and inflation 
ends when x/ f gets equal to vr. Furthermore, the right panel of Fig. |23l confirms that during 
the slow-roll ip is given by fl6.18p and 6 which is associated with each initial value of axion 
decreases as Xo/f gets closer to vr. Remarkably, regardless of the value of parameters and 
the initial value of fields, after slow-roll inflation ends all trajectories lead to an oscillatory 
ip, as in gauge-flation case. 

6.4 A quick treatment of chromo-natural cosmic perturbation theory 

So far, we showed that regardless of the initial axion field value, chromo-natural model can 
lead to a successful slow-roll inflation with enough number of e-folds in a vast part of its 
natural parameter space. On the other hand, almost all of the precision data of cosmology 
is related to perturbations, and we hence need to study the cosmic perturbation theory. Full 
analysis requires considering all the physical degrees of freedoms, i.e. the gauge field and 
metric perturbations as was discuss in section 15. 6[ and perturbation of the axion field 6x- 
A through and complete analysis which takes into account all the details of interactions, as 
we saw in section 15.61 has not yet appeared in the literature, see however |137l 1149] . We 
hope to discuss this in a later publication [148] . Here, we first show that vector and tensor 
perturbations in chromo-natural and gauge-flation models are exactly the same. Then, we 
review the cosmic perturbation analysis presented in |140j and |132j . 
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6.4.1 Vector and tensor perturbations in chromo-natural model 



Here we explicitly show that in vector and tensor sectors, the chromo-natural model is exactly 
the same as gauge-fiation model which has been studied in the previous section. 

The matter Lagrangian density of gauge-fiation and chromo-natural models (have been 
introduced in ( I5.2ip and (16. ip respectively) below 

= ^YM^'^Gf^ (6.30) 

'Cm = '^YAf+'^Cn' (6.31) 

where is the Yang-Mills term, = ^(F A F)^ and = A F. 

Since the axion-gauge field interaction term in the chromo-natural model is metric in- 
dependent, it it does not contribute to the energy- momentum tensor, and in this model we 

>t = (6,32) 



-9 

On the other hand, since the metric dependence of L^j, is only through detg, T'^ of the 
gauge-fiation model is given by 



n = -^ "^^~ -^.-^c., (6.33) 



which indicates that L^^ term can only contribute to the scalar sector of the perturbations 
and the vector and tensor parts are coming from the Yang-Mills term. Hence, the first order 
tensor and vector parts of the Einstein equations are exactly the same in the gauge-fiation 
and chromo-natural models. On the other hand, to have a complete set of equations, we need 
the field equation of the tensor mode perturbation of the gauge field hij which is provided 
by the tensor second order action. Similarly, enough number of equations for vector modes 
is provided by the field equation of as well as the vector parts of the Einstein equations. 
The former is obtained from the second order action of the vector perturbations. 

Considering the tensor perturbations, up to the linear order we have 5{y/—gF A F) = 0. 
Thus the contribution of and C^^. to the tensor second order action are both given by 
5'^{^—gF /\F)^ where 5^ denotes second order of perturbations. After an integration by part 
and removing the total derivative terms, contributions of these terms to the second order 
tensor action is obtained as 

2i / r.A'i 




where Xij = hij + \hij. Using (6.4b) and (5.33c), we find that they are identically the same. 
Thus, the second order action for tensor perturbations of the chromo-natural and gauge- 
flations model are exactly the same. The tensor mode perturbations of chromo-natural 
model has also appeared in |137] . 
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It is straightforward to prove, in a similar manner, that and has no contribution 
to the second order vector action, which indicates that chromo-natural and gauge-flation 
models has the same vector second order action and hence similar vector mode power spectra. 



6.4.2 Scalar modes 



Following |140l 1147] . we assume that the non-Abelian gauge field has a mass of order H 
during infiation and being very massive has a negligible effect on the quantum fluctuations. 
Therefore, the perturbations are predominantly coming from the axion field x- Oiig may then 
use the 5 A^-formulation (see Appendix [B] for a review) and after neglecting the contribution 
of the gauge field perturbations, from (1B.12P we have 



iN = ?^6x + -, (6,34) 

where N is given in fl6.25p . Then, this leads us to the following rough estimation of the 
scalar power spectrum |140] 

^'^^^ ~ ifJ' 

and scalar spectral tilt 

y" 

~ 1 - 2e + 2— . (6.36) 

In the regime where gauge fields are very massive one may also use the effective action 
obtained from integrating out the gauge fields f l6.10p and use K-infiation perturbation theory 
results 



Here for completeness we finish by summarizing two sets of possible values of parameters 
and initial fields of chromo-natural model which leads to successful model. 



Small axion model. This is the range of parameters discussed in [1401 1132] 



A = 200, ^ = 2xl0"^ e~1.8xl0~^ ^^'^^^ 
^7~1.4xlO^^ ^/^ ~ -1 X 10~^i/, V^3.3xl0^^ 

(For the above values of parameters to get enough number of e-folds, infiation should start 
when Xo// < 0.47r. 

Large axion region which is essentially the gauge-fiation region covered by parameters 

m 

Xo = vr/-5x 10-^ / = 10-^ A = 2xl0^ 

if ~ 3.3 X 10-^ /i^ = 1.6 X 10-^ g = 10~^ (6.38) 

~ 4 X 10"^ e ~ 4 X 10^^ ~ 1.6 X 10"^ 
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As we discussed in the previous subsection, however, it is possible to get slow-roll inflation 
with enough number of e-folds for all values of Xo- Note that all dimensionful quantities are 
measured in units of Mpi. 

7 Summary and outlook 

This article is devoted to reviewing various analysis and discussions appeared so far in the 
literature about the role gauge fields may have played in the dynamics of the early Universe, 
specifically in inflationary cosmology. Here we classified and discussed several such effects. 
We reviewed some models in which gauge fields are the main driving force for the inflationary 
background. Here we did not discuss the "vector inflation models" [5] , as these models involve 
vector fields and not gauge fields. It is an established common knowledge in quantum field 
theory that in the absence of gauge symmetry one should, in principle, worry about the 
consistency vector field theories and in particular absence of propagating ghosts. Presence 
of ghosts will spoil unitarity and/or positivity of the energy (Hamiltonian) of the model, and 
hence rendering the model inconsistent [6]. 

To avoid inconsistencies with the ghosts in dealing with vector fields, presence of gauge 
symmetry is essential. Gauge symmetry, among other things, forbids addition of arbitrary 
"potential" terms for the gauge fields. Moreover, one should in principle worry about the 
rotation asymmetry caused by turning vector gauge fields in the background. As we discussed 
these issues can be dealt with, if we consider non-Abelian gauge fields. This was done through 
the gauge-flation model reviewed in section [51 Here we discussed one specific inflationary 
model within the non-Abelian gauge inflation setup. The gauge invariant actions are not of 
course limited to the action we discussed; the other terms may be added to the action because 
of the loop effects and/or coupling to other fields. One such case, which is closely related to 
the gauge-flation model, was reviewed in section |6l There are many more possibilities which 
may be explored. 

One of the appealing features of the models with non-Abelian gauge fields turned on at 
the background level is the stability of the isotropic inflationary trajectories with respect to 
the choice of initial conditions as well as in the parameter space, despite the fact that we have 
vector gauge fields in the background. Providing a setup for a stable inflationary background, 
we discussed specific feature of the cosmic perturbation theory within the gauge-flation and 
chromo-natural model. It is desirable to analyze in further detail the perturbation theory of 
these models, e.g. studying non-Gaussianity and possible observable statistical anisotropic 
effects of gauge-flation or its extension. 

One of the motivations to study non-Abelian gauge field inflation was its presence in all 
high energy particle physics models and/or string theory. In fact, there is a general argument 
that all the symmetries in a theory of quantum gravity, like string theory, are gauged and 
these theories provide a wealth of gauge theories of various matter fields. The next natural 
task is hence try to embed our models in the particle physics or string theory setups. First 
steps in this direction has been taken in |133[ I15U] and a more detailed discussion is still 
needed. 
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Two other lines of research in gauge fields and infiation has been followed where gauge 
fields do not have the dominant contribution to the energy budget of the Universe during 
infiation, nonetheless, they have very interesting observational consequences. These two lines 
were discussed in sections |3] and HI As discussed quantum fiuctuations of gauge fields can, in 
principle, provide the seed for primordial magnetic fields, although no successful mechanism 
is known yet. Remarkably, gauge fields can leave observable effects on the power spectrum 
and bi-spectrum (and hence on non-Gaussianity) on the CMB anisotropies. Depending 
on various couplings that gauge fields can have with other matter fields, and in particular 
the infiaton field(s), these observable effects will have their own specific features. More 
interestingly, gauge fields with axion coupling can produce primordial gravitational waves 
even for the low scale infiation. 

The gauge fields, Abelian or non-Abelian, if they have certain kinetic coupling to infiaton 
fields can lead to anisotropic Bianchi-type quasi-de Sitter cosmology. This is interesting in 
two different ways: theoretical and observational. Theoretically, this provides examples of 
inflationary models with cosmic hair, evading the cosmic no hair conjecture [SS] and Wald's 
cosmic no hair theorem [7]. This latter, of course, happens because these models, although 
successful (quasi-de Sitter) inflationary models, do not strictly obey the assumptions of 
Wald's theorem. In fact, motivated by these examples, we proved an extension of Wald's 
theorem which all models of inflation are subject to it. This theorem is reviewed in Appendix 
iDl According to this theorem inflationary models in principle can have a cosmic hair, but 
the amount of anisotropy has an upper bound. This upper bound is stronger for slow-roll 
inflationary models and in fact models discussed in section H] saturate the bound. 

From the CMB observations viewpoint the anisotropic inflationary models have their 
own signiflcance. There are already bounds on the amount of statistical anisotropy of the 
CMB temperature fluctuations and these bounds will be improved by the upcoming CMB 
data. Hence, the anisotropy of the expansion have to be sufficiently small. However, the 
anisotropic coupling induced by vector-hair can generate the statistical anisotropy, the cross 
correlation between curvature perturbations and primordial gravitational waves, and the 
observable non-Gaussianity. As discussed in some detail in section |U the consistency relation 
among observables including non-Gaussianity can then provide a tool to discriminate between 
anisotropic inflationary models. Or conversely, if the consistency relations are conflrmed, it 
will be a strong case for the models discussed in section HJ motivating further studies of these 
models or building more inflationary models with anisotropic hair. 

Gauge fields in infiation is a subject which has received a lot of attention especially in 
recent years and as we discussed, it has many interesting theoretical, model building and 
observational aspects and we expect it to attract more attention in the coming years. This 
review contains a summary of the current results and we believe can set the stage for the 
researchers who want to work in this field. 
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A Conventions and notations 



Here we fix the conventions used in this article. 



Metric signature: We use mostly plus convention for metric, ( — h ++) signature. 

Volume form: In our conventions e^i,Q,^ is the Levi-Civita totally ant i- symmetric tensor, 
where 60123 = —a/— det g and hence e'^^^'^ = + ^^dctg " 

Natural units: Unless explicitly mentioned, we use units in which speed of light c and 
Planck constant h are set to one. In these units reduced Planck mass Mpi, SttGn = 

is 

Mpi = 2.43 X 10^^ GeV. 

It will be convenient to choose Mpi our energy/mass unit and set it to 1. That is, all 
the dimensionful quantities are measured in Planck units. 

In our conventions Einstein equation is the equation of motion for the four dimensional 
Einstein-Hilbert action: 



Seh = j (fx./d^g{]^R) , (A.l) 



where R is the Ricci scalar. 

For the flat FLRW metric (12.21) non-zero components of Ricci tensor are. 



- = ?>{H^ + H) , R, = C- + 2(^)2 + 2:^^ 5^, = {?,H^ + H + 2-^] 



and hence R = 6{2H^ + H 
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B A quick review on 5N formulation 



In this appendix, we review the 6N formahsm which is an efficient method for calculating 
power spectra and non-Gaussianity. This can be regarded as the leading order truncation 
of the long wavelength approximation which is valid on superhorizon scales. Here we will 
be very brief and for further readings we refer the reader to [16]. Let us start with FLRW 
metric in conformal time coordinate r with scalar perturbations. Recalling f l2.24p . that is 



ds^ = a2(r) [-(1 + 2A)dT^ + 2diBdx'dT + {(1 - 2C)6ij + 2dijE} dx'dx^] 
The spatial curvature can be calculated as 



(B.l) 



(B.2) 



Therefore, we can interpret C as the curvature perturbations. Using the unit normal vector 
with components Uq = — a(l + A) , rij = 0, we can compute the shear tensor as 



a. 



-adijW ag , 



where we defined the shear perturbations ag 
calculated as 



E' — B. The expansion 9 



9 = 3i!j 



A - -C 



+ - 

a 



(B.3) 



can be 



(B.4) 



On superhorizon scales, the shear perturbations ag can be negligible compared to other 
terms. Thus, using the proper time T defined by 



dT = a(l + A)dT 
we can compute the e-fold number [NT 



(B.5) 



1 
3 

N - I C'dT 





r'K - 


j OdT = ^ 






' a . 



A - 



a(l + A)dT = j pi- C] dr 



(B.6) 



The e-fold number perturbation 5N may then be defined as 

5N = y^- N{tf) = -C{rf) + Cin) , (B.7) 

where N{tf) = J. "Kdr = J. Hdt, is the background e-fold number at the final time tf. 

So far, we have not specified any time slicing. Now, let us take the flat slicing gauge C = 
on the initial hypersurface and fix the comoving gauge 50 = on the final hypersurface. 
Then, (IB.7P reads 

6N = -CiTf) . (B.8) 
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Recalling the definition of the comoving curvature perturbation 

X = C+^6<j), (B.9) 

•^c(t/) = C{Tf) on the final hypersurface, and hence 

6N = -Jlc{Tf) . (B.IO) 

Under the assumption made in the above derivation, we can read off the functional 
dependence 

l^ = y^[<i){n,x'),Tf) . (B.ll) 

In principle, we also need to take into account the ^'-dependence of 3\f. However, we can 
neglect this dependence for the slow-roll models. Thus, we finally obtain the useful relation 

dN 1 d'^N 

-'Jl^{rf) = 5N = —5ct>{n,x') + --^5ct>{n,x'f + -- - , (B.12) 

where we have used the fact that fluctuations of the e-fold number K comes from fluctuations 
of the infiaton. 



Given this formalism, we can calculate the power spectrum as 

( 

For example, for the quadratic potential V = m?(jP' /2, the e-fold number can be calculated 

as 



fdN\ ^ 

(3^,(k)3i,(kO)= hTT Mk)<50(k')) . (B.13) 



/>2 



Hence, we obtain 



N{ct>) = j^ + const. . (B.14) 



(dNV _ 02 _ 1 
where we used the slow roll parameter 



ev^ = 2^ . (B.16) 

Since we know the power spectrum of the canonically normalized scalar field at the horizon 
crossing 



{Sm^m) = (2vr)^5(k + k')^ ( ^ ) , (B.17) 
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we finally obtain 

(:R.(kWkO) = (2.mk + k')^^(Ay . (B.18) 

Similarly, we can calculate the bispectrum. In the case of the quadratic potential, we see 
that the parameter /nl can be estimated as 

/nl~t^ = ev, (B.19) 

which is the famous result obtained in [25] . 

It is straightforward to extend this formalism to multi-field cases. Even in the presence 
of other fields, as long as the expansion is isotropic, the 6N formalism would be useful. 



C Bianchi models, a quick review 

As it is standard textbook material e.g. see |1U81 1151^ 1152] . FLRW metrics are the only 
isotropic and homogeneous spacetimes. Anisotropic, but spatially homogeneous geometries 
have also been classified and known as Bianchi models0 The Bianchi family are geometries 
with spatially homogeneous (constant t) surfaces which are invariant under the action of a 
three dimensional Lie group, the symmetry group and can hence be foliated into the spatial 
homogeneous hypersurfaces St, 

M = M X St, 

here M is the time variable. In each of these spatially homogeneous surfaces there exist a set 
of basis vectors Cj (Killing vectors) that spans a Lie algebra as 

[ei,e,] = C%ek. (C.l) 

A Bianchi metric can always be written as (e.g. see |153j ) 

ds^ = -dt^ + hij(t)e' ® e\ (C.2) 

where z,j = 1,2,3 label the coordinates in homogeneous space-like hypersurfaces, hij is the 
homogeneous spatial metric and e* are one-forms dual to basis Cj. Recalling ( IC.ip . we find 
the following property for e*s 

de^ = --C^.^e-^' Ae^ (C.3) 



Moreover, the spatial scalar curvature '^'i? is given as 



'''r = -C\f\^ + \c^^,C\ ' - -C,,,C^=\ (C.4) 



■^^The Kantowski-Sachs model is also a spatially homogeneous model which has a four dimensional sym- 
metry group [151] However, since this solution is unstable we do not consider it here. 
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which indicates that each three dimensional Lie algebra corresponds to a spatially homoge- 
neous cosmological model. 

Bianchi classification categorized the family of three dimensional Lie algebras into 9 
different classes in which each algebra is labeled by a number I-IX. For instance, Bianchi 
type /, which generalized the flat FLRW metric with the symmetry group described by 
C\j = 0, corresponds to flat hypersurfaces. Meanwhfle, Bianch type IX corresponds to the 

so(3) Lie algebra with a positive '^'i?. In fact all Bianchi models, expect type IX, have a 
non-positive spatial scalar curvature [TJ 1153] 

^'^R < 0. (C.5) 

It is common to write the structure constants C\j in terms of a trace part and a trace-free 
part, the so cafled Behr decomposition [154j . 

C% = A[,5^^+e,,kM''\ (C.6) 

where M^'' is a symmetric tensor. The models with C^'^j = are called class A models, whiles 
modes with ^ are called class B models. 



D Cosmic no-hair theorem and its extensions 

By the current data the Universe as cosmological scales looks homogeneous and isotropic up 
to one in 10^ part. On the other hand, one of the theoretically attractive feature of inflation- 
ary cosmology has been detaching, to a good extent, the outcome of inflation from the initial 
conditions. These initial conditions, unless severely flne-tuned, are not homogeneous and 
isotropic. Therefore, it is very appealing if one can seek a "dynamical" explanation for this 
homogeneity and isotropy. That is, to show that isotropic and homogeneous Universe is an 
attractor of the cosmic evolution. Within the standard model of cosmology cosmic evolution 
is governed by the Einstein gravity coupled to the cosmic fluid, which is not necessarily a 
perfect fluid and hence if this latter idea is true, homogeneity and isotropy should be an 
outcome of Einstein equations. The flrst such attempt, which was dubbed as cosmic no-hair 
conjecture was made in [89] arguing that the late-time behavior of any accelerating Universe 
is an isotropic Universe. The flrst attempt to prove this conjecture was presented in Wald's 
seminal paper [7|, which is called (Wald's) cosmic no-hair theorem. 

Wald's cosmic no-hair theorem [7] states that Bianchi-type models (except Bianchi 
IX) with the total energy momentum tensor of the form 

and with a constant positive Aq and a T^y satisfying Strong and Dominant Energy Conditions, 
respectively SEC and DEC, approach de Sitter space exponentially fast, within a few Hubble 
times = ^S/Aq. The fate of Bianchi type IX is similar, if Aq is large enough compared 
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with spatial curvature terms. As a result, in these systems inflation never ends and Universe 
will appear to be an empty de Sitter space. For a more thorough historical review and other 
related works on cosmic no-hair conjecture/theorem see |152] and references therein. 

As stated and formulated, Wald's cosmic no-hair theorem is not primarily meant for 
inflationary models and in fact inflationary models better not obey the fate of Universe 
predicted by Wald's theorem, because to match the observations, inflation must end and 
we should enter a phase which is not (quasi) de Sitter at the end of inflation. This is of 
course no surprise, because inflationary models, by construct, do not satisfy the DEC and 
SEC conditions of Wald's theorem (to have a built-in mechanism for terminating inflation). 
So, Wald's theorem, per se, is not applicable to inflationary models. In [8], and motivated 
by the anisotropic inflationary models discussed in detail in section HI Wald's cosmic no-hair 
theorem was extended to include inflationary models too. This was done by modifying SEC 
and DEC to milder conditions respected by requirement of having accelerated expansion. 



Inflationary extended cosmic no-hair theorem [8] states that for general inflationary 
systems of all Bianchi type with energy momentum tensor of the form 

where A(t) is a cosmological term which decreases by time and 7^^ satisfies SEC and WEC, 
in contrast to the cosmic no-hair conjecture (and also Wald's cosmic no-hair theorem), 
anisotropy may grow nonetheless there is an upper bound on the growth of anisotropy. 
For slow-roll inflationary models this upper bound is of the order of slow-roll parameters. 

In what follows, we first present the general setup in which Wald's cosmic no-hair theorem 
and its inflationary extension can be formulated and then present the two theorems. 

D.l Cosmic no-hair theorems, a general setup 

Consider a general Bianchi model and let be the unit tangent vector field of the congruence 
of time- like geodesies orthogonal to the homogeneous space-like hypersurfaces St. Then, we 
obtain the following covariant form for the spatial metric h^^ (lC.2p 

hf,u = gf,u + n^,n^, (D.l) 

where g^y is the metric of the spacetime. Moreover, hij = hij(t) can be decomposed as 

/i^. = e2V^'^ (D.2) 

where e'^ is the isotropic scale factor and is a traceless matrix which describes the 
anisotropy. 

Extrinsic curvature of is defined as 
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where the dot represents derivative with respect to the time coordinate t. One can decompose 
Kfj^i, into trace, and traceless parts 

where 

K{t) = K^.h''-' = 3H{t) . (D.5) 

Here H{t) = d is the Hubble parameter corresponding to the homogeneous scale factor e° in 
flD.2l) . The shear of the time-like geodesic congruences is related to the time derivative 
of Pij and is a symmetric, traceless and purely spatial tensor: 



We now analyze Einstein's equation for Bianchi models which can be decomposed into 
four constraint equations 



T^.n'^n'' = ^^'R--a^,a^^ + -K\ (D.7) 

T.,/iV = K'^fl. + K'^fi'.,. (D.8) 
and six dynamical equations 

(T^A - ]^KxT)h'',h\ = 2naij + ^{Kh,, + K^hij + Ka,j) - 2a,xa) + . (D.9) 

(3) 

Here Rij is the spatial Ricci tensor and can be written in terms of the structure-constant 
tensor C*^ (1C.3P as 

(3)j- _ _ 



where i,j indices are raised and lowered with metric hij, and '^'i? = ^^^Rij0 is the spatial 
curvature of St ( 1C.4I) . As mentioned in the Appendix [Cl all but one (Bianchi IX) of Bianchi 
models have negative spatial curvature. 

Combining trace of f ID.Qp with (]D.7p . we obtain Raychaudhuri equation 

{T^u - ^g,,T)n^n'' = -K - - G^,a^\ (D.ll) 

where K = 2inK, and contracting ( ]D.9P with h^^ and removing the trace part, we obtain the 
equation for the shear tensor: 



a'^ + Ka'^ + = Tkih^'h'^ - ^Tuh^'h' ^ , (D.12) 



where '^^5** • is the anisotropic part of the spatial 3-curvature 



^''^^ . = ''^i?^ . - 1 '-'W^. (D.13) 
o 
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Integrating equation (1D.12|) . we obtain the following integral equation for cf^ At) 



^^^^{Tkih^'h^, - f^kih^'h^, - + c^e-^'^W , (D.14) 

with integration constants c*^. In order to determine cr^jit) we need more information about 
the energy momentum tensor T^y. Before that, let us remind three common energy conditions 
which are defined as: 

• Strong energy condition (SEC) is satisfied when for energy momentum tensor T^^ 
and for all time-like t'^ 

(V - l9>.u7)t''t'' > 0. (D.15) 

• Dominant energy condition (DEC) stipulates that for all future-directed causal 
vectors t'',t"' 

Tt.Xt"' > 0. (D.16) 

• Weak energy condition (WEC) is defined through the dominant energy condition 
for t'^' = That is, 

l^.f^t" > 0. (D.17) 

At this point, we consider two different forms for T^,^, one leads to the Wald's cosmic 
no-hair theorem and the other one is the most general form of the energy momentum tensor 
for an inflationary setting and investigate the dynamical evaluation of anisotropy for each 
case. 



D.2 Wald's cosmic no-hair theorem 

Following [7], here we consider initially expanding general Bianchi models with the following 
energy momentum tensor 

in which Aq is a positive cosmological constant and T^^ satisfies SEC and DEC. Putting 
(|DA8|) into ^Dl\\ and flmTj) . we have 



= 3Ao + ^a^X--^ '''i^ + 3Vn'^n^ (D.19) 
K = Ao-^K^ -a^''a^,-{7^,-^g^,7)n^n\ (D.20) 



After using (ID.lSp . f lD.17p and fID.lOp . from the combination of (ID.19P and f lD.20p we obtain 



the following inequalities for all Bianchi types except IX 



K'^ > 3Ao and K <A- -K^ < 0. (D.21^ 

3 
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Then integrating the inequahty as well as using the lower bound on gives 



(3Ao)^<i^< ^^^'}\ . (D.22) 
tanh((f 

As we see, K is sandwiched between its lower bound (SAq)^ and an upper bound which 

exponentially approaches (3Ao)2 with a time scale Applying this result on (1D.19P 
again, we obtain 

< , and T^,n^'n'' < ^\ , . (D.23) 

' sm\i\{^)n) ^ " sinh2((f ^ ' 

Thus, the shear of the homogeneous hypersurfaces as well as the energy density of T^j, rapidly 
approach zero. Since satisfies DEC, 7^^n^n'^ is larger or equal than the rest of orthonor- 
mal frame elements of 7^^. Hence, the homogeneous hypersurfaces and all components of 
7^^ exponentially approach zero, as such the geometry, regardless of the initial conditions, 
evolve toward the de Sitter solution. 



D.3 Inflationary extended cosmic no-hair theorem 

Following [Hj, here we consider inflationary Bianchi models with the following energy mo- 
mentum tensor 

T^, = -Ait)g^, + V, (D.24) 

in which A(t) is a positive but decreasing in time cosmological term and T^,y satisfies SEC 
and WEC. Energy momentum tensors of inflationary models fall into the class we consider 
here. 



Without loss of generality T^,^ in (ID.24P may be decomposed as 



= p{t)n^n, + P(t) V + ^f^At), (D-25) 
where tt^i, is the anisotropic stress tensor which is symmetric, traceless, and purely spatial: 

TT^^/i/^'^ = & vr^^n'^ = . (D.26) 

We also define Aq = 3Hq, where subscript denotes the initial value (at the beginning of 
inflation) and 

A(t) = Ao + 6A{t) . (D.27) 

Moreover, the spatial curvature '^'i? has a time dependence proportional to e~^" and this 
term is damped quickly in inflationary systems. Therefore, using (lD.24p and (ID.25P in (ID.7P 
and fID.lip . after a few e- folds '^'i? is negligible and we approximately have 



6 Ait) ~p{t)-P{t) e{t) , H, 



(D.28) 



3i/2(t) QH^(t) 3 H^{t) 

e{t). (D.29) 



P{t) + P{t) , a^^a^^'' 



2H^{t) 2H\t) 
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On the other hand, T^i, satisfies SEC which choosing = n^+s^, where is a normahzed 
arbitrary space-hke 4- vector orthogonal to n^, imphes 

p{t) + P{t) +7rij{t)s's^ >0, i,j = 1,2,3, (D.30) 

where s* = s^h^^. Recalhng the fact that vTjj is traceless and the above inequahty should 



hold for all s*, one obtain: 
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H\t) - H^t) 



From the combination of fiDiM and (10129]) and noting that (iDlSOl) implies p{t) + P{t) > 
we obtain 



a^,a^' < 2H\t)e{t) 



M(t)<-(p(t) + ^^)<0, 



(D.32) 



The above already implies a weak upper bound on anisotropy 

\a'j\ < H{t)V2^. (D.33) 

As we will see later, assuming slow-roll leads to a stronger bound. 

During inflation (on the average) e is an increasing quantity and this result among other 
things shows the possibility of growth for the anisotropic part of the stress tensor. This is 
in contrast with behavior of systems described by the Wald's theorem, in which all elements 
of 7fj,u, including 7r*j(t)'s, are damped exponentially with a time scale Hq^ [7j. 

Inserting (ID. 251) into (ID. 141) . we have the following form for the Hubble- normalized shear 
tensor 



H{t) 3Hit)eMt)J,^ (1 + j^) ^'(^')' 

where to is a few e-folds after the beginning of inflation in which the spatial curvature terms 

has become negligible. As we see 'jjjjjj acts as a source term for anisotropy cr*^ and its flve 
degrees of freedom may be determined for speciflc models of inflation. 

For models with perfect fluid T^,y, in which vr*j(t) is identically zero, shear tensor cr^jit) 
is damped exponentially fast. Thus, in this kind of systems inflation washes away any initial 
anisotropy and isotropizes the system. These models hence comply with the cosmic no-hair 
conjuncture. However, this is not necessarily the case in a general inflationary model with 
non-zero anisotropic stress tensor vr*j(t). In fact, in presence of anisotropic stress tensor vr*^- 
elements of shear tensor a^j (t) are allowed to grow by inflationary dynamics, which is against 
the cosmic no-hair conjuncture. Although j{t) can increase, as one expects from ( ID. 311) . 
inflation enforces an upper bound on how large they can grow. 



•^^ Although we have not used here, we remind that WEC on T^jy imply p > 0. 
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Quasi-de Sitter expansion 



In the following we focus on the slow-roll models. We investigate the dynamics of system 
assuming that 7r^j{t)/H'^{t) always saturates its maximum value, and determine the possible 
upper bound value of anisotropies during slow-roll inflation. 

From (1D.32P and (lD.3ip . we find that 7^^, cr^^cr^^'^ and '^'i? are at most of the order eH'^. 
Hence recalling flD.27|) and using fID.Tp . we obtain 



H{t) ~ Ho{l-eoHot), (D.35) 

where ~ means to the first order in slow-roll parameters, subscript denotes the initial value 
and Hq = ^ Ao/3. Inserting the above relation in flD.28P and flD.29p . we obtain 

^W^-ifi - .„(1.2(.„-,„W). (D.3T) 

As we see in (lD.37p . o^^Op^y is not directly related to the slow-roll dynamics, but the com- 
bination of ^^2H^t)^ ~^ "iH^lt) should be slow varying during inflation. As a result, (jf^y is not 
fully determined by the slow-roll dynamics, while ( ID. 361) implies that the isotropic part of 
the system is governed by the slow-roll. Thus, assuming slow-roll inflation, anisotropies can 
still evolve quickly in time. 



Applying slow-roll approximation in f lD.34p . we obtain the following form for the late 



time behavior of the Hubble-normalized shear tensor 



H{t) 



^ - ifo(l + e,Hot)e-'^^' f dt'e'"^^^'{l - 2eM)^^. (D.38) 

J to -n ) 



Hit) 



In order to determine 7^^, we need 7r*j(t)'s five degrees of freedom which should be provided 
by independent equations. However, regardless of the particular inflationary model which 



we consider, recalling ( ]D.3ip . slow- roll enforces an upper bound on the anisotropy growth. 



Rewriting (ID. 310 in terms of slow-roll parameters, we obtain 



< 4e(t) - 4eo(l + 2(eo - Vo)Hot)) . (D.39) 



As mentioned before, from the combination of (ID. 291) and (ID. 341) . we realize that during 
the slow-roll inflation, 'ir^j{t)/H'^{t) can saturate its upper bound and grow in time. Then, 
assuming that TT^j{t)/H'^{t) always saturates its upper bound during the slow- roll, we can 
find an upper bound value on the value of At)/ H{t) enforced by the slow-roll dynamics. 



''^Note that we are considering simple slow-roll models where e(t) is an always increasing function during 
inflation, e > 0. From the definition of e (|2.11[) we learn that ~ 2(e — 77). 
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Inserting (1D.39|) into ( ]D.38|) . we obtain the upper bound value on the anisotropies at the 

<?(eo-^o). (D.40) 



end of slow-roll inflation 



where tsi is the end of slow-roll inflation and approximately is Hoeotgi — 1. Thus, in a general 
model of slow-roll inflation we end up with an almost isotropic Universe, with a shear of the 
order eo- Note that slow-roll dynamics has improved the bound (1D.33P to (ID. 400 . 

To summarize, we have proved an inflationary extended no-hair theorem, according which 
despite the fact that anisotropies can in principle grow during inflation (when the anisotropic 
part of energy momentum tensor of the inflaton sector is non-zerc@), having quasi-de Sit- 
ter expansion puts an upper bound on anisotropy: Hubble-normalized anisotropic stress 
71^ j{t)/ H'^{t), proportional to the slow-roll parameter e{t) and consequently an upper on the 
dimensionless anisotropies aj{t)/H(t) < 2y/e. Slow-roll assumption improves this bound to 

— ^ eo- 

Finally, we end this Appendix with some comments: 

1) Note that e is generically a growing function as inflation progresses, and hence in principle, 
anisotropy can grow, signaling the possibility of violation of cosmic no-hair conjecture. In 
fact, as reviewed in section m explicit models of inflation have been constructed exhibiting 
this [9l]. All the models discussed in section H] nicely comply with our inflationary extended 
no-hair theorem. 

2) Even assuming slow-roll inflation, the dynamics of anisotropy does not necessarily follow 
the slow-roll dynamics and they can generally evolve quickly, and 

3) our analysis is general and our upper bounds is expected to hold for all Bianchi models, 
including type IX. 
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